Spatial and spectral wavefront analysis and measurement

ABSTRACT

A method and apparatus for wavefront analysis including obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed which has an amplitude and a phase, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed.

REFERENCE TO RELATED APPLICATIONS

[0001] This application is based on provisional application U.S. Ser. No. 60/196,862, filed on Apr. 12, 2000.

FIELD OF THE INVENTION

[0002] The present invention relates to wavefront analysis generally and to various applications of wavefront analysis.

BACKGROUND OF THE INVENTION

[0003] The following patents and publications are believed to represent the current state of the art:

[0004] U.S. PATENTS:

[0005] U.S. Pat. Nos. 5,969,855; 5,969,853; 5,936,253; 5,870,191; 5,814,815; 5,751,475; 5,619,372; 5,600,440; 5,471,303; 5,446,540; 5,235,587; 4,407,569; 4,190,366;

[0006] NON-US PATENTS:

[0007] JP 9230247 (Abstract); JP 9179029 (Abstract); JP 8094936 (Abstract); JP 7261089 (Abstract); JP 7225341 (Abstract); JP 6186504 (Abstract);

[0008] OTHER PUBLICATIONS:

[0009] Phillion D. W. “General methods for generating phase-shifting interferometry algorithms”—Applied Optics, Vol. 36, 8098 (1997).

[0010] Pluta M. “Stray-light problem in phase contrast microscopy or toward highly sensitive phase contrast devices: a review”—Optical Engineering, Vol. 32, 3199 (1993).

[0011] Noda T. and Kawata S. “Separation of phase and absorption images in phase-contrast microscopy”—Journal of the Optical Society of America A, Vol. 9., 924 (1992).

[0012] Creath K. “Phase measurement interferometry techniques”—Progress in Optics XXVI, 348 (1988).

[0013] Greivenkamp J. E. “Generalized data reduction for heterodyne interferometry”—Optical Engineering, Vol. 23, 350 (1984).

[0014] Morgan C. J. “Least-squares estimation in phase-measurement interferometry”—Optics Letters, Vol. 7, 368 (1982).

[0015] Golden L. J. “Zernike test. 1: Analytical aspects”—Applied Optics, Vol. 16, 205 (1977).

[0016] Bruning J. H. et al. “Digital wavefront measuring interferometer for testing optical surfaces and lenses”—Applied Optics, Vol. 13, 2693 (1974).

SUMMARY OF THE INVENTION

[0017] The present invention seeks to provide methodologies and systems for wavefront analysis as well as for surface mapping, phase change analysis, spectral analysis, object inspection, stored data retrieval, three-dimensional; imaging and other suitable applications utilizing wavefront analysis.

[0018] There is thus provided in accordance with a preferred embodiment of the present invention a method of wavefront analysis. The method includes obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed which has an amplitude and a phase, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed.

[0019] There is also provided in accordance with a preferred embodiment of the present an apparatus for wavefront analysis including a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed which has an amplitude and a phase, an intensity map generator operating to provide a plurality of intensity maps of the plurality of phase changed transformed wavefronts and an intensity map utilizer, employing the plurality of intensity maps for providing an output indicating the amplitude and phase of the wavefront being analyzed.

[0020] Further in accordance with a preferred embodiment of the present invention the plurality of intensity maps are employed to provide an analytical output indicating the amplitude and phase.

[0021] Still further in accordance with a preferred embodiment of the present invention the plurality of differently phase changed transformed wavefronts are obtained by interference of the wavefront being analyzed along a common optical path.

[0022] Additionally in accordance with a preferred embodiment of the present invention the plurality of differently phase changed transformed wavefronts are realized in a manner substantially different from performing a delta-function phase change to the transformed wavefront.

[0023] Further in accordance with a preferred embodiment of the present invention the plurality of intensity maps are employed to obtain an output indicating the phase which is substantially free from halo and shading off distortions.

[0024] Preferably, the plurality of differently phase changed transformed wavefronts include a plurality of wavefronts resulting from at least one of application of spatial phase changes to a transformed wavefront and transforming of a wavefront following application of spatial phase changes thereto.

[0025] Additionally in accordance with a preferred embodiment of the present invention, the step of obtaining a plurality of differently phase changed transformed wavefronts includes applying a transform to the wavefront being analyzed thereby to obtain a transformed wavefront and applying a plurality of different phase changes to the transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts. Preferably, the plurality of different phase changes includes spatial phase changes and the plurality of different spatial phase changes are effected by applying a time-varying spatial phase change to part of the transformed wavefront.

[0026] Further in accordance with a preferred embodiment of the present invention the plurality of different spatial phase changes are effected by applying a spatially uniform, time-varying spatial phase change to part of the transformed wavefront. Preferably, the transform applied to the wavefront being analyzed is a Fourier transform and wherein the step of obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts includes applying a Fourier transform to the plurality of differently phase changed transformed wavefronts.

[0027] Further in accordance with a preferred embodiment of the present invention the transform applied to the wavefront being analyzed is a Fourier transform and the plurality of different spatial phase changes includes at least three different phase changes. Preferably, the plurality of intensity maps includes at least three intensity maps and the step of employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed includes: expressing the wavefront being analyzed as a first complex function which has an amplitude and phase identical to the amplitude and phase of the wavefront being analyzed, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change, defining a second complex function, having an absolute value and a phase, as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change. Expressing each of the plurality of intensity maps as a third function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function, a difference between the phase of the wavefront being analyzed and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes corresponding to one of the at least three intensity maps, solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the wavefront being analyzed by adding the phase of the second complex function to the difference between the phase of the wavefront being analyzed and the phase of the second complex function.

[0028] Further in accordance with a preferred embodiment of the present invention the absolute value of the second complex function is obtained by approximating the absolute value to a polynomial of a given degree.

[0029] Still further in accordance with a preferred embodiment of the present invention the second complex function is obtained by expressing the second complex function as an eigen-value problem where the complex function is an eigen-vector obtained by an iterative process.

[0030] Preferably the second complex function is obtained by: approximating the Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change to a polynomial and approximating the second complex function to a polynomial.

[0031] Preferably, the wavefront being analyzed, the absolute value of the second complex function, and the difference between the phase of the second complex function and the phase of the wavefront being analyzed, are obtained by a least-square method, which has increased accuracy as the number of the plurality of intensity maps increases.

[0032] Further in accordance with a preferred embodiment of the present invention the plurality of different phase changes includes at least four different phase changes, the plurality of intensity maps includes at least four intensity maps and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed and includes: expressing each of the plurality of intensity maps as a third function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function, a difference between the phase of the wavefront being analyzed and the phase of the second complex function, a known phase delay produced by one of the at least four different phase changes in which each corresponds to one of the at least four intensity maps and at least one additional unknown relating to the wavefront analysis, where the number of the at least one additional unknown is no greater than the number by which the plurality intensity maps exceeds three and solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function, the difference between the phase of the wavefront being analyzed and the phase of the second complex function and the additional unknown.

[0033] Further in accordance with a preferred embodiment of the present invention the phase changes are chosen as to maximize contrast in the intensity maps and to minimize effects of noise on the phase of the wavefront being analyzed.

[0034] Preferably, expressing each of the plurality of intensity maps as a third function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function, a difference between the phase of the wavefront being analyzed and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes which corresponds to one of the at least three intensity maps and includes: defining fourth, fifth and sixth complex functions, none of which being a function of any of the plurality of intensity maps or of the time-varying spatial phase change, each of the fourth, fifth and sixth complex functions being a function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function and expressing each of the plurality of intensity maps as a sum of the fourth complex function, the fifth complex function multiplied by the sine of the known phase delay corresponding to each one of the plurality of intensity maps and the sixth complex function multiplied by the cosine of the known phase delay corresponding to each one of the plurality of intensity maps.

[0035] Preferably, the step of solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function includes: obtaining two solutions for each of the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function, the two solutions being a higher value solution and a lower value solution, combining the two solutions into an enhanced absolute value solution for the absolute value of the second complex function, by choosing at each spatial location either the higher value solution or the lower value solution of the two solutions in a way that the enhanced absolute value solution satisfies the second complex function. Preferably, combining the two solutions of the amplitude of the wavefront being analyzed into enhanced amplitude solution, by choosing at each spatial location the higher value solution or the lower value solution of the two solutions of the amplitude in the way that at each location where the higher value solution is chosen for the absolute value solution, the higher value solution is chosen for the amplitude solution and at each location where the lower value solution is chosen for the absolute value solution, the lower value solution is chosen for the amplitude solution, combining the two solutions of the difference between the phase of the wavefront being analyzed and the phase of the second complex function into an enhanced difference solution, by choosing at each spatial location the higher value solution or the lower value solution of the two solutions of the difference in the way that at each location where the higher value solution is chosen for the absolute value solution, the higher value solution is chosen for the difference solution and at each location where the lower value solution is chosen for the absolute value solution, the lower value solution is chosen for the difference solution.

[0036] Further in accordance with a preferred embodiment of the present invention the spatially uniform, time-varying spatial phase change is applied to a spatially central part of the transformed wavefront.

[0037] Preferably, the transform applied to the wavefront being analyzed is a Fourier transform and wherein the step of obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts includes applying a Fourier transform to the plurality of differently phase changed transformed wavefronts.

[0038] Still further in accordance with a preferred embodiment of the present invention the method also includes adding a phase component including relatively high frequency components to the wavefront being analyzed prior to applying the transform thereto in order to increase the high-frequency content of the transformed wavefront prior to the applying the spatially uniform, time-varying spatial phase change to part of the transformed wavefront.

[0039] Preferably, the spatially uniform, time-varying spatial phase change is applied to a spatially centered generally circular region of the transformed wavefront and the spatially uniform, time-varying spatial phase change is applied to approximately one half of the transformed wavefront.

[0040] Additionally in accordance with a preferred embodiment of the present invention the transformed wavefront includes a DC region and a non-DC region and the spatially uniform, time-varying spatial phase change is applied to at least part of both the DC region and the non-DC region.

[0041] Further in accordance with a preferred embodiment of the present invention the plurality of differently phase changed transformed wavefronts include a plurality of wavefronts whose phase has been changed by employing an at least time varying phase change function. Alternatively, the plurality of differently phase changed transformed wavefronts include a plurality of wavefronts whose phase has been changed by applying an at least time varying phase change function to the wavefront being analyzed.

[0042] Preferably, the at least time varying phase change function is applied to the wavefront being analyzed prior to transforming thereof. Alternatively, the at least time varying phase change function is applied to the wavefront being analyzed subsequent to transforming thereof.

[0043] Further in accordance with a preferred embodiment of the present invention the plurality of differently phase changed transformed wavefronts include a plurality of wavefronts whose phase has been changed by employing an at least time varying phase change function.

[0044] Additionally or alternatively, the plurality of differently phase changed transformed wavefronts include a plurality of wavefronts whose phase has been changed by applying an at least time varying phase change function to the wavefront to be analyzed.

[0045] Preferably, the at least time varying phase change function is a spatially uniform spatial function.

[0046] Additionally in accordance with a preferred embodiment of the present invention the transformed wavefront includes a plurality of different wavelength components and the plurality of different spatial phase changes are effected by applying a phase change to the plurality of different wavelength components of the transformed wavefront.

[0047] Preferably, the phase change applied to the plurality of different wavelength components of the transformed wavefront is a time-varying spatial phase change.

[0048] Further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components of the transformed wavefront is effected by passing the transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.

[0049] Still further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components of the transformed wavefront is effected by reflecting the transformed wavefront from a spatially varying surface.

[0050] Further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components of the transformed wavefront is selected to be different to a predetermined extent for at least some of the plurality of different wavelength components.

[0051] Additionally in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components of the transformed wavefront is identical for at least some of the plurality of different wavelength components.

[0052] Further in accordance with a preferred embodiment of the present invention the wavefront being analyzed includes a plurality of different wavelength components. Preferably, the plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to the plurality of different wavelength components of the wavefront being analyzed.

[0053] Preferably, the phase change is applied to the plurality of different wavelength components of the wavefront being analyzed prior to transforming thereof.

[0054] Further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components is effected by passing the wavefront being analyzed through an object, at least one of whose thickness and refractive index varies spatially.

[0055] Further in accordance with a preferred embodiment of the present invention the step of obtaining a plurality of intensity maps is performed simultaneously for all of the plurality of different wavelength components and obtaining a plurality of intensity maps includes dividing the plurality of phase changed transformed wavefronts into separate wavelength components.

[0056] Still further in accordance with a preferred embodiment of the present invention the step of dividing the plurality of phase changed transformed wavefronts is effected by passing the plurality of phase changed transformed wavefronts through a dispersion element.

[0057] Additionally in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components is effected by passing the wavefront being analyzed through an object, at least one of whose thickness and refractive index varies spatially, following transforming of the wavefront being analyzed.

[0058] Preferably, the phase change which is applied to the plurality of different wavelength components is effected by reflecting the wavefront being analyzed from a spatially varying surface, following transforming of the wavefront being analyzed.

[0059] Further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components is selected to be different to a predetermined extent for at least some of the plurality of different wavelength components. Preferably, the phase change which is applied to the plurality of different wavelength components is identical for at least some of the plurality of different wavelength components.

[0060] Further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different wavelength components is effected by passing the wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.

[0061] Preferably, the phase change applied to the plurality of different wavelength components is effected by passing the wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially, following transforming of the wavefront being analyzed.

[0062] Further in accordance with a preferred embodiment of the present invention the wavefront being analyzed includes a plurality of different polarization components and the plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to the plurality of different polarization components of the wavefront being analyzed prior to transforming thereof.

[0063] Still further in accordance with a preferred embodiment of the present invention the transformed wavefront includes a plurality of different polarization components and the plurality of different spatial phase changes are effected by applying a phase change to the plurality of different polarization components of the transformed wavefront.

[0064] Additionally in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different polarization components of the transformed wavefront is different for at least some of the plurality of different polarization components.

[0065] Further in accordance with a preferred embodiment of the present invention the phase change applied to the plurality of different polarization components of the transformed wavefront is identical for at least some of the plurality of different polarization components.

[0066] Additionally in accordance with a preferred embodiment of the present invention the step of obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts includes applying a transform to the plurality of differently phase changed transformed wavefronts.

[0067] Preferably, the plurality of phase changed transformed wavefronts are reflected from a reflecting surface so that the transform applied to the plurality of differently phase changed transformed wavefronts is identical to the transform applied to the wavefront being analyzed.

[0068] Further in accordance with a preferred embodiment of the present invention the transform applied to the wavefront being analyzed is a Fourier transform.

[0069] Still further in accordance with a preferred embodiment of the present invention the plurality of intensity maps are obtained by reflecting the plurality of differently phase changed transformed wavefronts from a reflecting surface so as to transform the plurality of differently phase changed transformed wavefronts.

[0070] Additionally in accordance with a preferred embodiment of the present invention the method of obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts includes applying a transform to the plurality of differently phase changed transformed wavefronts.

[0071] Further in accordance with a preferred embodiment of the present invention the method of employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed includes expressing the plurality of intensity maps as at least one mathematical function of phase and amplitude of the wavefront being analyzed and employing the at least one mathematical function to obtain an output indicating the phase and amplitude.

[0072] Preferably, the method of employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed includes expressing the plurality of intensity maps as at least one mathematical function of phase and amplitude of the wavefront being analyzed and of the plurality of different phase changes, wherein the phase and amplitude are unknowns and the plurality of different phase changes are known and employing the at least one mathematical function to obtain an output indicating the phase and amplitude.

[0073] Further in accordance with a preferred embodiment of the present invention the plurality of intensity maps includes at least four intensity maps and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed, includes employing a plurality of combinations, each of at least three of the plurality of intensity maps, to provide a plurality of indications of the amplitude and phase of the wavefront being analyzed.

[0074] Preferably, the method also includes employing the plurality of indications of the amplitude and phase of the wavefront being analyzed to provide an enhanced indication of the amplitude and phase of the wavefront being analyzed.

[0075] Further in accordance with a preferred embodiment of the present invention at least some of the plurality of indications of the amplitude and phase are at least second order indications of the amplitude and phase of the wavefront being analyzed.

[0076] Further in accordance with a preferred embodiment of the present invention the step of obtaining a plurality of differently phase changed transformed wavefronts includes applying a transform to the wavefront being analyzed, thereby obtaining a transformed wavefront and applying a plurality of different phase and amplitude changes to the transformed wavefront, thereby obtaining a plurality of differently phase and amplitude changed transformed wavefronts.

[0077] Further in accordance with a preferred embodiment of the present invention the plurality of different phase and amplitude changes includes at least three different phase and intensity changes, the plurality of different phase and amplitude changes are effected by applying at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial amplitude change to at least part of the transformed wavefront, the plurality of intensity maps includes at least three intensity maps. Preferably, the step of employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed includes: expressing the wavefront being analyzed as a first complex function which has an amplitude and phase identical to the amplitude and phase of the wavefront being analyzed, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial amplitude change, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change, expressing each of the plurality of intensity maps as a third function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function and a difference between the phase of the wavefront being analyzed and the phase of the second complex function. Preferably, the spatial function governing at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial amplitude change includes: defining fourth, fifth, sixth and seventh complex functions, none of which is a function of any of the plurality of intensity maps or of the time-varying spatial phase change. Preferably, each of the fourth, fifth, sixth and seventh complex functions being a function of at least one of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function, defining an eighth function of a phase delay and of an amplitude change, both produced by one of the at least three different phase and amplitude changes, corresponding to the at least three intensity maps and expressing each of the plurality of intensity maps as a sum of the fourth complex function, the fifth complex function multiplied by the absolute value squared of the eighth function, the sixth complex function multiplied by the eighth function and the seventh complex function multiplied by the complex conjugate of the eighth function, solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the wavefront being analyzed by adding the phase of the second complex function to the difference between the phase of the wavefront being analyzed and phase of the second complex function.

[0078] Further in accordance with a preferred embodiment of the present invention the wavefront being analyzed includes at least two wavelength components. Preferably, the step of obtaining a plurality of intensity maps also includes dividing the phase changed transformed wavefronts according to the at least two wavelength components in order to obtain at least two wavelength components of the phase changed transformed wavefronts and in order to obtain at least two sets of intensity maps, each set corresponding to a different one of the at least two wavelength components of the phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed, obtaining an output indicative of the phase of the wavefront being analyzed from each of the at least two sets of intensity maps and combining the outputs to provide an enhanced indication of phase of the wavefront being analyzed, in which enhanced indication, there is no 2π ambiguity.

[0079] Additionally in accordance with a preferred embodiment of the present invention the wavefront being analyzed is an acoustic radiation wavefront.

[0080] Still further in accordance with a preferred embodiment of the present invention the wavefront being analyzed includes at least one one-dimensional component, the transform applied to the wavefront being analyzed is a one-dimensional Fourier transform, performed in a dimension perpendicular to a direction of propagation of the wavefront being analyzed, thereby to obtain at least one one-dimensional component of the transformed wavefront in the dimension perpendicular to the direction of propagation. The plurality of differently phase changed transformed wavefronts are obtained by applying the plurality of different phase changes to each of the at least one one-dimensional component, thereby obtaining at least one one-dimensional component of the plurality of phase changed transformed wavefronts and the plurality of intensity maps are employed to obtain an output indicating amplitude and phase of the at least one one-dimensional component of the wavefront being analyzed.

[0081] Preferably, the plurality of different phase changes is applied to each of the one-dimensional component by providing a relative movement between the wavefront being analyzed and an element. Preferably, the element generates spatially varying, time-constant phase changes, the relative movement being in an additional dimension which is perpendicular both to the direction of propagation and to the dimension perpendicular to the direction of propagation.

[0082] Further in accordance with a preferred embodiment of the present invention the wavefront being analyzed includes a plurality of different wavelength components, the plurality of different phase changes are applied to the plurality of different wavelength components of each of the plurality of one-dimensional components of the wavefront being analyzed and the step of obtaining a plurality of intensity maps includes dividing the plurality of one-dimensional components of the plurality of phase changed transformed wavefronts into separate wavelength components.

[0083] Still further in accordance with a preferred embodiment of the present invention dividing the plurality of one-dimensional components of the plurality of phase changed transformed wavefronts into separate wavelength components is achieved by passing the plurality of phase changed transformed wavefronts through a dispersion element.

[0084] Further in accordance with a preferred embodiment of the present invention the transform applied to the wavefront being analyzed includes an additional Fourier transform to minimize cross-talk between different one-dimensional components of the wavefront being analyzed.

[0085] There is provided in accordance with another preferred embodiment of the present invention a method of surface mapping. The method includes obtaining a surface mapping wavefront having an amplitude and a phase, by reflecting radiation from a surface and analyzing the surface mapping wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to the surface mapping wavefront, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the surface mapping wavefront.

[0086] There is further provided in accordance with a preferred embodiment of the present invention an apparatus for surface mapping. The apparatus includes a wavefront obtainer operating to obtain a surface mapping wavefront having an amplitude and a phase, by reflecting radiation from a surface, a wavefront analyzer, analyzing the surface mapping wavefront and including a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to the surface mapping wavefront, an intensity map generator operating to provide a plurality of intensity maps of the plurality of phase changed transformed wavefronts and an intensity map utilizer, the plurality of intensity maps provide an output indicating the amplitude and phase of the surface mapping wavefront.

[0087] Further in accordance with a preferred embodiment of the present invention the radiation reflected from the surface has a narrow band about a given wavelength, causing the phase of the surface mapping wavefront to be proportional to geometrical variations in the surface, the proportion being an inverse linear function of the wavelength.

[0088] Still further in accordance with a preferred embodiment of the present invention the radiation reflected from the surface has at least two narrow bands, each centered about a different wavelength, providing at least two wavelength components in the surface mapping wavefront and at least two indications of the phase of the surface mapping wavefront, thereby enabling an enhanced mapping of the surface to be obtained by avoiding an ambiguity in the mapping which exceeds the larger of the different wavelengths about which the two narrow bands are centered.

[0089] Additionally in accordance with a preferred embodiment of the present invention the step of obtaining a plurality of differently phase changed transformed wavefronts includes applying a transform to the surface mapping wavefront, thereby to obtain a transformed wavefront and applying a plurality of different phase changes, including spatial phase changes, to the transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.

[0090] Further in accordance with a preferred embodiment of the present invention the transform applied to the surface mapping wavefront is a Fourier transform, the plurality of different phase changes includes at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of the transformed wavefront, the plurality of intensity maps includes at least three intensity maps. Preferably, the step of employing the plurality of intensity maps to obtain an output indicates the amplitude and phase of the surface mapping wavefront and includes expressing the surface mapping wavefront as a first complex function which has an amplitude and phase identical to the amplitude and phase of the surface mapping wavefront, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change. Preferably, expressing each of the plurality of intensity maps as a third function of: the amplitude of the surface mapping wavefront, the absolute value of the second complex function, a difference between the phase of the surface mapping wavefront and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes, corresponding to one of the at least three intensity maps, solving the third function to obtain the amplitude of the surface mapping wavefront, the absolute value of the second complex function and the difference between the phase of the surface mapping wavefront and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the surface mapping wavefront by adding the phase of the second complex function to the difference between the phase of the surface mapping wavefront and phase of the second complex function.

[0091] Preferably, the surface mapping wavefront includes a plurality of different wavelength components. The plurality of differently phase changed transformed wavefronts are preferably obtained by: transforming the surface mapping wavefront thereby obtaining a transformed wavefront including a plurality of different wavelength components and applying a phase change to the plurality of different wavelength components of the transformed wavefront by passing the transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.

[0092] There is also provided in accordance with yet another preferred embodiment of the present invention a method of inspecting an object. The method includes obtaining an object inspection wavefront which has an amplitude and a phase, by transmitting radiation through the object and analyzing the object inspection wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to the object inspection wavefront, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the object inspection wavefront.

[0093] There is further provided in accordance with a preferred embodiment of the present invention an apparatus for inspecting an object. The apparatus includes a wavefront obtainer operating to obtain an object inspection wavefront which has an amplitude and a phase, by transmitting radiation through the object, a wavefront analyzer, analyzing the object inspection wavefront and including a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to the object inspection wavefront, an intensity map generator operating to provide a plurality of intensity maps of the plurality of phase changed transformed wavefronts and an intensity map utilizer, employing the plurality of intensity maps to provide an output indicating the amplitude and phase of the object inspection wavefront.

[0094] Preferably, when the object is substantially uniform in material and other optical properties, the phase of the object inspection wavefront is proportional to the object thickness.

[0095] Additionally, when the object is substantially uniform in thickness, the phase of the object inspection wavefront is proportional to optical properties of the object.

[0096] Further in accordance with a preferred embodiment of the present invention the radiation has at least two narrow bands, each centered about a different wavelength, providing at least two wavelength components in the object inspection wavefront and at least two indications of the phase of the object inspection wavefront, thereby enabling an enhanced mapping of thickness of the object to be inspected by avoiding an ambiguity in the mapping which exceeds the larger of the different wavelengths about which the two narrow bands are centered.

[0097] Still further in accordance with a preferred embodiment of the present invention the method of obtaining a plurality of differently phase changed transformed wavefronts includes applying a transform to the object inspection wavefront, thereby obtaining a transformed wavefront and applying a plurality of different phase changes, including spatial phase changes, to the transformed wavefront, thereby obtaining a plurality of differently phase changed transformed wavefronts.

[0098] Further in accordance with a preferred embodiment of the present invention the transform applied to the object inspection wavefront is a Fourier transform, the plurality of different phase changes includes at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of the transformed wavefront. Preferably, the plurality of intensity maps includes at least three intensity maps and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the object inspection wavefront and includes: expressing the object inspection wavefront as a first complex function which has an amplitude and phase identical to the amplitude and phase of the object inspection wavefront, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change, expressing each of the plurality of intensity maps as a third function of: the amplitude of the object inspection wavefront, the absolute value of the second complex function, a difference between the phase of the object inspection wavefront and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes, corresponding to one of the at least three intensity maps, solving the third function to obtain the amplitude of the object inspection wavefront, the absolute value of the second complex function and the difference between the phase of the object inspection wavefront and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the object inspection wavefront by adding the phase of the second complex function to the difference between the phase of the object inspection wavefront and phase of the second complex function.

[0099] Still further in accordance with a preferred embodiment of the present invention the object inspection wavefront includes a plurality of different wavelength components. The plurality of differently phase changed transformed wavefronts are preferably obtained by: transforming the object inspection wavefront thereby obtaining a transformed wavefront including a plurality of different wavelength components and applying a phase change to the plurality of different wavelength components of the transformed wavefront by reflecting the transformed wavefront from a spatially varying surface.

[0100] There is also provided in accordance with yet another preferred embodiment of the present invention a method of spectral analysis. The method includes obtaining a spectral analysis wavefront having an amplitude and a phase, by causing radiation to impinge on an object, analyzing the spectral analysis wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to the spectral analysis wavefront which has an amplitude and a phase, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the spectral analysis wavefront and employing the output indicating the amplitude and phase to obtain an output indicating spectral content of the radiation.

[0101] There is provided in accordance with a further preferred embodiment of the present invention an apparatus for spectral analysis. The apparatus includes a wavefront obtainer operating to obtain a spectral analysis wavefront having an amplitude and a phase, by causing radiation to impinge on an object, a wavefront analyzer, analyzing the spectral analysis wavefront, including a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to the spectral analysis wavefront which has an amplitude and a phase, an intensity map generator operating to provide a plurality of intensity maps of the plurality of phase changed transformed wavefronts, an intensity map utilizer, employing the plurality of intensity maps to provide an output indicating the amplitude and phase of the spectral analysis wavefront and a phase and amplitude utilizer, employing the output indicating the amplitude and phase to obtain an output indicating spectral content of the radiation.

[0102] Further in accordance with a preferred embodiment of the present invention and wherein obtaining the spectral analysis wavefront is effected by reflecting the radiation from the object.

[0103] Still further in accordance with a preferred embodiment of the present invention and wherein obtaining the spectral analysis wavefront is effected by transmitting the radiation through the object.

[0104] Additionally in accordance with a preferred embodiment of the present invention the radiation is substantially of a single wavelength, the phase of the spectral analysis wavefront is inversely proportional to the single wavelength, and is related to at least one of a surface characteristic and thickness of the impinged object.

[0105] Still further in accordance with a preferred embodiment of the present invention the step of employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the spectral analysis wavefront includes: expressing the plurality of intensity maps as at least one mathematical function of phase and amplitude of the spectral analysis wavefront and of the plurality of different phase changes, wherein at least the phase is unknown and a function generating the plurality of phase changed transformed wavefronts is known and employing the at least one mathematical function to obtain an output indicating at least the phase.

[0106] Additionally in accordance with a preferred embodiment of the present invention the step of obtaining a plurality of differently phase changed transformed wavefronts includes applying a transform to the spectral analysis wavefront, thereby obtaining a transformed wavefront and applying a plurality of different phase changes, including spatial phase changes, to the transformed wavefront, thereby obtaining a plurality of differently phase changed transformed wavefronts.

[0107] Further in accordance with a preferred embodiment of the present invention the transform applied to the spectral analysis wavefront is a Fourier transform, the plurality of different phase changes includes at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of the transformed wavefront. Preferably, the plurality of intensity maps includes at least three intensity maps and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the spectral analysis wavefront includes: expressing the spectral analysis wavefront as a first complex function which has an amplitude and phase identical to the amplitude and phase of the spectral analysis wavefront, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change, expressing each of the plurality of intensity maps as a third function of: the amplitude of the spectral analysis wavefront, the absolute value of the second complex function, a difference between the phase of the spectral analysis wavefront and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes, corresponding to one of the at least three intensity maps, solving the third function to obtain the amplitude of the spectral analysis wavefront, the absolute value of the second complex function and the difference between the phase of the spectral analysis wavefront and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the spectral analysis wavefront by adding the phase of the second complex function to the difference between the phase of the spectral analysis wavefront and phase of the second complex function.

[0108] Further in accordance with a preferred embodiment of the present invention the spectral analysis wavefront includes a plurality of different wavelength components and the plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to the plurality of different wavelength components of the spectral analysis wavefront.

[0109] There is further provided in accordance with a preferred embodiment of the present invention a method of phase change analysis. The method includes obtaining a phase change analysis wavefront which has an amplitude and a phase, applying a transform to the phase change analysis wavefront thereby to obtain a transformed wavefront, applying a plurality of different phase changes to the transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indication of differences between the plurality of different phase changes applied to the transformed phase change analysis wavefront.

[0110] There is also provided in accordance with yet another preferred embodiment of the present invention an apparatus for phase change analysis. The apparatus includes a wavefront obtainer, operating to obtain a phase change analysis wavefront which has an amplitude and a phase, a transform applier, applying a transform to the phase change analysis wavefront thereby to obtain a transformed wavefront, a phase change applier, applying at least one phase change to the transformed wavefront, thereby to obtain at least one phase changed transformed wavefront, an intensity map generator operating to provide at least one intensity map of the phase changed transformed wavefront and an intensity map utilizer, employing the plurality of intensity maps to provide an output indication of differences between the plurality of different phase changes applied to the transformed phase change analysis wavefront.

[0111] Typically, when lateral shifts appear in the plurality of different phase changes, corresponding changes appear in the plurality of intensity maps and the step of employing the plurality of intensity maps results in obtaining an indication of the lateral shifts.

[0112] Still further in accordance with a preferred embodiment of the present invention the step of employing the plurality of intensity maps to obtain an output indication of differences between the plurality of different phase changes applied to the transformed phase change analysis wavefront includes: expressing the plurality of intensity maps as at least one mathematical function of phase and amplitude of the phase change analysis wavefront and of the plurality of different phase changes, where at least the phase and amplitude are known and the plurality of different phase changes are unknown and employing the mathematical function to obtain an output indicating the differences between the plurality of different phase changes.

[0113] There is further provided in accordance with yet a further preferred embodiment of the present invention a method of phase change analysis. The method includes obtaining a phase change analysis wavefront which has an amplitude and a phase, applying a transform to the phase change analysis wavefront thereby to obtain a transformed wavefront, applying at least one phase change to the transformed wavefront, thereby to obtain at least one phase changed transformed wavefront, obtaining at least one intensity map of the at least one phase changed transformed wavefront and employing the intensity map to obtain an output indication of the at least one phase change applied to the transformed phase change analysis wavefront.

[0114] There is also provided in accordance with yet another preferred embodiment of the present invention an apparatus for phase change analysis. The apparatus includes a wavefront obtainer, operating to obtain a phase change analysis wavefront which has an amplitude and a phase, a transform applier, applying a transform to the phase change analysis wavefront thereby to obtain a transformed wavefront, a phase change applier, applying at least one phase change to the transformed wavefront, thereby to obtain at least one phase changed transformed wavefront, an intensity map generator operating to provide at least one intensity map of the phase changed transformed wavefront and an intensity map utilizer, employing the intensity map to provide an output indication of the phase change applied to the transformed phase change analysis wavefront.

[0115] Preferably, the phase change is a phase delay, having a value selected from a plurality of pre-determined values, and the output indication of the phase change includes the value of the phase delay.

[0116] There is also provided in accordance with a preferred embodiment of the present invention a method of stored data retrieval. The method includes obtaining a stored data retrieval wavefront which has an amplitude and a phase, by reflecting radiation from the media in which information is encoded, by selecting the height of the media at each of a multiplicity of different locations on the media. Preferably, analyzing the stored data retrieval wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to the stored data retrieval wavefront, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an indication of the amplitude and phase of the stored data retrieval wavefront and employing the indication of the amplitude and phase to obtain the information.

[0117] There is further provided in accordance with yet another preferred embodiment of the present invention an apparatus for stored data retrieval. The apparatus includes a wavefront obtainer operating to obtain a stored data retrieval wavefront which has an amplitude and a phase, by reflecting radiation from the media in which information is encoded by selecting the height of the media at each of a multiplicity of different locations on the media, a wavefront analyzer, analyzing the stored data retrieval wavefront and including a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to the stored data retrieval wavefront, an intensity map generator operating to obtain a plurality of intensity maps of the plurality of phase changed transformed wavefronts and an intensity map utilizer, employing the plurality of intensity maps to provide an indication of the amplitude and phase of the stored data retrieval wavefront and a phase and amplitude utilizer employing the indication of the amplitude and phase to provide the information.

[0118] Preferably, the step of obtaining a plurality of differently phase changed transformed wavefronts includes: applying a transform to the stored data retrieval wavefront thereby to obtain a transformed wavefront and applying a plurality of different phase changes to the transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.

[0119] Further in accordance with a preferred embodiment of the present invention the transform applied to the stored data retrieval wavefront is a Fourier transform, the plurality of different phase changes includes at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of the transformed wavefront, the plurality of intensity maps includes at least three intensity maps and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the stored data retrieval wavefront includes: expressing the stored data retrieval wavefront as a first complex function which has an amplitude and phase identical to the amplitude and phase of the stored data retrieval wavefront, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change, expressing each of the plurality of intensity maps as a third function of: the amplitude of the stored data retrieval wavefront, the absolute value of the second complex function, a difference between the phase of the stored data retrieval wavefront and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes, corresponding to one of the at least three intensity maps, solving the third function to obtain the amplitude of the stored data retrieval wavefront, the absolute value of the second complex function and the difference between the phase of the stored data retrieval wavefront and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the stored data retrieval wavefront by adding the phase of the second complex function to the difference between the phase of the stored data retrieval wavefront and phase of the second complex function.

[0120] Still further in accordance with a preferred embodiment of the present invention the stored data retrieval wavefront includes at least one one-dimensional component, the transform applied to the data retrieval wavefront is a one-dimensional Fourier transform, performed in a dimension perpendicular to a direction of propagation of the data retrieval wavefront, thereby to obtain at least one one-dimensional component of the transformed wavefront in the dimension perpendicular to the direction of propagation, the plurality of differently phase changed transformed wavefronts are obtained by applying the plurality of different phase changes to each of the one-dimensional component, thereby obtaining at least one one-dimensional component of the plurality of phase changed transformed wavefronts and the plurality of intensity maps are employed to obtain an output indicating amplitude and phase of the one-dimensional component of the data retrieval wavefront.

[0121] Preferably, the plurality of different phase changes is applied to each of the at least one one-dimensional component by providing a relative movement between the media and a component generating spatially varying, time-constant phase changes, the relative movement being in a dimension perpendicular to the direction of propagation and to the dimension perpendicular to the direction of propagation.

[0122] Additionally in accordance with a preferred embodiment of the present invention the information is encoded on the media whereby: an intensity value is realized by reflection of light from each location on the media to lie within a predetermined range of values, the range corresponding an element of the information stored at the location and by employing the plurality of intensity maps, multiple intensity values are realized for each location, providing multiple elements of information for each location on the media.

[0123] Preferably, the plurality of differently phase changed transformed wavefronts include a plurality of wavefronts whose phase has been changed by applying an at least time varying phase change function to the stored data retrieval wavefront.

[0124] Further in accordance with a preferred embodiment of the present invention the stored data retrieval wavefront includes a plurality of different wavelength components and the plurality of differently phase changed transformed wavefronts are obtained by applying at least one phase change to the plurality of different wavelength components of the stored data retrieval wavefront.

[0125] Further in accordance with a preferred embodiment of the present invention the radiation which is reflected from the media includes a plurality of different wavelength components, resulting in the stored data retrieval wavefront including a plurality of different wavelength components and the plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to the plurality of different wavelength components of the stored data retrieval wavefront.

[0126] Still further in accordance with a preferred embodiment of the present invention the information encoded by selecting the height of the media at each of a multiplicity of different locations on the media is also encoded by selecting the reflectivity of the media at each of a plurality of different locations on the media and employing the indication of the amplitude and phase to obtain the information includes employing the indication of the phase to obtain the information encoded by selecting the height of the media and employing the indication of the amplitude to obtain the information encoded by selecting the reflectivity of the media.

[0127] There is provided in accordance with another preferred embodiment of the present invention a method of 3-dimensional imaging. The method includes obtaining a 3-dimensional imaging wavefront, which has an amplitude and a phase, by reflecting radiation from an object to be viewed and analyzing the 3-dimensional imaging wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to the 3-dimensional imaging wavefront, obtaining a plurality of intensity maps of the plurality of differently phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the 3-dimensional imaging wavefront.

[0128] There is further provided in accordance with a preferred embodiment of the present invention an apparatus for 3-dimensional imaging. The apparatus includes a wavefront obtainer operating to obtain a 3-dimensional imaging wavefront, which has an amplitude and a phase, by reflecting radiation from an object to be viewed, a wavefront analyzer, analyzing the 3-dimensional imaging wavefront including a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to the 3-dimensional imaging wavefront, an intensity map generator operative to provide a plurality of intensity maps of the plurality of differently phase changed transformed wavefronts and an intensity map utilizer, employing the plurality of intensity maps to provide an output indicating the amplitude and phase of the 3-dimensional imaging wavefront.

[0129] Further in accordance with a preferred embodiment of the present invention the radiation reflected from the object has a narrow band about a given wavelength, causing the phase of the 3-dimensional imaging wavefront to be proportional to geometrical variations in the object, the proportion being an inverse linear function of the wavelength.

[0130] Additionally in accordance with a preferred embodiment of the present invention the step of obtaining a plurality of differently phase changed transformed wavefronts includes applying a transform to the 3-dimensional imaging wavefront, thereby to obtain a transformed wavefront and applying a plurality of different phase changes, including spatial phase changes, to the transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.

[0131] Still further in accordance with a preferred embodiment of the present invention the 3-dimensional imaging wavefront includes a plurality of different wavelength components and the plurality of differently phase changed transformed wavefronts are obtained by: transforming the 3-dimensional imaging wavefront, thereby obtaining a transformed wavefront including a plurality of different wavelength components and applying phase changes to the plurality of different wavelength components of the transformed wavefront by passing the transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.

[0132] There is also provided in accordance with yet another preferred embodiment of the present invention a method of wavefront analysis. The method includes obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output indicating at least the phase of the wavefront being analyzed by combining the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, obtaining at least an output indicative of the phase of the wavefront being analyzed from each of the second plurality of combined intensity maps and combining the outputs to provide at least an enhanced indication of phase of the wavefront being analyzed.

[0133] There is also provided in accordance with yet another preferred embodiment of the present invention an apparatus wavefront analysis. The apparatus includes a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, an intensity map generator operating to obtain a plurality of intensity maps of the plurality of phase changed transformed wavefronts and an intensity map utilizer, employing the plurality of intensity maps to obtain an output indicating at least amplitude of the wavefront being analyzed and including an intensity combiner operating to combine the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, an indication provider operating to provide at least an output indicative of the amplitude of the wavefront being analyzed from each of the second plurality of combined intensity maps and an enhanced indication provider, combining the outputs to provide at least an enhanced indication of amplitude of the wavefront being analyzed.

[0134] There is provided in accordance with a further preferred embodiment of the present invention a method of wavefront analysis. The method includes obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefront and employing the plurality of intensity maps to obtain an output indicating at least amplitude of the wavefront being analyzed by combining the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, obtaining at least an output indicative of the amplitude of the wavefront being analyzed from each of the second plurality of combined intensity maps and combining the outputs to provide at least an enhanced indication of amplitude of the wavefront being analyzed.

[0135] There is provided in accordance with a preferred embodiment of the present invention an apparatus for wavefront analysis. The apparatus includes a wavefront transformer operating to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, an intensity map generator operating to provide a plurality of intensity maps of the plurality of phase changed transformed wavefronts and an intensity map utilizer, employing the plurality of intensity maps to provide an output indicating at least the phase of the wavefront being analyzed. Preferably, the apparatus also includes an intensity map expresser, expressing the plurality of intensity maps as a function of: amplitude of the wavefront being analyzed, phase of the wavefront being analyzed and a phase change function characterizing the plurality of differently phase changed transformed wavefronts, a complex function definer, defining a complex function of: the amplitude of the wavefront being analyzed, the phase of the wavefront being analyzed and the phase change function characterizing the plurality of differently phase changed transformed wavefronts, the complex function being characterized in that the intensity at each location in the plurality of intensity maps is a function predominantly of a value of the complex function at the location and of the amplitude and the phase of the wavefront being analyzed at the location. The apparatus also typically, includes complex function expresser, expressing the complex function as a function of the plurality of intensity maps and a phase obtainer, obtaining values for the phase by employing the complex function expressed as a function of the plurality of intensity maps.

[0136] There is also provided in accordance with another preferred embodiment of the present invention a method of wavefront analysis. The method includes obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to provide an output indicating at least the phase of the wavefront being analyzed by: expressing the plurality of intensity maps as a function of: amplitude of the wavefront being analyzed, phase of the wavefront being analyzed and a phase change function characterizing the plurality of differently phase changed transformed wavefronts. Additionally, defining a complex function of: the amplitude of the wavefront being analyzed, the phase of the wavefront being analyzed and the phase change function characterizing the plurality of differently phase changed transformed wavefronts, the complex function being characterized in that the intensity at each location in the plurality of intensity maps is a function predominantly of a value of the complex function at the location and of the amplitude and the phase of the wavefront being analyzed at the location, expressing the complex function as a function of the plurality of intensity maps and obtaining values for the phase by employing the complex function expressed as a function of the plurality of intensity maps.

[0137] There is further provided in accordance with yet a further preferred embodiment of the present invention a method of wavefront analysis. The method includes applying a Fourier transform to a wavefront being analyzed which has an amplitude and a phase, thereby obtaining a transformed wavefront, applying a spatially uniform time-varying spatial phase change to part of the transformed wavefront, thereby to obtain at least three differently phase changed transformed wavefronts, applying a second Fourier transform to obtain at least three intensity maps of the at least three phase changed transformed wavefronts and employing the at least three intensity maps to obtain an output indicating at least one of the phase and the amplitude of the wavefront being analyzed by: expressing the wavefront being analyzed as a first complex function which has an amplitude and phase identical to the amplitude and phase of the wavefront being analyzed, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change, expressing each of the plurality of intensity maps as a third function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function, a difference between the phase of the wavefront being analyzed and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes, which each correspond to one of the at least three intensity maps, solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function, solving the second complex function to obtain the phase of the second complex function and obtaining the phase of the wavefront being analyzed by adding the phase of the second complex function to the difference between the phase of the wavefront being analyzed and phase of the second complex function.

[0138] There is further provided in accordance with yet a further preferred embodiment of the present invention an apparatus for wavefront analysis. The apparatus includes a first transform applier, applying a Fourier transform to a wavefront being analyzed which has an amplitude and a phase thereby to obtain a transformed wavefront, a phase change applier, applying a spatially uniform time-varying spatial phase change to part of the transformed wavefront, thereby obtaining at least three differently phase changed transformed wavefronts, a second transform applier, applying a second Fourier transform to the at least three differently phase changed transformed wavefronts, thereby obtaining at least three intensity maps. The apparatus also typically includes an intensity map utilizer, employing the at least three intensity maps to provide an output indicating the phase and the amplitude of the wavefront being analyzed and a wavefront expresser, expressing the wavefront being analyzed as a first complex function which has an amplitude and phase identical to the amplitude and phase of the wavefront being analyzed, a first intensity map expresser, expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial phase change. Preferably, the apparatus also includes a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial phase change, a second intensity map expresser, expressing each of the plurality of intensity maps as a third function of: the amplitude of the wavefront being analyzed, the absolute value of the second complex function, a difference between the phase of the wavefront being analyzed and the phase of the second complex function and a known phase delay produced by one of the at least three different phase changes, which each correspond to one of the at least three intensity maps. The apparatus further typically includes a first function solver, solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function, a second function solver, solving the second complex function to obtain the phase of the second complex function and a phase obtainer, obtaining the phase of the wavefront being analyzed by adding the phase of the second complex function to the difference between the phase of the wavefront being analyzed and the phase of the second complex function.

BRIEF DESCRIPTION OF THE DRAWINGS

[0139] The present invention will be understood and appreciated more fully from the following detailed description, taken in conjunction with the drawings in which:

[0140]FIG. 1A is a simplified partially schematic, partially pictorial illustration of wavefront analysis functionality operative in accordance with a preferred embodiment of the present invention;

[0141]FIG. 1B is a simplified partially schematic, partially block diagram illustration of a wavefront analysis system suitable for carrying out the functionality of FIG. 1A in accordance with a preferred embodiment of the present invention;

[0142]FIG. 2 is a simplified functional block diagram illustration of the functionality of FIG. 1A where time-varying phase changes are applied to a transformed wavefront;

[0143]FIG. 3 is a simplified functional block diagram illustration of the functionality of FIG. 1A where time-varying phase changes are applied to a wavefront prior to transforming thereof;

[0144]FIG. 4 is a simplified functional block diagram illustration of the functionality of FIG. 2 where time-varying, non-spatially varying spatial phase changes are applied to a transformed wavefront;

[0145]FIG. 5 is a simplified functional block diagram illustration of the functionality of FIG. 3 where time-varying, non-spatially varying spatial phase changes are applied to a wavefront prior to transforming thereof;

[0146]FIG. 6 is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different wavelength components of a transformed wavefront;

[0147]FIG. 7 is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different wavelength components of a wavefront prior to transforming thereof;

[0148]FIG. 8 is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different polarization components of a transformed wavefront;

[0149]FIG. 9 is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different polarization components of a wavefront prior to transforming thereof;

[0150]FIG. 10A is a simplified functional block diagram illustration of the functionality of FIG. 1A where a wavefront being analyzed comprises at least one one-dimensional component;

[0151]FIG. 10B is a simplified partially schematic, partially pictorial illustration of a wavefront analysis system suitable for carrying out the functionality of FIG. 10A in accordance with a preferred embodiment of the present invention;

[0152]FIG. 11 is a simplified functional block diagram illustration of the functionality of FIG. 1A where an additional transform is applied following the application of spatial phase changes;

[0153]FIG. 12 is a simplified functional block diagram illustration of the functionality of FIG. 1A, wherein intensity maps are employed to provide information about a wavefront being analyzed, such as indications of amplitude and phase of the wavefront;

[0154]FIG. 13 is a simplified functional block diagram illustration of part of the functionality of FIG. 1A, wherein the transform applied to the wavefront being analyzed is a Fourier transform, wherein at least three different spatial phase changes are applied to a transformed wavefront, and wherein at least three intensity maps are employed to obtain indications of at least the phase of a wavefront;

[0155]FIG. 14 is a simplified partially schematic, partially pictorial illustration of part of one preferred embodiment of a wavefront analysis system of the type shown in FIG. 1B;

[0156]FIG. 15 is a simplified partially schematic, partially pictorial illustration of a system for surface mapping employing the functionality and structure of FIGS. 1A and 1B;

[0157]FIG. 16 is a simplified partially schematic, partially pictorial illustration of a system for object inspection employing the functionality and structure of FIGS. 1A and 1B;

[0158]FIG. 17 is a simplified partially schematic, partially pictorial illustration of a system for spectral analysis employing the functionality and structure of FIGS. 1A and 1B;

[0159]FIG. 18 is a simplified partially schematic, partially pictorial illustration of a system for phase-change analysis employing the functionality and structure of FIGS. 1A and 1B;

[0160]FIG. 19 is a simplified partially schematic, partially pictorial illustration of a system for stored data retrieval employing the functionality and structure of FIGS. 1A and 1B;

[0161]FIG. 20 is a simplified partially schematic, partially pictorial illustration of a system for 3-dimensional imaging employing the functionality and structure of FIGS. 1A and 1B;

[0162]FIG. 21A is a simplified partially schematic, partially pictorial illustration of wavefront analysis functionality operative in accordance with another preferred embodiment of the present invention;

[0163]FIG. 21B is a simplified partially schematic, partially block diagram illustration of a wavefront analysis system suitable for carrying out the functionality of FIG. 21A in accordance with another preferred embodiment of the present invention; and

[0164]FIG. 22 is a simplified partially schematic, partially pictorial illustration of a system for surface mapping employing the functionality and structure of FIGS. 21A and 21B.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

[0165] Reference is now made to FIG. 1A, which is a simplified partially schematic, partially pictorial illustration of wavefront analysis functionality operative in accordance with a preferred embodiment of the present invention. The functionality of FIG. 1A can be summarized as including the following sub-functionalities:

[0166] A. obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, which has an amplitude and a phase;

[0167] B. obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts; and

[0168] C. employing the plurality of intensity maps to obtain an output indicating at least one and possibly both of the phase and the amplitude of the wavefront being analyzed.

[0169] As seen in FIG. 1A, the first sub-functionality, designated “A” may be realized by the following functionalities:

[0170] A wavefront, which may be represented by a plurality of point sources of light, is generally designated by reference numeral 100. Wavefront 100 has a phase characteristic which is typically spatially non-uniform, shown as a solid line and indicated generally by reference numeral 102. Wavefront 100 also has an amplitude characteristic which is also typically spatially non-uniform, shown as a dashed line and indicated generally by reference numeral 103. Such a wavefront may be obtained in a conventional manner by receiving light from any object, such as by reading an optical disk, for example a DVD or compact disk 104.

[0171] A principal purpose of the present invention is to measure the phase characteristic, such as that indicated by reference numeral 102, which is not readily measured. Another purpose of the present invention is to measure the amplitude characteristic, such as that indicated by reference numeral 103 in an enhanced manner. A further purpose of the present invention is to measure both the phase characteristic 102 and the amplitude characteristic 103. While there exist various techniques for carrying out such measurements, the present invention provides a methodology which is believed to be superior to those presently known, inter alia due to its relative insensitivity to noise.

[0172] A transform, indicated here symbolically by reference numeral 106, is applied to the wavefront being analyzed 100, thereby to obtain a transformed wavefront. A preferred transform is a Fourier transform. The resulting transformed wavefront is symbolically indicated by reference numeral 108.

[0173] A plurality of different phase changes, preferably spatial phase changes, represented by optical path delays 110, 112 and 114 are applied to the transformed wavefront 108, thereby to obtain a plurality of differently phase changed transformed wavefronts, represented by reference numerals 120, 122 and 124 respectively. It is appreciated that the illustrated difference between the individual ones of the plurality of differently phase changed transformed wavefronts is that portions of the transformed wavefront are delayed differently relative to the remainder thereof. The difference in the phase changes, which are applied to the transformed wavefront 108, is represented in FIG. 1A by the change in thickness of the optical path delays 110, 112 and 114.

[0174] As seen in FIG. 1A, the second sub-functionality, designated “B”, may be realized by applying a transform, preferably a Fourier transform, to the plurality of differently phase changed transformed wavefronts. Alternatively, the sub-functionality B may be realized without the use of a Fourier transform, such as by propagation of the differently phase changed transformed wavefronts over an extended space. Finally, functionality B requires detection of the intensity characteristics of plurality of differently phase changed transformed wavefronts. The outputs of such detection are the intensity maps, examples of which are designated by reference numerals 130, 132 and 134.

[0175] As seen in FIG. 1A, the third sub-functionality. designated “C” may be realized by the following functionalities:

[0176] expressing, such as by employing a computer 136, the plurality of intensity maps, such as maps 130, 132 and 134, as at least one mathematical function of phase and amplitude of the wavefront being analyzed and of the plurality of different phase changes, wherein at least one and possibly both of the phase and the amplitude are unknown and the plurality of different phase changes, typically represented by optical path delays 110, 112 and 114 to the transformed wavefront 108, are known; and

[0177] employing, such as by means of the computer 136, the at least one mathematical function to obtain an indication of at least one and possibly both of the phase and the amplitude of the wavefront being analyzed, here represented by the phase function designated by reference numeral 138 and the amplitude function designated by reference numeral 139, which, as can be seen, respectively represent the phase characteristics 102 and the amplitude characteristics 103 of the wavefront 100. In this example, wavefront 100 may represent the information contained in the compact disk or DVD 104.

[0178] In accordance with an embodiment of the present invention, the plurality of intensity maps comprises at least four intensity maps. In such a case, employing the plurality of intensity maps to obtain an output indicating at least the phase of the wavefront being analyzed includes employing a plurality of combinations, each of at least three of the plurality of intensity maps, to provide a plurality of indications at least of the phase of the wavefront being analyzed.

[0179] Preferably, the methodology also includes employing the plurality of indications of at least the phase of the wavefront being analyzed to provide an enhanced indication at least of the phase of the wavefront being analyzed.

[0180] Also in accordance with an embodiment of the present invention, the plurality of intensity maps comprises at least four intensity maps. In such a case, employing the plurality of intensity maps to obtain an output indicating at least the amplitude of the wavefront being analyzed includes employing a plurality of combinations, each of at least three of the plurality of intensity maps, to provide a plurality of indications at least of the amplitude of the wavefront being analyzed.

[0181] Preferably, the methodology also includes employing the plurality of indications of at least the amplitude of the wavefront being analyzed to provide an enhanced indication at least of the amplitude of the wavefront being analyzed.

[0182] It is appreciated that in this manner, enhanced indications of both phase and amplitude of the wavefront may be obtained.

[0183] In accordance with a preferred embodiment of the present invention, at least some of the plurality of indications of the amplitude and phase are at least second order indications of the amplitude and phase of the wavefront being analyzed.

[0184] In accordance with one preferred embodiment of the present invention, the plurality of intensity maps are employed to provide an analytical output indicating the amplitude and phase.

[0185] Preferably, the phase changed transformed wavefronts are obtained by interference of the wavefront being analyzed along a common optical path.

[0186] In accordance with one preferred embodiment of the present invention, the plurality of differently phase changed transformed wavefronts are realized in a manner substantially different from performing a delta-function phase change to the transformed wavefront, whereby a delta-function phase change is applying a uniform phase delay to a small spatial region, having the characteristics of a delta-function, of the transformed wavefront.

[0187] In accordance with another preferred embodiment of the present invention, the plurality of intensity maps are employed to obtain an output indicating the phase of the wavefront being analyzed, which is substantially free from halo and shading off distortions, which are characteristic of many of the existing ‘phase-contrast’ methods.

[0188] In accordance with another embodiment of the present invention the output indicating the phase of the wavefront being analyzed may be processed to obtain the polarization mode of the wavefront being analyzed.

[0189] In accordance with still another embodiment of the present invention, the plurality of intensity maps may be employed to obtain an output indicating the phase of the wavefront being analyzed by combining the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, obtaining at least an output indicative of the phase of the wavefront being analyzed from each of the second plurality of combined intensity maps and combining the outputs to provide an enhanced indication of the phase of the wavefront being analyzed.

[0190] In accordance with yet another embodiment of the present invention, the plurality of intensity maps may be employed to obtain an output indicating amplitude of the wavefront being analyzed by combining the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, obtaining at least an output indicative of the amplitude of the wavefront being analyzed from each of the second plurality of combined intensity maps and combining the outputs to provide an enhanced indication of the amplitude of the wavefront being analyzed.

[0191] Additionally in accordance with a preferred embodiment of the present invention, the foregoing methodology may be employed for obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output of an at least second order indication of phase of the wavefront being analyzed.

[0192] Additionally or alternatively in accordance with a preferred embodiment of the present invention, the foregoing methodology may be employed for obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts and employing the plurality of intensity maps to obtain an output of an at least second order indication of amplitude of the wavefront being analyzed.

[0193] In accordance with yet another embodiment of the present invention, the obtaining of the plurality of differently phase changed transformed wavefronts comprises applying a transform to the wavefront being analyzed, thereby to obtain a transformed wavefront, and then applying a plurality of different phase and amplitude changes to the transformed wavefront, where each of these changes can be a phase change, an amplitude change or a combined phase and amplitude change, thereby to obtain a plurality of differently phase and amplitude changed transformed wavefronts.

[0194] In accordance with yet another embodiment of the present invention, a wavefront being analyzed comprises at least two wavelength components. In such a case, obtaining a plurality of intensity maps also includes dividing the phase changed transformed wavefronts according to the at least two wavelength components in order to obtain at least two wavelength components of the phase changed transformed wavefronts and in order to obtain at least two sets of intensity maps, each set corresponding to a different one of the at least two wavelength components of the phase changed transformed wavefronts.

[0195] Subsequently, the plurality of intensity maps are employed to provide an output indicating the amplitude and phase of the wavefront being analyzed by obtaining an output indicative of the phase of the wavefront being analyzed from each of the at least two sets of intensity maps and combining the outputs to provide an enhanced indication of phase of the wavefront being analyzed. In the enhanced indication, there is no 2π ambiguity once the value of the phase exceeds 2π, which conventionally results when detecting a phase of a single wavelength wavefront.

[0196] It is appreciated that the wavefront being analyzed may be an acoustic radiation wavefront.

[0197] It is also appreciated that the wavefront being analyzed may be an electromagnetic radiation wavefront, of any suitable wavelength, such as visible light, infrared, ultra-violet and X-ray radiation.

[0198] It is further appreciated that wavefront 100 may be represented by a relatively small number of point sources and defined over a relatively small spatial region. In such a case, the detection of the intensity characteristics of the plurality of differently phase changed transformed wavefronts may be performed by a detector comprising only a single detection pixel or several detection pixels. Additionally, the output indicating at least one and possibly both of the phase and amplitude of the wavefront being analyzed, may be provided by computer 136 in a straight-forward manner.

[0199] Reference is now made to FIG. 1B, which is a simplified partially schematic, partially block diagram illustration of a wavefront analysis system suitable for carrying out the functionality of FIG. 1A in accordance with a preferred embodiment of the present invention. As seen in FIG. 1B, a wavefront, here designated by reference numeral 150 is focused, as by a lens 152, onto a phase manipulator 154, which is preferably located at the focal plane of lens 152. The phase manipulator 154 generates phase changes, and may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects.

[0200] A second lens 156 is arranged so as to image wavefront 150 onto a detector 158, such as a CCD detector. Preferably the second lens 156 is arranged such that the detector 158 lies in its focal plane. The output of detector 158 is preferably supplied to data storage and processing circuitry 160, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A.

[0201] Reference is now made to FIG. 2, which is a simplified functional block diagram illustration of the functionality of FIG. 1A where time-varying phase changes are applied to a transformed wavefront. As seen in FIG. 2, and as explained hereinabove with reference to FIG. 1A, a wavefront 200 is preferably transformed to provide a transformed wavefront 208.

[0202] A first phase change, preferably a spatial phase change, is applied to the transformed wavefront 208 at a first time T1, as indicated by reference numeral 210, thereby producing a phase changed transformed wavefront 212 at time T1. This phase changed transformed wavefront 212 is detected, as by detector 158 (FIG. 1B), producing an intensity map, an example of which is designated by reference numeral 214, which map is stored as by circuitry 160 (FIG. 1B).

[0203] Thereafter, a second phase change, preferably a spatial phase change, is applied to the transformed wavefront 208 at a second time T2, as indicated by reference numeral 220, thereby producing a phase changed transformed wavefront 222 at time T2. This phase changed transformed wavefront 222 is detected, as by detector 158 (FIG. 1B), producing an intensity map, an example of which is designated by reference numeral 224, which map is stored as by circuitry 160 (FIG. 1B).

[0204] Thereafter, a third phase change, preferably a spatial phase change, is applied to the transformed wavefront 208 at a third time T3, as indicated by reference numeral 230, thereby producing a phase changed transformed wavefront 232 at time T3. This phase changed transformed wavefront 232 is detected, as by detector 158 (FIG. 1B), producing an intensity map, an example of which is designated by reference numeral 234, which map is stored as by circuitry 160 (FIG. 1B).

[0205] It is appreciated that any suitable number of spatial phase changes may be made at successive times and stored for use in accordance with the present invention.

[0206] In accordance with a preferred embodiment of the present invention, at least some of the phase changes 210, 220 and 230, are spatial phase changes effected by applying a spatial phase change to part of the transformed wavefront 208.

[0207] In accordance with another preferred embodiment of the present invention, at least some of the phase changes 210, 220 and 230, are spatial phase changes, effected by applying a time-varying spatial phase change to part of the transformed wavefront 208.

[0208] In accordance with another preferred embodiment of the present invention, at least some of the phase changes 210, 220 and 230, are spatial phase changes, effected by applying a non time-varying spatial phase change to part of transformed wavefront 208, producing spatially phase changed transformed wavefronts 212, 222 and 232, which subsequently produce spatially varying intensity maps 214, 224 and 234 respectively.

[0209] Reference is now made to FIG. 3, which is a simplified functional block diagram illustration of the functionality of FIG. 1A where time-varying phase changes are applied to a wavefront prior to transforming thereof. As seen in FIG. 3, a first phase change, preferably a spatial phase change, is applied to a wavefront 300 at a first time T1, as indicated by reference numeral 310. Following application of the first phase change to wavefront 300, a transform, preferably a Fourier transform, is applied thereto, thereby producing a phase changed transformed wavefront 312 at time T1. This phase changed transformed wavefront 312 is detected, as by detector 158 (FIG. 1B), producing an intensity map, an example of which is designated by reference numeral 314, which map is stored as by circuitry 160 (FIG. 1B).

[0210] Thereafter, a second phase change, preferably a spatial phase change, is applied to wavefront 300 at a second time T2, as indicated by reference numeral 320. Following application of the second phase change to wavefront 300, a transform, preferably a Fourier transform, is applied thereto, thereby producing a phase changed transformed wavefront 322 at time T2. This phase changed transformed wavefront 322 is detected, as by detector 158 (FIG. 1B), producing an intensity map, an example of which is designated by reference numeral 324, which map is stored as by circuitry 160 (FIG. 1B).

[0211] Thereafter, a third phase change, preferably a spatial phase change, is applied to wavefront 300 at a third time T3, as indicated by reference numeral 330. Following application of the third phase change to wavefront 300, a transform, preferably a Fourier transform, is applied thereto, thereby producing a phase changed transformed wavefront 332 at time T3. This phase changed transformed wavefront 332 is detected, as by detector 158 (FIG. 1B), producing an intensity map, an example of which is designated by reference numeral 334, which map is stored as by circuitry 160 (FIG. 1B).

[0212] It is appreciated that any suitable number of spatial phase changes may be made at successive times and stored for use in accordance with the present invention.

[0213] In accordance with a preferred embodiment of the present invention, at least some of the phase changes 310, 320 and 330, are spatial phase changes effected by applying a spatial phase change to part of wavefront 300.

[0214] In accordance with another preferred embodiment of the present invention, at least some of the phase changes 310, 320 and 330, are spatial phase changes, effected by applying a time-varying spatial phase change to part of wavefront 300.

[0215] In accordance with another preferred embodiment of the present invention, at least some of the phase changes 310, 320 and 330, are spatial phase changes, effected by applying a non time-varying spatial phase charge to part of wavefront 300, producing spatially phase changed transformed wavefronts 312, 322 and 332, which subsequently produce spatially varying intensity maps 314, 324 and 334 respectively.

[0216] Reference is now made to FIG. 4, which is a simplified functional block diagram illustration of the functionality of FIG. 2, specifically in a case where time-varying, non-spatially varying, spatial phase changes are applied to a transformed wavefront. As seen in FIG. 4, and as explained hereinabove with reference to FIG. 1A, a wavefront 400 is preferably transformed to provide a transformed wavefront 408. A preferred transform is a Fourier transform.

[0217] A first spatial phase change is applied to the transformed wavefront 408 at a first time T1, as indicated by reference numeral 410. This phase change preferably is effected by applying a spatially uniform spatial phase delay D, designated by reference ‘D=D1’, to a given spatial region of the transformed wavefront 408. Thus, at the given spatial region of the transformed wavefront, the value of the phase delay at time T1 is D1, while at the remainder of the transformed wavefront, where no phase delay is applied, the value of the phase delay is D=0.

[0218] The first spatial phase change 410 thereby produces a spatially phase changed transformed wavefront 412 at time T1. This spatially phase changed transformed wavefront 412 is detected, as by detector 158 (FIG. 1B), producing a spatially varying intensity map, an example of which is designated by reference numeral 414, which map is stored as by circuitry 160 (FIG. 1B).

[0219] Thereafter, a second spatial phase change is applied to the transformed wavefront 408 at a second time T2, as indicated by reference numeral 420. This phase change preferably is effected by applying a spatially uniform spatial phase delay D, designated by reference ‘D=D2’, to a given spatial region of the transformed wavefront 408. Thus, at the given spatial region of the transformed wavefront, the value of the phase delay at time T2 is D2, while at the remainder of the transformed wavefront, where no phase delay is applied, the value of the phase delay is D=0.

[0220] The second spatial phase change 420 thereby produces a spatially phase changed transformed wavefront 422 at time T2. This spatially phase changed transformed wavefront 422 is detected, as by detector 158 (FIG. 1B), producing a spatially varying intensity map, an example of which is designated by reference numeral 424, which map is stored as by circuitry 160 (FIG. 1B).

[0221] Thereafter, a third spatial phase change is applied to the transformed wavefront 408 at a third time T3, as indicated by reference numeral 430. This phase change preferably is effected by applying a spatially uniform spatial phase delay D, designated by reference ‘D=D3’, to a given spatial region of the transformed wavefront 408. Thus, at the given spatial region of the transformed wavefront, the value of the phase delay at time T3 is D3, while at the remainder of the transformed wavefront, where no phase delay is applied, the value of the phase delay is D=0.

[0222] The third spatial phase change 430 thereby produces a spatially phase changed transformed wavefront 432 at time T3. This spatially phase changed transformed wavefront 432 is detected, as by detector 158 (FIG. 113), producing a spatially varying intensity map, an example of which is designated by reference numeral 434, which map is stored as by circuitry 160 (FIG. 1B).

[0223] It is appreciated that any suitable number of spatial phase changes may be made at successive times and stored for use in accordance with the present invention.

[0224] In accordance with a preferred embodiment of the present invention, the transform applied to the wavefront 400 is a Fourier transform, thereby providing a Fourier-transformed wavefront 408. In addition, the plurality of phase changed transformed wavefronts 412, 422 and 432 may be further transformed, preferably by a Fourier transform, prior to detection thereof.

[0225] In accordance with a preferred embodiment of the present invention, the spatial region of the transformed wavefront 408 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a spatially central region of the transformed wavefront 408.

[0226] In accordance with an embodiment of the present invention, a phase component comprising relatively high frequency components may be added to the wavefront 400 prior to applying the transform thereto, in order to increase the high-frequency content of the transformed wavefront 408 prior to applying the spatially uniform, spatial phase delays to a spatial region thereof.

[0227] Additionally, in accordance with a preferred embodiment of the present invention, the spatial region of the transformed wavefront 408 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a spatially central region of the transformed wavefront 408, the transform applied to the wavefront 400 is a Fourier transform, and the plurality of phase changed transformed wavefronts 412, 422 and 432 are Fourier transformed prior to detection thereof.

[0228] In accordance with another embodiment of the present invention, the region of the transformed wavefront 408 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a spatially centered generally circular region of the transformed wavefront 408.

[0229] In accordance with yet another embodiment of the present invention, the region of the transformed wavefront 408 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a region covering approximately one half of the entire region in which transformed wavefront 408 is defined.

[0230] In accordance with a preferred embodiment of the present invention, the transformed wavefront 408 includes a non-spatially modulated region, termed a DC region, which represents an image of a light source generating the wavefront 400, and a non-DC region. The region of the transformed wavefront 408 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively includes at least parts of both the DC region and the non-DC region.

[0231] Reference is now made to FIG. 5, which is a simplified functional block diagram illustration of the functionality of FIG. 3, where time-varying, non-spatially varying, spatial phase changes are applied to a wavefront prior to transforming thereof.

[0232] As seen in FIG. 5, a first spatial phase change is applied to a wavefront 500 at a first time T1, as indicated by reference numeral 510. This phase change preferably is effected by applying a spatially uniform spatial phase delay D, designated by reference ‘D=D1’, to a given spatial region of the wavefront 500. Thus, at the given spatial region of the wavefront, the value of the phase delay at time T1 is D1, while at the remainder of the wavefront, where no phase delay is applied, the value of the phase delay is D=0.

[0233] Following application of the first spatial phase change to wavefront 500, a transform, preferably a Fourier transform, is applied thereto, thereby producing a spatially phase changed transformed wavefront 512 at time T1. This spatially phase changed transformed wavefront 512 is detected, as by detector 158 (FIG. 1B), producing a spatially varying intensity map, an example of which is designated by reference numeral 514, which map is stored as by circuitry 160 (FIG. 1B).

[0234] Thereafter, a second spatial phase change is applied to wavefront 500 at a second time T2, as indicated by reference numeral 520. This phase change preferably is effected by applying a spatially uniform spatial phase delay D, designated by reference ‘D=D2’, to a given spatial region of the wavefront 500. Thus, at the given spatial region of the wavefront, the value of the phase delay at time T2 is D2, while at the remainder of the wavefront, where no phase delay is applied, the value of the phase delay is D=0.

[0235] Following application of the second spatial phase change to wavefront 500, a transform, preferably a Fourier transform, is applied thereto, thereby producing a spatially phase changed transformed wavefront 522 at time T2. This spatially phase changed transformed wavefront 522 is detected, as by detector 158 (FIG. 1B), producing a spatially varying intensity map, an example of which is designated by reference numeral 524, which map is stored as by circuitry 160 (FIG. 1B).

[0236] Thereafter, a third spatial phase change is applied to wavefront 500 at a third time T3, as indicated by reference numeral 530. This phase change preferably is effected by applying a spatially uniform spatial phase delay D, designated by reference ‘D=D3’, to a given spatial region of the wavefront 500. Thus, at the given spatial region of the wavefront, the value of the phase delay at time T3 is D3, while at the remainder of the wavefront, where no phase delay is applied, the value of the phase delay is D=0.

[0237] Following application of the third spatial phase change to wavefront 500, a transform, preferably a Fourier transform, is applied thereto, thereby producing a spatially phase changed transformed wavefront 532 at time T3. This spatially phase changed transformed wavefront 532 is detected, as by detector 158 (FIG. 1B), producing a spatially varying intensity map, an example of which is designated by reference numeral 534, which map is stored as by circuitry 160 (FIG. 1B).

[0238] It is appreciated that any suitable number of spatial phase changes may be made at successive times and stored for use in accordance with the present invention.

[0239] In accordance with a preferred embodiment of the present invention, the spatial region of the wavefront 500 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a spatially central region of the wavefront 500.

[0240] In accordance with an embodiment of the present invention, a phase component comprising relatively high frequency components may be added to the wavefront 500 prior to applying the spatial phase changes thereto, in order to increase the high-frequency content of the wavefront 500.

[0241] Additionally, in accordance with a preferred embodiment of the present invention, the spatial region of the wavefront 500 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a spatially central region of the wavefront 500, the transforms are Fourier transforms, and the plurality of phase changed transformed wavefronts 512, 522 and 532 are Fourier transformed prior to detection thereof.

[0242] In accordance with another embodiment of the present invention, the region of the wavefront 500 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a spatially centered generally circular region of the wavefront 500.

[0243] In accordance with yet another embodiment of the present invention, the region of the wavefront 500 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively is a region covering approximately one half of the entire region in which wavefront 500 is defined.

[0244] In accordance with a preferred embodiment of the present invention, the wavefront 500 includes a non-spatially modulated region, termed a DC region, which represents an image of a light source generating the wavefront 500, and a non-DC region. The region of the wavefront 500 to which the spatially uniform, spatial phase delays D1, D2 and D3 are applied at times T1, T2 and T3 respectively includes at least parts of both the DC region and the non-DC region.

[0245] Reference is now made to FIG. 6, which is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different wavelength components of a transformed wavefront. As seen in FIG. 6, a wavefront 600, which comprises a plurality of different wavelength components, is preferably transformed to obtain a transformed wavefront 602. The transform is preferably a Fourier transform.

[0246] Similarly to wavefront 600, the transformed wavefront 602 also includes a plurality of different wavelength components, represented by reference numerals 604, 606 and 608. It is appreciated that both the wavefront 600 and the transformed wavefront 602 can include any suitable number of wavelength components.

[0247] A plurality of phase changes, preferably spatial phase changes, represented by reference numerals 610, 612 and 614 are applied to respective wavelength components 604, 606 and 608 of the transformed wavefront, thereby providing a plurality of differently phase changed transformed wavefront components, represented by reference numerals 620, 622 and 624 respectively.

[0248] The phase changed transformed wavefront components 620, 622, and 624 may be transformed, preferably by a Fourier transform, and are subsequently detected, as by detector 158 (FIG. 1B), producing spatially varying intensity maps, examples of which are designated by reference numerals 630, 632 and 634 respectively. These intensity maps are subsequently stored as by circuitry 160 (FIG. 1B).

[0249] In accordance with an embodiment of the present invention, phase changes 610, 612 and 614 are effected by passing the transformed wavefront 602 through an object, at least one of whose thickness and refractive index varies spatially, thereby applying a different spatial phase delay to each of the wavelength components 604, 606 and 608 of the transformed wavefront.

[0250] In accordance with another embodiment of the present invention, the phase changes 610, 612 and 614 are effected by reflecting the transformed wavefront 602 from a spatially varying surface, thereby applying a different spatial phase delay to each of the wavelength components 604, 606 and 608 of the transformed wavefront.

[0251] In accordance with yet another embodiment of the present invention, the phase changes 610, 612 and 614 are realized by passing the transformed wavefront 602 through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially. The spatial variance of the thickness or of the refractive index of the plurality of objects is selected in a way such that the phase changes 610, 612 and 614 differ to a selected predetermined extent for at least some of the plurality of different wavelength components 604, 606 and 608.

[0252] Alternatively, the spatial variance of the thickness or refractive index of the plurality of objects is selected in a way such that the phase changes 610, 612 and 614 are identical for at least some of the plurality of different wavelength components 604, 606 and 608.

[0253] Additionally, in accordance with an embodiment of the present invention, the phase changes 610, 612 and 614 are time-varying spatial phase changes. In such a case, the plurality of phase changed transformed wavefront components 620, 622 and 624 include a plurality of differently phase changed transformed wavefronts for each wavelength component thereof, and the intensity maps 630, 632 and 634 include a time-varying intensity map for each such wavelength component.

[0254] In accordance with an embodiment of the present invention, termed a “white light” embodiment, all the wavelength components may be detected by a single detector, resulting in a time-varying intensity map representing several wavelength components.

[0255] In accordance with another embodiment of the present invention, the plurality of phase changed transformed wavefront components 620, 622 and 624 are broken down into separate wavelength components, such as by a spatial separation effected, for example, by passing the phase changed transformed wavefront components through a dispersion element. In such a case, the intensity maps 630, 632 and 634 are provided simultaneously for all of the plurality of different wavelength components.

[0256] Reference is now made to FIG. 7, which is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different wavelength components of a wavefront, prior to transforming thereof. As seen in FIG. 7, a wavefront 700 comprises a plurality of different wavelength components 704, 706 and 708. It is appreciated that the wavefront can include any suitable number of wavelength components.

[0257] A plurality of phase changes, preferably spatial phase changes, represented by reference numerals 710, 712 and 714, are applied to the respective wavelength components 704, 706 and 708 of the wavefront.

[0258] Following application of the spatial phase changes to wavefront components 704, 706 and 708, a transform, preferably a Fourier transform, is applied thereto, thereby providing a plurality of different phase changed transformed wavefront components, represented by reference numerals 720, 722 and 724 respectively.

[0259] These phase changed transformed wavefront components 720, 722 and 724 are subsequently detected, as by detector 158 (FIG. 113), producing spatially varying intensity maps, examples of which are designated by reference numerals 730, 732 and 734. These intensity maps are subsequently stored as by circuitry 160 (FIG. 1B).

[0260] In accordance with an embodiment of the present invention, phase changes 710, 712 and 714 are effected by passing the wavefront 700 through an object, at least one of whose thickness and refractive index varies spatially, thereby applying a different spatial phase delay to each of the wavelength components 704, 706 and 708 of the wavefront.

[0261] In accordance with another embodiment of the present invention, the phase changes 710, 712 and 714 are effected by reflecting the wavefront 700 from a spatially varying surface, thereby applying a different spatial phase delay to each of the wavelength components 704, 706 and 708 of the wavefront.

[0262] In accordance with yet another embodiment of the present invention phase changes 710, 712 and 714 are realized by passing the wavefront 700 through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially. The spatial variance of the thickness or refractive index of these objects is selected in a way such that the phase changes 710, 712 and 714 differ to a selected predetermined extent for at least some of the plurality of different wavelength components 704, 706 and 708.

[0263] Alternatively, the spatial variance of the thickness or refractive index of these objects is selected in a way that the phase changes 710, 712 and 714 are identical for at least some of the plurality of different wavelength components 704, 706 and 708.

[0264] Reference is now made to FIG. 8, which is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different polarization components of a transformed wavefront. As seen in FIG. 8, a wavefront 800, which comprises a plurality of different polarization components, is preferably transformed to obtain a transformed wavefront 802. The transform is preferably a Fourier transform. Similarly to wavefront 800, the transformed wavefront 802 also includes a plurality of different polarization components, represented by reference numerals 804 and 806. It is appreciated that the polarization components 804 and 806 can be either spatially different or spatially identical, but are each of different polarization. It is further appreciated that both the wavefront 800 and the transformed wavefront 802 preferably each include two polarization components but can include any suitable number of polarization components.

[0265] A plurality of phase changes, preferably spatial phase changes, represented by reference numerals 810 and 812, are applied to the respective polarization components 804 and 806 of the transformed wavefront 802, thereby providing a plurality of differently phase changed transformed wavefront components, represented by reference numerals 820 and 822 respectively.

[0266] It is appreciated that phase changes 810 and 812 can be different for at least some of the plurality of different polarization components 804 and 806. Alternatively, phase changes 810 and 812 can be identical for at least some of the plurality of different polarization components 804 and 806.

[0267] The phase changed transformed wavefront components 820 and 822 are detected, as by detector 158 (FIG. 1B), producing spatially varying intensity maps, examples of which are designated by reference numerals 830 and 832. These intensity maps are subsequently stored as by circuitry 160 (FIG. 1B).

[0268] Reference is now made to FIG. 9, which is a simplified functional block diagram illustration of the functionality of FIG. 1A where phase changes are applied to a plurality of different polarization components of a wavefront prior to transforming thereof. As seen in FIG. 9, a wavefront 900 comprises a plurality of different polarization components 904 and 906. It is appreciated that the wavefront preferably includes two polarization components but can include any suitable number of polarization components.

[0269] A plurality of phase changes, preferably spatial phase changes, represented by reference numerals 910 and 912, are applied to the respective polarization components 904 and 906 of the wavefront.

[0270] It is appreciated that phase changes 910 and 912 can be different for at least some of the plurality of different polarization components 904 and 906. Alternatively, phase changes 910 and 912 can be set to be identical for at least some of the plurality of different polarization components 904 and 906.

[0271] Following application of the spatial phase changes to wavefront components 904 and 906, a transform, preferably a Fourier transform, is applied thereto, thereby providing a plurality of different phase changed transformed wavefront components, designated by reference numerals 920 and 922 respectively.

[0272] Phase changed transformed wavefront components 920 and 922 are subsequently detected, as by detector 158 (FIG. 1B), producing spatially varying intensity maps, examples of which are designated by reference numeral 930 and 932. These intensity maps are subsequently stored as by circuitry 160 (FIG. 1B).

[0273] Reference is now made to FIG. 10A, which is a simplified functional block diagram illustration of the functionality of FIG. 1A, where a wavefront being analyzed comprises at least one one-dimensional component. In the embodiment of FIG. 10A, a one-dimensional Fourier transform is applied to the wavefront. Preferably, the transform is performed in a dimension perpendicular to a direction of propagation of the wavefront being analyzed, thereby to obtain at least one one-dimensional component of the transformed wavefront in the dimension perpendicular to the direction of propagation.

[0274] A plurality of different phase changes are applied to each of the at least one one-dimensional components, thereby obtaining at least one one-dimensional component of the plurality of phase changed transformed wavefronts.

[0275] A plurality of intensity maps are employed to obtain an output indicating amplitude and phase of the at least one one-dimensional component of the wavefront being analyzed.

[0276] As seen in FIG. 10A, a plurality of different phase changes are applied to at least one one-dimensional component of a transformed wavefront. In the illustrated embodiment, typically five one-dimensional components of a wavefront are shown and designated by reference numerals 1001, 1002, 1003, 1004 and 1005. The wavefront is transformed, preferably by a Fourier transform. It is thus appreciated that due to transform of the wavefront, the five one-dimensional components 1001, 1002, 1003, 1004 and 1005 are transformed into five corresponding one-dimensional components of the transformed wavefront, respectively designated by reference numerals 1006, 1007, 1008, 1009 and 1010.

[0277] Three phase changes, respectively designated 1011, 1012 & 1013 are each applied to the one-dimensional components 1006, 1007, 1008, 1009 and 1010 of transformed wavefront to produce three phase changed transformed wavefronts, designated generally by reference numerals 1016, 1018 and 1020.

[0278] In the illustrated embodiment, phase changed transformed wavefront 1016 includes five one-dimensional components, respectively designated by reference numerals 1021, 1022, 1023, 1024 and 1025.

[0279] In the illustrated embodiment, phase changed transformed wavefront 1018 includes five one-dimensional components, respectively designated by reference numerals 1031, 1032, 1033, 1034 and 1035.

[0280] In the illustrated embodiment, phase changed transformed wavefront 1020 includes five one-dimensional components, respectively designated by reference numerals 1041, 1042, 1043, 1044 and 1045.

[0281] The phase changed transformed wavefronts 1016, 1018 and 1020 are detected, as by detector 158 (FIG. 1B), producing three intensity maps, designated generally by reference numerals 1046, 1048 and 1050.

[0282] In the illustrated embodiment, intensity map 1046 includes five one-dimensional intensity map components, respectively designated by reference numerals 1051, 1052, 1053, 1054 and 1055.

[0283] In the illustrated embodiment, intensity map 1048 includes five one-dimensional intensity map components, respectively designated by reference numerals 1061, 1062, 1063, 1064 and 1065.

[0284] In the illustrated embodiment, intensity map 1050 includes five one-dimensional intensity map components, respectively designated by reference numerals 1071, 1072, 1073, 1074 and 1075.

[0285] The intensity maps 1046, 1048 and 1050 are stored as by circuitry 160 (FIG. 1B).

[0286] In accordance with an embodiment of the present invention, the wavefront being analyzed, illustrated in FIG. 10A by the one-dimensional components 1001, 1002, 1003, 1004 and 1005, may comprise a plurality of different wavelength components and the plurality of different phase changes, 1011, 1012 and 1013, are applied to the plurality of different wavelength components of each of the plurality of one-dimensional components of the wavefront being analyzed. Preferably, obtaining a plurality of intensity maps 1046, 1048 and 1050, includes dividing the plurality of one-dimensional components of the plurality of phase changed transformed wavefronts 1016, 1018 and 1020 into separate wavelength components.

[0287] Preferably, dividing the plurality of one-dimensional components of the plurality of phase changed transformed wavefronts into separate wavelength components is achieved by passing the plurality of phase changed transformed wavefronts 1016, 1018 and 1020 through a dispersion element.

[0288] Reference is now made to FIG. 10B, which is a simplified partially schematic, partially pictorial illustration of a wavefront analysis system suitable for carrying out the functionality of FIG. 10A in accordance with a preferred embodiment of the present invention.

[0289] As seen in FIG. 10B, a wavefront, here designated by reference numeral 1080, and here including five one-dimensional components 1081, 1082, 1083, 1084 and 1085 is focused, as by a cylindrical lens 1086 onto a single axis displaceable phase manipulator 1087, which is preferably located at the focal plane of lens 1086. Lens 1086 preferably produces a one-dimensional Fourier transform of each of the one-dimensional wavefront components 1081, 1082, 1083, 1084 and 1085 along the Y-axis.

[0290] As seen in FIG. 10B, the phase manipulator 1087 preferably comprises a multiple local phase delay element, such as a spatially non-uniform transparent object, typically including five different phase delay regions, each arranged to apply a phase delay to one of the one-dimensional components at a given position of the object along an axis, here designated as the X-axis, extending perpendicularly to the direction of propagation of the wavefront along a Z-axis and perpendicular to the axis of the transform produced by lens 1086, here designated as the Y-axis.

[0291] A second lens 1088, preferably a cylindrical lens, is arranged so as to image the one-dimensional components 1081, 1082, 1083, 1084 and 1085 onto a detector 1089, such as a CCD detector. Preferably the second lens 1088 is arranged such that the detector 1089 lies in its focal plane. The output of detector 1089 is preferably supplied to data storage and processing circuitry 1090, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A.

[0292] There is provided relative movement between the optical system comprising phase manipulator 1087, lenses 1086 and 1088 and detector 1089 and the one-dimensional wavefront components 1081, 1082, 1083, 1084 and 1085 along the X-axis. This relative movement sequentially matches different phase delay regions with different wavefront components, such that preferably each wavefront component passes through each phase delay region of the phase manipulator 1087.

[0293] It is a particular feature of the embodiment of FIGS. 10A and 10B, that each of the one dimensional components of the wavefront is separately processed. Thus, in the context of FIG. 10B, it can be seen that the five one-dimensional wavefront components 1081, 1082, 1083, 1084 and 1085 are each focused by a separate portion of the cylindrical lens 1086, are each imaged by a corresponding separate portion of the cylindrical lens 1088 and each pass through a distinct region of the phase manipulator 1087. The images of each of the five one-dimensional wavefront components 1081, 1082, 1083, 1084 and 1085 at detector 1089 are thus seen to be separate and distinct images, as designated respectively by reference numerals 1091, 1092, 1093, 1094 and 1095. It is appreciated that these images may appear on separate detectors together constituting detector 1089 instead of on a monolithic detector.

[0294] In accordance with an embodiment of the present invention, the transform applied to the wavefront includes an additional Fourier transform. This additional Fourier transform may be performed by lens 1086 or by an additional lens and is operative to minimize cross-talk between different one-dimensional components of the wavefront. In such a case, preferably a further transform is applied to the phase changed transformed wavefront. This further transform may be performed by lens 1088 or by an additional lens.

[0295] Reference is now made to FIG. 11, which is a simplified functional block diagram illustration of the functionality of FIG. 1A, where an additional transform is applied following the application of spatial phase changes. As seen in FIG. 11, and as explained hereinabove with reference to FIG. 1A, a wavefront 1100 is transformed, preferably by a Fourier transform and a plurality of phase changes are applied to the transformed wavefront, thereby to provide a plurality of differently phased changed transformed wavefronts, represented by reference numerals 1120, 1122, and 1124.

[0296] The phase changed transformed wavefronts are subsequently transformed, preferably by a Fourier transform, and then detected, as by detector 158 (FIG. 1B), producing spatially varying intensity maps, examples of which are designated by reference numerals 1130, 1132 and 1134. These intensity maps are subsequently stored as by circuitry 160 (FIG. 1B).

[0297] It is appreciated that any suitable number of differently phased changed transformed wavefronts can be obtained, and subsequently transformed to a corresponding plurality of intensity maps to be stored for use in accordance with the present invention.

[0298] Reference in now made to FIG. 12, which 1s a simplified functional block diagram illustration of the functionality of FIG. 1A, wherein intensity maps are employed to provide information about a wavefront being analyzed, such as indications of amplitude and phase of the wavefront. As seen in FIG. 12, and as explained hereinabove with reference to FIG. 1A, a wavefront 1200 is transformed, preferably by a Fourier transform, and phase changed by a phase-change function to obtain several, preferably at least three, differently phase-changed transformed wavefronts, respectively designated by reference numerals 1210, 1212 and 1214. The phase changed transformed wavefronts 1210, 1212 and 1214 are subsequently detected, as by detector 158 (FIG. 1B), producing spatially varying intensity maps, examples of which are designated by reference numerals 1220, 1222 and 1224.

[0299] In parallel to producing the plurality of intensity maps, such as intensity maps 1220, 1222 and 1224, the expected intensity maps are expressed as a first function of the amplitude of wavefront 1200, of the phase of wavefront 1200, and of the phase change function characterizing the differently phase changed transformed wavefronts 1210, 1212 and 1214, as indicated at reference numeral 1230.

[0300] In accordance with a preferred embodiment of the present invention, at least one of the phase and the amplitude of the wavefront is unknown or both the phase and the amplitude are unknown. The phase-change function is known.

[0301] The first function of the phase and amplitude of the wavefront and of the phase change function is subsequently solved as indicated at reference numeral 1235, such as by means of a computer 136 (FIG. 1A), resulting in an expression of at least one and possibly both of the amplitude and phase of wavefront 1200 as a second function of the intensity maps 1220, 1222 and 1224, as indicated at reference numeral 1240.

[0302] The second function is then processed together with the intensity maps 1220, 1222 and 1224 as indicated at reference numeral 1242. As part of this processing, detected intensity maps 1220, 1222 and 1224 are substituted into the second function. The processing may be carried out by means of a computer 136 (FIG. 1A) and provides information regarding wavefront 1200, such as indications of at least one and possibly both of the amplitude and the phase of the wavefront.

[0303] In accordance with a further embodiment of the present invention, the plurality of intensity maps comprises at least four intensity maps. In such a case, employing the plurality of intensity maps to obtain an indication of at least one of the phase and the amplitude of the wavefront 1200 includes employing a plurality of combinations, each of the combinations being a combination of at least three of the plurality of intensity maps, to provide a plurality of indications of at least one of the phase and the amplitude of wavefront 1200. Preferably, this methodology also includes employing the plurality of indications of at least one of the phase and the amplitude of the wavefront 1200 to provide an enhanced indication at least one of the phase and the amplitude of the wavefront 1200.

[0304] In accordance with a preferred embodiment of the present invention, at least some of the plurality of indications of the amplitude and phase are at least second order indications of the amplitude and phase of the wavefront 1200.

[0305] In accordance with another embodiment of the present invention, the first function may be solved as a function of some unknowns to obtain the second function by expressing, as indicated by reference numeral 1240, some unknowns, such as at least one of the amplitude and phase of wavefront 1200, as a second function of the intensity maps.

[0306] Accordingly, solving the first function may include:

[0307] defining a complex function of the amplitude of wavefront 1200, of the phase of wavefront 1200, and of the phase change function characterizing the differently phase changed transformed wavefronts 1210, 1212 and 1214. This complex function is characterized in that intensity at each location in the plurality of intensity maps is a function predominantly of a value of the complex function at that location and of the amplitude and the phase of wavefront 1200 at the same location;

[0308] expressing the complex function as a third function of the plurality of intensity maps 1220, 1222 and 1224; and

[0309] obtaining values for the unknowns, such as at least one of phase and amplitude of wavefront 1200, by employing the complex function expressed as a function of the plurality of intensity maps.

[0310] In accordance with this embodiment, preferably the complex function is a convolution of another complex function, which has an amplitude and phase identical to the amplitude and phase of wavefront 1200, and of a Fourier transform of the phase change function characterizing the differently phase changed transformed wavefronts 1210, 1212 and 1214.

[0311] Reference in now made to FIG. 13, which is a simplified functional block diagram illustration of part of the functionality of FIG. 1A, wherein the transform applied to the wavefront being analyzed is a Fourier transform, wherein at least three different spatial phase changes are applied to the thus transformed wavefront, and wherein at least three intensity maps are employed to obtain indications of at least one of the phase and the amplitude of the wavefront.

[0312] As explained hereinabove with reference to FIG. 1A, a wavefront 100 (FIG. 1A) being analyzed, is transformed and phase changed by at least three different spatial phase changes, all governed by a spatial function, to obtain at least three differently phase-changed transformed wavefronts, represented by reference numerals 120, 122 and 124 (FIG. 1A) which are subsequently detected, as by detector 158 (FIG. 1B), producing spatially varying intensity maps, examples of which are designated by reference numerals 130, 132 and 134 (FIG. 1A). As seen in FIG. 13, and designated as sub-functionality “C” hereinabove with reference in FIG. 1A, the intensity maps are employed to obtain an output indication of at least one and possibly both of the phase and the amplitude of the wavefront being analyzed.

[0313] Turning to FIG. 13, it is seen that the wavefront being analyzed is expressed as a first complex function f(x)=A(x)_(e) ^(iφ(x)), where ‘x’ is a general indication of a spatial location. The complex function has an amplitude distribution A(x) and a phase distribution φ(x) identical to the amplitude and phase of the wavefront being analyzed. The first complex function f(x)=A(x)_(eiφ(x)) is indicated by reference numeral 1300.

[0314] As noted hereinabove with reference to FIG. 1A, each of the plurality of different spatial phase changes is applied to the transformed wavefront preferably by applying a spatially uniform spatial phase delay having a known value to a given spatial region of the transformed wavefront. As seen in FIG. 13, the spatial function governing these different phase changes is designated by ‘G’ and an example of which, for a phase delay value of θ, is designated by reference numeral 1304.

[0315] Function ‘G’ is a spatial function of the phase change applied in each spatial location of the transformed wavefront. In the specific example designated by reference numeral 1304, the spatially uniform spatial phase delay, having a value of θ, is applied to a spatially central region of the transformed wavefront, as indicated by the central part of the function having a value of θ, which is greater than the value of the function elsewhere.

[0316] A plurality of expected intensity maps, indicated by spatial functions I₁(x), I₂(x) and I₃(x), are each expressed as a function of the first complex function f(x) and of the spatial function G, as indicated by reference numeral 1308.

[0317] Subsequently, a second complex function S(x), which has an absolute value |S(x)| and a phase α(x), is defined as a convolution of the first complex function f(x) and of a Fourier transform of the spatial function ‘G’. This second complex function, designated by reference numeral 1312, is indicated by the equation S(x)=f(x)*{tilde over (ℑ)}(G)=|S(x)|_(e) ^(iα(x)), where the symbol ‘*’ indicates convolution and {tilde over (ℑ)}(G) is the Fourier transform of the function ‘G’.

[0318] The difference between φ(x), the phase of the wavefront, and α(x), the phase of the second complex function, is indicated by ψ(x), as designated by reference numeral 1316.

[0319] The expression of each of the expected intensity maps as a function off(x) and G, as indicated by reference numeral 1308, the definition of the absolute value and the phase of S(x), as indicated by reference numeral 1312 and the definition of ψ(x), as indicated by reference numeral 1316, enables expression of each of the expected intensity maps as a third function of the amplitude of the wavefront A(x), the absolute value of the second complex function ″S(x)|, the difference between the phase of the wavefront and the phase of the second complex function ψ(x), and the known phase delay produced by one of the at least three different phase changes which each correspond to one of the at least three intensity maps.

[0320] This third function is designated by reference numeral 1320 and includes three functions, each preferably having the general form I_(n)(x) = A(x) + (^(θ_(n)) − 1)S(x)^(−  ψ(x))²

[0321] where I_(n)(x) are the expected intensity maps and n=1, 2 or 3. In the three functions, θ₁, θ₂ and θ₃ are the known values of the uniform spatial phase delays, each applied to a spatial region of the transformed wavefront, thus effecting the plurality of different spatial phase changes which produce the intensity maps I₁(x), I₂(x) and I₃(x), respectively.

[0322] It is appreciated that preferably the third function at any given spatial location x₀ is a function of A, ψ and |S| only at the same spatial location x₀.

[0323] The intensity maps are designated by reference numeral 1324.

[0324] The third function is solved for each of the specific spatial locations x₀, by solving at least three equations, relating to at least three intensity values I₁(x₀), I₂(x₀) and I₃(x₀) at at least three different phase delays θ₁, θ₂ and θ₃, thereby to obtain at least part of three unknowns A(x₀), ″S(x₀)| and ψ(x₀). This process is typically repeated for all spatial locations and results in obtaining the amplitude of the wavefront A(x), the absolute value of the second complex function |S(x)| and the difference between the phase of the wavefront and the phase of the second complex function ψ(x), as indicated by reference numeral 1328.

[0325] Thereafter, once A(x), |S(x)| and ψ(x) are known, the equation defining the second complex function, represented by reference numeral 1312, is typically solved globally for a substantial number of spatial locations ‘x’ to obtain α(x), the phase of the second complex function, as designated by reference numeral 1332.

[0326] Finally, the phase φ(x) of the wavefront being analyzed is obtained by adding the phase α(x) of the second complex function to the difference ψ(x) between the phase of the wavefront and the phase of the second complex function, as indicated by reference numeral 1336.

[0327] In accordance with an embodiment of the present invention, the absolute value |S| of the second complex function is obtained preferably for every specific spatial location x₀ by approximating the absolute value to a polynomial of a given degree in the spatial location x.

[0328] In accordance with another preferred embodiment of the present invention, the phase α(x) of the second complex function is obtained by expressing the second complex function S(x) as an eigen-value problem, such as S=S·M where M is a matrix, and the complex function is an eigen-vector of the matrix obtained by an iterative process. An example of such an iterative process is S₀=|S|, S_(n+1)=S_(n)M /∥S_(n)M∥,where n is the iterative step number.

[0329] In accordance with yet another preferred embodiment of the present invention, the phase α(x) of the second complex function is obtained by approximating the Fourier transform of the spatial function ‘G’, governing the spatial phase change, to a polynomial in the location x, by approximating the second complex function S(x) to a polynomial in the location x, and by solving, according to these approximations, the equation defining the second complex function: ${{S(x)} = {\left( {\frac{{A(x)}^{{\psi}{(x)}}}{{S(x)}}{S(x)}} \right)*{\lbrack G\rbrack}}},$

[0330] where the function $\frac{{A(x)}^{{\psi}{(x)}}}{{S(x)}}$

[0331] is known.

[0332] In accordance with still another preferred embodiment of the present invention, at any location x the amplitude A(x) of the wavefront being analyzed, the absolute value |S(x)| of the second complex function, and the difference ψ(x) between the phase of the second complex function and the phase of the wavefront are obtained by a best-fit method, such as a least-square method, preferably a linear least-square method, from the values of the intensity maps at this location I_(n)(x), where n=1,2, . . . , N and N is the number of intensity maps. The accuracy of this process increases as the number N of the plurality of intensity maps increases.

[0333] In accordance with one preferred embodiment of the present invention, the plurality of different phase changes comprises at least four different phase changes, the plurality of intensity maps comprises at least four intensity maps, and the function designated by reference numeral 1320 can express each of the expected intensity maps as a third function of:

[0334] the amplitude of the wavefront A(x);

[0335] the absolute value of the second complex function |S(x)|;

[0336] the difference between the phase of the wavefront and the phase of the second complex function ψ(x);

[0337] the known phase delay produced by one of the at least four different phase changes each of which corresponds to one of the at least four intensity maps; and

[0338] at least one additional unknown relating to the wavefront analysis,

[0339] where the number of the at least one additional unknown is no greater than the number by which the plurality intensity maps exceeds three.

[0340] The third function 1320, is then solved by solving at least four equations, resulting from at least four intensity values at at least four different phase delays, thereby to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function, the difference between the phase of the wavefront and the phase of the second complex function and the at least one additional unknown.

[0341] In accordance with another preferred embodiment of the present invention, the values of the uniform spatial phase delays θ₁, θ₂, . . . , θ_(N) applied to a spatial region of the transformed wavefront, thus effecting the plurality of different spatial phase changes, producing the intensity maps I₁(x), I₂(x), . . . , I_(N)(X) respectively, are chosen as to maximize contrast in the intensity maps and to minimize effects of noise on the phase of the wavefront being analyzed.

[0342] In accordance with one more preferred embodiment of the present invention, the function designated by reference numeral 1320, expressing each of the expected intensity maps as a third function of the amplitude of the wavefront A(x), the absolute value of the second complex function |S(x)|, the difference between the phase of the wavefront and the phase of the second complex function ψ(x), and the known phase delay θ₁ produced by one of the at least three different phase changes which each correspond to one of the at least three intensity maps, comprises several functionalities:

[0343] defining fourth, fifth and sixth complex functions, designated as β₀(x), β_(s)(x) and β_(c)(x) respectively, none of which is a function of any of the plurality of intensity maps or of the spatial function ‘G’ governing the phase change. Each of the fourth, fifth and sixth complex functions is preferably a function of the amplitude of the wavefront A(x), the absolute value of the second complex function |S(x)|, the difference between the phase of the wavefront and the phase of the second complex function ψ(x); and

[0344] expressing each of the plurality of intensity maps I_(n)(x) as I_(n)(x)=β₀(x)+β_(c)(x)cos(θ_(n))+β_(s)(x)sin(θ_(n)), where θ_(n) is the value of the phase delay corresponding to intensity map I_(n)(x). Each intensity map I_(n)(x), where n 1,2, . . . N, preferably expressed as I_(n)(x) = A(x) + (^(θ_(n)) − 1)S(x)^(−ψ(x))²,

[0345] can be subsequently expressed as I_(n)(x)=β₀(x)+β_(c)(x)cos(θ_(n))+β_(s)(x)sin(θ_(n)), where $\left\{ \begin{matrix} {{\beta_{0}(x)} = {{A(x)}^{2} + {2{{S(x)}}^{2}} - {2{A(x)}{{S(x)}}{\cos (\psi)}}}} \\ {{\beta_{C}(x)} = {{2{A(x)}{{S(x)}}{\cos (\psi)}} - {2{{S(x)}}^{2}}}} \\ {{\beta_{S}(x)} = {2{A(x)}{{S(x)}}{\sin (\psi)}}} \end{matrix}\quad \right.$

[0346] Preferably the foregoing methodology also includes solving the third function 1320 by using a linear least-square method to compute from the different intensities I(θ₁), . . . , I(θ_(N)), the values of β₀, β_(c) and β_(s) best fitting to I(θ_(n))=β₀+β_(c) cos θ_(n)+β_(s) sin θ_(n). Subsequently the amplitude A(x) is found by

[0347] the absolute value |S(x)| of the second complex function is found by solving the second degree equation ${{{S(x)}}^{4} - {{\beta_{0}(x)}{{S(x)}}^{2}} + \frac{{\beta_{c}(x)}^{2} + {\beta_{s}(x)}^{2}}{4}} = 0$

[0348] for |S(x)|², and ψ(x) is found by ψv(x)=arg(β_(c)(x)+2|S(x)|²+iβ_(s)(x)).

[0349] In accordance with yet another preferred embodiment of the present invention, solving of the third function, designated by reference numeral 1320, to obtain, as designated by reference numeral 1328, the amplitude of the wavefront A(x), the absolute value of the second complex function |S(x)| and the difference between the phase of the wavefront and the phase of the second complex function ψ(x), includes several functionalities:

[0350] obtaining two solutions for the absolute value |S(x)| of the second complex function, these two solutions, being designated by |S_(h)(x)| and |S_(l)(x)|, namely a higher value solution and a lower value solution respectively; and

[0351] combining the two solutions into an enhanced absolute value solution |S(x)| for the absolute value of the second complex function, by choosing at each spatial location ‘x₀’ either the higher value solution |S_(h)(x₀)| or the lower value solution |S_(l)(x₀)| such that the enhanced absolute value solution satisfies the second complex function, designated by reference numeral 1312.

[0352] Preferably the methodology also includes:

[0353] obtaining two solutions for each of the amplitude A(x) of the wavefront being analyzed and the difference ψ(x) between the phase of the wavefront and the phase of the second complex function, these two solutions being higher value solutions A_(h)(X) and ψ_(h)(x) and lower value solutions A_(l)(x) and ψ_(l)(x); and

[0354] combining the two solutions A_(h)(x) and A_(l)(x) for the amplitude into an enhanced amplitude solution A(x) by choosing at each spatial location ‘x₀’ either the higher value solution A_(h)(x₀) or the lower value solution A_(l)(x₀) in a way that at each spatial location ‘x₀’ if |S_(h)(x₀)| is chosen for the absolute value solution, then A_(h)(X₀) is chosen for the amplitude solution and at each location ‘x₁’ if |S_(l)(x)| is chosen for the absolute value solution, then A_(l)(x₁) is chosen for the amplitude solution; and

[0355] combining the two solutions ψ_(h)(x) and ψ_(l)(x) of the difference between the phase of the wavefront and the phase of the second complex function into an enhanced difference solution ψ(x), by choosing at each spatial location ‘x₀’ either the higher value solution ψ_(h)(x₀) or the lower value solution ψ_(l)(x₀) in a way that at each spatial location ‘x₀’ if |S_(h)(x₀)| is chosen for the absolute value solution, then ψ_(h)(x₀) is chosen for the difference solution and at each location ‘x₁’ if |S_(l)(x₁) is chosen for the absolute value solution, then ψ_(l)(x₁) is chosen for the difference solution.

[0356] Additionally, in accordance with an embodiment of the present invention, the plurality of different phase changes applied to the transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts, also include amplitude changes, resulting in a plurality of differently phase and amplitude changed transformed wavefronts. These amplitude changes are preferably known amplitude attenuations applied to the same spatial region of the transformed wavefront to which the uniform phase delays θ₁, θ₂, . . . , θ_(N), are applied, the spatial region being defined by the spatial function ‘G’.

[0357] The amplitude attenuations are designated by σ₁, σ₂, . . . , σ_(N) where the n-th change, where n=1, 2, . . . N, applied to the transformed wavefront includes a phase change θ_(n) and an amplitude attenuation σ_(n). It is appreciated that some of the phase changes may be equal to zero, indicating no phase-change and that some of the amplitude attenuations may be equal to unity, indicating no amplitude attenuation.

[0358] In this embodiment, the function designated by reference numeral 1320, expressing each of the expected intensity maps I_(n)(x) as a third function of the amplitude of the wavefront A(x), the absolute value of the second complex function |S(x)|, the difference between the phase of the wavefront and the phase of the second complex function ψ(x), and the phase delay θ_(n), also expresses each of the expected intensity maps also as a function of the amplitude attenuation σ_(n) and comprises several functionalities:

[0359] defining fourth, fifth, sixth and seventh complex functions, designated by β₀(x), β₁(x), β₂(x) and β₃(x) respectively, none of which is a function of any of the plurality of intensity maps or of the spatial function ‘G’ governing the phase and amplitude changes. Each of the fourth, fifth, sixth and seventh complex functions is preferably a function of the amplitude of the wavefront A(x), the absolute value of the second complex function |S(x)|, the difference between the phase of the wavefront and the phase of the second complex function ψ(x);

[0360] defining an eighth function, designated μ, as a combination of the phase delay and of the amplitude attenuation, where for the n-th change applied to the transformed wavefront, including a phase change θ_(n) and an amplitude attenuation σ_(n), this eighth function is designated by μ_(n). Preferably the combination μ_(n) is defined by μ_(n)=σ_(n)e^(lθ) ^(_(n)) −1; and

[0361] expressing each of the plurality of intensity maps I_(n)(x) as I_(n)(x)=β₀(X)+β₁(x)|μ_(n)|²+β₂(x)μ_(n)+β₃(x){overscore (μ)}_(n), where β₀(x)=A²(x); β₁(x)=|S(x)²|; β₂ (x)=A(x)⊕S(x)|e^(−lψ(x)) and β₃(x)=A(x)|S(x)|e^(iψ(x)).

[0362] Preferably the foregoing methodology also includes solving the third function 1320 by computing from the different intensities I_(n)(x), the values of β₀(x), β₁(x), β₂(x) and β₃(x) best fitting to the equation I_(n)(x)=β₀(x)+β₁(x)|μ_(n)|²+β₂(x)μ_(n)β₃(x){overscore (μ)}_(n). Subsequently the amplitude A(x) is found by A(x)={square root}{square root over (β₀(x))}, the absolute value |S(x)| of the second complex function is found by |S(x)|={square root}{square root over (β₁(x))} and ψ(x) is found by solving e^(iψ(x))=angle(β₃(x)).

[0363] It is appreciated that the amplitude attenuations σ₁, σ₂, . . . , σ_(N), may be unknown. In such a case, additional intensity maps are obtained, where the number of the unknowns is no greater than the number by which the plurality of intensity maps exceeds three. The unknowns are obtained in a manner similar to that described hereinabove, where there exists at least one unknown relating to the wavefront analysis.

[0364] Reference is now made to FIG. 14, which is a simplified partially schematic, partially pictorial illustration of part of one preferred embodiment of a wavefront analysis system of the type shown in FIG. 1B. As seen in FIG. 14, a wavefront, here designated by reference numeral 1400 is partially transmitted through a beam splitter 1402 and subsequently focused, as by a lens 1404 onto a phase manipulator 1406, which is preferably located at the focal plane of lens 1404. The phase manipulator 1406 may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects.

[0365] A reflecting surface 1408 is arranged so as to reflect wavefront 1400 after it passes through the phase manipulator 1406. The reflected wavefront is imaged by lens 1404 onto a detector 1410, such as a CCD detector via beam splitter 1402. Preferably the beam splitter 1402 and the detector 1410 are arranged such that the detector 1410 lies in the focal plane of lens 1404. The output of detector 1410 is preferably supplied to data storage and processing circuitry 1412, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A.

[0366] It is appreciated that adding the reflecting surface 1408 to an imaging system, doubles the phase delay generated by phase manipulator 1406, enables imaging with a single lens 1404, and generally enables realization of a more compact system.

[0367] Reference is now made to FIG. 15, which is a simplified partially schematic, partially pictorial illustration of a system for surface mapping employing the functionality and structure of FIGS. 1A and 1B. As seen in FIG. 15, a beam of radiation, such as light or acoustic energy, is supplied from a radiation source 1500, optionally via a beam expander 1502, onto a beam splitter 1504, which reflects at least part of the radiation onto a surface 1506 to be inspected. The radiation reflected from the inspected surface 1506, is a surface mapping wavefront, which has an amplitude and a phase, and which contains information about the surface 1506. At least part of the radiation incident on surface 1506 is reflected from the surface 1506 and transmitted via the beam splitter 1504 and focused via a focusing lens 1508 onto a phase manipulator 1510, which is preferably located at the image plane of radiation source 1500.

[0368] The phase manipulator 1510 may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects. It is appreciated that phase manipulator 1510 can be configured such that a substantial part of the radiation focused thereonto is reflected therefrom. Alternatively the phase manipulator 1510 can be configured such that a substantial part of the radiation focused thereonto is transmitted therethrough.

[0369] A second lens 1512 is arranged so as to image surface 1506 onto a detector 1514, such as a CCD detector. Preferably the second lens 1512 is arranged such that the detector 1514 lies in its focal plane. The output of detector 1514, an example of which is a set of intensity maps designated by reference numeral 1515, is preferably supplied to data storage and processing circuitry 1516, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A, providing an output indicating at least one and possibly both of the phase and the amplitude of the surface mapping wavefront. This output is preferably further processed to obtain information about the surface 1506, such as geometrical variations and reflectivity of the surface.

[0370] In accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1500 has a narrow wavelength band about a given central wavelength, causing the phase of the radiation reflected from surface 1506 to be proportional to geometrical variations in the surface 1506, the proportion being an inverse linear function of the central wavelength of the radiation.

[0371] In accordance with another preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1500 has at least two narrow wavelength bands, each centered about a different wavelength, designated λ₁, . . . , λ_(n). In such a case, the radiation reflected from the surface 1506 has at least two wavelength components, each centered around a wavelength λ₁, . . . , λ_(n) and at least two indications of the phase of the surface mapping wavefront are obtained. Each such indication corresponds to a different wavelength component of the reflected radiation. These at least two indications may be subsequently combined to enable enhanced mapping of the surface 1506, by avoiding ambiguity in the mapping, known as 2π ambiguity, when the value of the mapping at a given spatial location in the surface exceeds the value of the mapping at a different spatial location in the surface by the largest of the different wavelengths λ₁, . . . , λ_(n). A proper choice of the wavelengths λ₁, . . . , λ_(n), may lead to elimination of this ambiguity when the difference in values of the mapping at different locations is smaller than the multiplication product of all the wavelengths.

[0372] In accordance with still another preferred embodiment of the present invention, the phase manipulator 1510 applies a plurality of different spatial phase changes to the radiation wavefront reflected from surface 1506 and Fourier transformed by lens 1508. Application of the plurality of different spatial phase changes provides a plurality of differently phase changed transformed wavefronts which may be subsequently detected by detector 1514.

[0373] In accordance with yet another preferred embodiment of the present invention, at least three different spatial phase changes are applied by phase manipulator 1510, resulting in at least three different intensity maps 1515. The at least three intensity maps are employed by the data storage and processing circuitry 1516 to obtain an output indicating at least the phase of the surface mapping wavefront. In such a case, the data storage and processing circuitry 1516, carries out functionality “C” described hereinabove with reference to FIG. 1A, preferably in a manner described hereinabove with reference to FIG. 13, where the wavefront being analyzed (FIG. 13) is the surface mapping wavefront.

[0374] Additionally, in accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1500 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the surface mapping wavefront and subsequently in the transformed wavefront impinging on phase manipulator 1510. In this case the phase manipulator may be an object, at least one of whose thickness, refractive index and surface geometry varies spatially. This spatial variance of the phase manipulator generates a different spatial phase change for each of the wavelength components, thereby providing a plurality of differently phase changed transformed wavefronts to be subsequently detected by detector 1514.

[0375] Reference is now made to FIG. 16, which is a simplified partially schematic, partially pictorial illustration of a system for object inspection employing the functionality and structure of FIGS. 1A and 1B. As seen in FIG. 16, a beam of radiation, such as light or acoustic energy, is supplied from a radiation source 1600, optionally via a beam expander, onto at least partially transparent object to be inspected 1602. The radiation transmitted through the inspected object 1602, is an object inspection wavefront, which has an amplitude and a phase, and which contains information about the object 1602. At least part of the radiation transmitted through object 1602 is focused via a focusing lens 1604 onto a phase manipulator 1606, which is preferably located at the image plane of radiation source 1600.

[0376] The phase manipulator 1606 may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects. It is appreciated that phase manipulator 1606 can be configured such that a substantial part of the radiation focused thereonto is reflected therefrom. Alternatively the phase manipulator 1606 can be configured such that a substantial part of the radiation focused thereonto is transmitted therethrough.

[0377] A second lens 1608 is arranged so as to image object 1602 onto a detector 1610, such as a CCD detector. Preferably, the second lens 1608 is arranged such that the detector 1610 lies in its focal plane. The output of detector 1610, an example of which is a set of intensity maps designated by reference numeral 1612, is preferably supplied to data storage and processing circuitry 1614, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A, providing an output indicating at least one and possibly both of the phase and the amplitude of the object inspection wavefront. This output is preferably further processed to obtain information about the object 1602, such as a mapping of the object's thickness, refractive index or transmission.

[0378] In accordance with one preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1600 has a narrow wavelength band about a given central wavelength, and the object 1602 is substantially uniform in material and other optical properties, causing the phase of the radiation transmitted through object 1602 to be proportional to thickness of the object 1602.

[0379] In accordance with one more preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1600 has a narrow wavelength band about a given central wavelength, and the object 1602 is substantially uniform in thickness, causing the phase of the radiation transmitted through object 1602 to be proportional to optical properties, such as refraction index or density, of the object 1602. It is appreciated that object 1602 may be any optical conduction element, such as an optical fiber.

[0380] In accordance with another preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1600 has at least two narrow wavelength bands, each centered about a different wavelength, designated λ₁, . . . , λ_(n). In such a case, the radiation transmitted through object 1602 has at least two wavelength components, each centered around a wavelength λ₁, . . . , λ_(n) and at least two indications of the phase of the object inspection wavefront are obtained. Each such indication corresponds to a different wavelength component of the transmitted radiation. These at least two indications may be subsequently combined to enable enhanced mapping of the properties, such as thickness, of object 1602, by avoiding ambiguity in the mapping, known as 2π ambiguity, when the value of the mapping at a given spatial location in the object exceeds the value of the mapping at a different spatial location in the object by the largest of the different wavelengths λ₁, . . . , λ_(n). A proper choice of the wavelengths λ₁, . . . , λ_(n), may lead to elimination of this ambiguity when the difference in values of the mapping at different locations is smaller than the multiplication product of all the wavelengths.

[0381] In accordance with still another preferred embodiment of the present invention, the phase manipulator 1606 applies a plurality of different spatial phase changes to the radiation wavefront transmitted through object 1602 and Fourier transformed by lens 1604. Application of the plurality of different spatial phase changes produces a plurality of differently phase changed transformed wavefronts which may be subsequently detected by detector 1610.

[0382] In accordance with yet another preferred embodiment of the present invention, at least three different spatial phase changes are applied by phase manipulator 1606, resulting in at least three different intensity maps 1612. The at least three intensity maps 1612 are employed by the data storage and processing circuitry 1614 to obtain an output indicating at least the phase of the object inspection wavefront. In such a case, the data storage and processing circuitry 1614, carries out functionality “C” described hereinabove with reference to FIG. 1A, preferably in a manner described hereinabove with reference to FIG. 13, where the wavefront being analyzed (FIG. 13) is the object inspection wavefront.

[0383] Additionally, in accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1600 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the object inspection wavefront and subsequently in the transformed wavefront impinging on phase manipulator 1606. In this case the phase manipulator 1606 may be an object, at least one of whose thickness, refractive index and surface geometry varies spatially. This spatial variance of the phase manipulator generates a different spatial phase change for each of the wavelength components, thereby providing a plurality of differently phase changed transformed wavefronts to be subsequently detected by detector 1610.

[0384] Reference is now made to FIG. 17, which is a simplified partially schematic, partially pictorial illustration of a system for spectral analysis employing the functionality and structure of FIGS. 1A and 1B. As seen in FIG. 17, a beam of radiation, such as light or acoustic energy, is supplied from a radiation source to be tested 1700, optionally via a beam expander, onto a known element 1702, such as an Etalon or a plurality of Etalons. Element 1702 is intended to generate an input wavefront, having at least varying phase or intensity. The radiation transmitted through the element 1702, is a spectral analysis wavefront, which has an amplitude and a phase, and which contains information about the spectrum of the radiation source 1700. At least part of the radiation transmitted through element 1702 is focused via a focusing lens 1704 onto a phase manipulator 1706, which is preferably located at the image plane of radiation source 1700.

[0385] The phase manipulator 1706 may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects. It is appreciated that phase manipulator 1706 can be configured such that a substantial part of the radiation focused thereonto is reflected therefrom. Alternatively the phase manipulator 1706 can be configured such that a substantial part of the radiation focused thereonto is transmitted therethrough.

[0386] A second lens 1708 is arranged so as to image element 1702 onto a detector 1710, such as a CCD detector. Preferably, the second lens 1708 is arranged such that the detector 1710 lies in its focal plane. The output of detector 1710, an example of which is a set of intensity maps designated by reference numeral 1712, is preferably supplied to data storage and processing circuitry 1714, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A, providing an output indicating at least one and possibly both of the phase and the amplitude of the spectral analysis wavefront. This output is preferably further processed to obtain information about the radiation source 1700, such as the spectrum of the radiation supplied from radiation source 1700.

[0387] In accordance with a preferred embodiment of the present invention, the spectral analysis wavefront is obtained by reflecting the radiation supplied from radiation source 1700 from element 1702.

[0388] In accordance with another preferred embodiment of the present invention, the spectral analysis wavefront is obtained by transmitting the radiation supplied from radiation source 1700 through element 1702.

[0389] In accordance with one more preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1700 has a narrow wavelength band about a central wavelength, causing the phase of the radiation impinged on the object 1702 to be inversely proportional to the central wavelength supplied from radiation source 1700 and related to at least one of a surface characteristic and thickness of element 1702.

[0390] In accordance with another preferred embodiment of the present invention, the plurality of intensity maps 1712 are employed by the data storage and processing circuitry 1714, to obtain an output indicating at least one and possibly both of the phase and amplitude of the spectral analysis wavefront by expressing the plurality of intensity maps as at least one mathematical function of phase and amplitude of the spectral analysis wavefront and of plurality of different phase changes applied by phase manipulator 1706, wherein at least one and possibly both of the phase and amplitude is unknown and a function generating the different phase changes is known. This at least one mathematical function is subsequently employed to obtain an output indicating at least the phase of the spectral analysis wavefront.

[0391] In accordance with still another preferred embodiment of the present invention, the phase manipulator 1706 applies a plurality of different spatial phase changes to the radiation wavefront transmitted through element 1702 and Fourier transformed by lens 1704. Application of the plurality of different spatial phase changes produces a plurality of differently phase changed transformed wavefronts which may be subsequently detected by detector 1710.

[0392] In accordance with yet another preferred embodiment of the present invention, at least three different spatial phase changes are applied by phase manipulator 1706, resulting in at least three different intensity maps 1712. The at least three intensity maps are employed by the data storage and processing circuitry 1714 to obtain an output indicating at least the phase of the spectral analysis wavefront. In such a case, the data storage and processing circuitry 1714, carries out functionality “C” described hereinabove with reference to FIG. 1A, preferably in a manner described hereinabove with reference to FIG. 13, where the wavefront being analyzed (FIG. 13) is the spectral analysis wavefront.

[0393] Additionally, in accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1700 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the spectral analysis wavefront and subsequently in the transformed wavefront impinging on phase manipulator 1706. In this case the phase manipulator may be an object, at least one of whose thickness, refractive index and surface geometry varies spatially. This spatial variance of the phase manipulator generates a different spatial phase change for each of the wavelength components, thereby providing a plurality of differently phase changed transformed wavefronts to be subsequently detected by detector 1710.

[0394] In accordance with an embodiment of the present invention, the phase manipulator 1706 comprises a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially. The spatial variance of the thickness or of the refractive index of the plurality of objects may be selected in a way such that the phase changes applied by phase manipulator 1706 differ to a selected predetermined extent for at least some of the wavelength components supplied by radiation source 1700.

[0395] A specific selection of the objects is such that the phase change applied to an expected wavelength of radiation source differs substantially from the phase change applied to an actual wavelength of the radiation source. Alternatively, the spatial variance of the thickness or refractive index of the plurality of objects may be selected in a way such that the phase changes applied by phase manipulator 1706 are identical for at least some of the plurality of different wavelength components wavelength components supplied by radiation source 1700.

[0396] In accordance with another embodiment of the present invention, the known element 1702 comprises a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially. The spatial variance of the thickness or of the refractive index of the plurality of objects may be selected in a way such that the wavelength components of the input wavefront, generated by passing the wavelength components of the radiation supplied by radiation source 1700 through the element 1702, differ to a selected predetermined extent for at least some of the wavelength components supplied by radiation source 1700.

[0397] A specific selection of the objects is such that the wavelength component of the input wavefront generated by an expected wavelength of radiation source differs substantially from the wavelength component of the input wavefront generated by an actual wavelength of the radiation source. Alternatively, the spatial variance of the thickness or refractive index of the plurality of objects may be selected in a way such that the wavelength components of the input wavefront, generated by passing the wavelength components of the radiation supplied by radiation source 1700 through the element 1702, are identical for at least some of the wavelength components supplied by radiation source 1700.

[0398] Reference is now made to FIG. 18, which is a simplified partially schematic, partially pictorial illustration of a system for phase-change analysis employing the functionality and structure of FIGS. 1A and 1B. As seen in FIG. 18, a known wavefront 1800, which is a phase change analysis wavefront, having an amplitude and a phase, is focused via a focusing lens 1802, preferably performing a Fourier transform to wavefront 1800, onto a phase manipulator 1804, which is preferably located at the focal plane of lens 1802. The phase manipulator applies a plurality of different phase changes to the transformed phase change analysis wavefront.

[0399] The phase manipulator 1804 may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects. It is appreciated that phase manipulator 1804 can be configured such that a substantial part of the radiation focused thereonto is reflected therefrom. Alternatively the phase manipulator 1804 can be configured such that a substantial part of the radiation focused thereonto is transmitted therethrough.

[0400] A second lens 1806 is arranged so as to image wavefront 1800 onto a detector 1808, such as a CCD detector. Preferably, the second lens 1806 is arranged such that the detector 1808 lies in its focal plane. The output of detector 1808, an example of which is a set of intensity maps designated by reference numeral 1810, is preferably supplied to data storage and processing circuitry 1812, which employs the plurality of intensity maps to obtain an output indication of differences between the plurality of different phase changes applied by the phase manipulator 1804.

[0401] In accordance with one preferred embodiment of the present invention, lateral shifts appear in the plurality of different phase changes. These may be produced, for example, by vibrations of the phase manipulator or by impurities in the phase manipulator. Consequently, corresponding changes appear in the plurality of intensity maps 1810, and result in obtaining an indication of these lateral shifts.

[0402] In accordance with another preferred embodiment of the present invention, the plurality of intensity maps 1810 are employed by the data storage and processing circuitry 1812 to obtain an output indicating the differences between the plurality of different phase changes applied by the phase manipulator 1804, by expressing the plurality of intensity maps as at least one mathematical function of phase and amplitude of the phase change analysis wavefront and of the plurality of different phase changes applied by phase manipulator 1804, where at least the phase and amplitude of the wavefront 1800 are known and the plurality of different phase changes are unknown. This at least one mathematical function is subsequently employed to obtain an output indicating at least the differences between the plurality of different phase changes.

[0403] In accordance with still another preferred embodiment of the present invention, the phase manipulator 1804 applies a plurality of different spatial phase changes to the wavefront 1800 Fourier transformed by lens 1802. Application of the plurality of different spatial phase changes provides a plurality of differently phase changed transformed wavefronts which may be subsequently detected by detector 1808.

[0404] In accordance with yet another preferred embodiment of the present invention, at least three different spatial phase changes are applied by phase manipulator 1804, resulting in at least three different intensity maps 1810. The at least three intensity maps are employed by the data storage and processing circuitry 1812 to obtain an output indicating at least the differences between the plurality of different phase changes. In such a case, the data storage and processing circuitry 1814, carries out functionality “C” described hereinabove with reference to FIG. 1A, preferably in a manner similar to the manner described hereinabove with reference to FIG. 13, where the wavefront being analyzed (FIG. 13) is the known phase change analysis wavefront, and the spatial phase changes applied by phase manipulator 1804 are unknown.

[0405] Additionally, in accordance with a preferred embodiment of the present invention, the wavefront 1800 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the transformed wavefront impinging on phase manipulator 1804. In this case the phase manipulator may be an object, at least one of whose thickness, refractive index and surface geometry varies spatially. This spatial variance of the phase manipulator generates a different spatial phase change for each of the wavelength components, thereby providing a plurality of differently phase changed transformed wavefronts to be subsequently detected by detector 1808.

[0406] Additionally, in accordance with another embodiment of the present invention, phase manipulator 1804 applies one phase change to the radiation focused onto each spatial location thereon, resulting in one intensity map 1810, as an output of detector 1808. In such a case, the data storage and processing circuitry 1812 employs the intensity map and the known wavefront 1800 to obtain at least an output indicating the phase change applied by phase manipulator 1804.

[0407] In accordance with the foregoing methodology, the phase change applied by the phase manipulator may be a phase delay, having a value selected from one of a plurality of pre-determined values, including a possible value of zero phase delay and the output indication of the phase change obtained by data storage and processing circuitry 1812 is the value of the phase delay. In such a case, the phase manipulator may be media which stores information by different values of the phase delays at each of a multiplicity of different locations thereon, where the value of the phase delay constitutes the stored information. The stored information, encoded in the different values of the phase delays, is retrieved by data storage and processing circuitry 1812. It is appreciated that in such a case, wavefront 1800 may also comprise a plurality of different wavelength components, resulting in a plurality of intensity maps and consequently in an increase of the information encoded on the phase manipulator at each of the multiplicity of different locations.

[0408] Reference is now made to FIG. 19, which is a simplified partially schematic, partially pictorial illustration of a system for stored data retrieval employing the functionality and structure of FIGS. 1A and 1B. As seen in FIG. 19, optical storage media 1900, such as a DVD or compact disk, has information encoded thereon by selecting the height of the media at each of a multiplicity of different locations thereon, as shown in enlargement and designated by reference numeral 1902. At each location on the media, the height of the media can be one of several given heights or levels. The specific level of the media at that location determines the information stored at that location.

[0409] A beam of radiation, such as light or acoustic energy, is supplied from a radiation source 1904, such as a laser or a LED, optionally via a beam expander, onto a beam splitter 1906, which reflects at least part of the radiation onto the surface of the media 1900. The radiation reflected from an area 1908 on the media, onto which the radiation impinges, is a stored data retrieval wavefront, which has an amplitude and a phase, and which contains information stored in area 1908. At least part of the radiation incident on area 1908 is reflected from the area 1908 and transmitted via the beam splitter 1906 onto an imaging system 1910, which may include a phase manipulator or other device which generates a varying phase function.

[0410] Imaging system 1910 preferably carries out functionalities “A” and “B” described hereinabove with reference to FIG. 1A, obtaining a plurality of differently phase changed transformed wavefronts corresponding to the stored data retrieval wavefront and obtaining a plurality of intensity maps of the plurality of phase changed transformed wavefronts.

[0411] Preferably, imaging system 1910 comprises a first lens 1508 (FIG. 15), a phase manipulator 1510 (FIG. 15), a second lens 1512 (FIG. 15) and a detector 1514 (FIG. 15). The outputs of imaging system 1910 are supplied to data storage and processing circuitry 1912, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A, providing an output indicating at least one and possibly both of the phase and amplitude of the stored data retrieval wavefront. This output is preferably further processed to read out the information encoded in area 1908 of media 1900 and displayed on display unit 1914.

[0412] In accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1904 has a narrow wavelength band about a given central wavelength, causing the phase of the radiation reflected from media 1900 to be proportional to geometrical variations in the media 1900, which contain the encoded information, the proportion being an inverse linear function of the central wavelength of the radiation.

[0413] In accordance with another preferred embodiment of the present invention, the beam of radiation supplied from radiation source 1904 has at least two narrow wavelength bands, each centered about a different wavelength, designated λ₁, . . . , λ_(n). In such a case, the radiation reflected from area 1908 in media 1900 has at least two wavelength components, each centered around a wavelength λ₁, . . . , λ_(n).

[0414] At least two indications of the phase of the stored data retrieval wavefront are obtained, each such indication corresponding to a different wavelength component of the reflected radiation. These at least two indications may be subsequently combined to enhance mapping of the surface of area 1908 of media 1900 and therefore enhance retrieval of the information, by avoiding an ambiguity in the mapping, known as 2π ambiguity, when the value of the height of the media at a given location exceeds the largest of the different wavelengths λ₁, . . . , λ_(n).

[0415] In such a case, the range of possible heights at each location in area 1908 can exceed the value of the largest of the different wavelengths, without ambiguity in the reading of the heights. This extended dynamic range enables storing more information on media 1900 than would otherwise be possible. A proper choice of the wavelengths λ₁, . . . , λ_(n), may lead to elimination of this ambiguity when the difference of heights of the media in area 1908 at different locations is smaller than the multiplication product of all the wavelengths.

[0416] In accordance with still another preferred embodiment of the present invention, a phase manipulator incorporated in imaging system 1910 applies a plurality of different spatial phase changes to the radiation wavefront reflected from media 1900 and Fourier transformed by a lens, also incorporated in imaging system 1910. Application of the plurality of different spatial phase changes provides a plurality of differently phase changed transformed wavefronts which may be subsequently detected by a detector incorporated in imaging system 1910.

[0417] In accordance with yet another preferred embodiment of the present invention, at least three different spatial phase changes are applied by a phase manipulator incorporated in imaging system 1910, resulting in an output from imaging system 1910 of at least three different intensity maps. The at least three intensity maps are employed by the data storage and processing circuitry 1912 to obtain an output indicating at least the phase of the stored data retrieval wavefront. In such a case, the data storage and processing circuitry 1912, carries out functionality “C” described hereinabove with reference to FIG. 1A, preferably in a manner described hereinabove with reference to FIG. 13, where the wavefront being analyzed (FIG. 13) is the stored data retrieval wavefront.

[0418] Additionally, in accordance with an embodiment of the present invention, the beam of radiation supplied from radiation source 1904 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the stored data retrieval wavefront and subsequently in the transformed wavefront impinging on a phase manipulator incorporated into imaging system 1910. In this case the phase manipulator may be an object, at least one of whose thickness, refractive index and surface geometry varies spatially. This spatial variance of the phase manipulator generates a different spatial phase change for each of the wavelength components, thereby providing a plurality of differently phase changed transformed wavefronts to be subsequently detected by a detector incorporated in imaging system 1910.

[0419] In accordance with another embodiment of the present invention, information is encoded on media 1900 by selecting the height of the media at each given location to be such that the intensity value of the intensity map resulting from light reflected from the location and passing through imaging system 1910 lies within a predetermined range of values. This range corresponds to an element of the information stored at the location. By employing the plurality of intensity maps, multiple intensity values are realized for each location, each intensity value lying within a specific range of values. The resulting plurality of ranges of intensity values provide multiple elements of information for each location on the media 1900.

[0420] It is appreciated that in such a case, retrieving the information stored at area 1908 on the media from the outputs of imaging system 1910 may be performed by data storage and processing circuitry 1912 in a straight-forward manner, as by mapping each of the resulting intensity values at each location to their corresponding ranges, and subsequently to the information stored at the location.

[0421] Preferably, the foregoing methodology also includes use of a phase manipulator incorporated in imaging system 1910, that applies an at least time-varying phase change function to the transformed data retrieval wavefront impinging thereon. This time-varying phase change function provides the plurality of intensity maps.

[0422] Alternatively or additionally, the beam of radiation supplied from radiation source 1904 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the stored data retrieval wavefront. The plurality of differently phase changed transformed wavefronts are obtained in imaging system 1910 by applying at least one phase change to the plurality of different wavelength components of the stored data retrieval wavefront. The phase changed transformed stored data retrieval wavefront can be detected by a single detector or alternatively separated, as by a dispersion element, into its separate plurality of different wavelength components, each component being detected by a different detector.

[0423] In accordance with yet another embodiment of the present invention, media 1900 may have different reflectivity coefficients for the radiation supplied from light source 1904 at each of a multiplicity of different locations on the media. At each location on the media, a different percentage of the radiation may be reflected. The percentage may have one of several given values, where the specific value may at least partially determine the information stored at that location.

[0424] In such a case, the information encoded on media 1900 may be encoded by selecting the height of the media at each of a multiplicity of different locations on the media and by selecting the reflectivity of the media at each of a multiplicity of different locations on the media. In such a case, more information can be stored at each location on the media, than could otherwise be stored. Furthermore, in such a case, employing an indication of the amplitude and phase of the stored data retrieval wavefront to obtain the encoded information includes employing the indication of the phase to obtain the information encoded by selecting the height of the media and employing the indication of the amplitude to obtain said information encoded by selecting the reflectivity.

[0425] In accordance with still another embodiment of the present invention, the information is encoded onto media 1900 at several layers in the media. The information is encoded on the media by selecting the height of the media at each of multiplicity of different locations on each layer of the media. Each of these layers, placed one on top of the other in media 1900, is partially reflecting and partially transmitting. The beam of radiation from source 1904 impinging onto media 1900 is partially reflected from the top, first layer of the media, and partially transmitted to the layers lying therebelow. The energy transmitted by the second layer is partially reflected and partially transmitted to the layers lying therebelow, and so on, until the radiation transmitted through all the layers is partially reflected from the undermost layer. In such a case, radiation source 1904 preferably includes a focusing system that focuses the radiation onto each one of the layers of media 1900 in order to retrieve the information stored on that layer. Alternatively or additionally, imaging system 1910 may include confocal microscopy elements operative to differentiate between the different layers.

[0426] It is appreciated that area 1908 of media 1900 may be a relatively small area, comprising a single location on which information is encoded and possibly also neighboring locations. In such a case, the detector incorporated in imaging system 1910 may define only a single or several detection pixels. Additionally, the output indicating at least one and possibly both of the phase and amplitude of the stored data retrieval wavefront provided by circuitry 1912, includes at least one and possibly both of the height of the media and the reflectivity of the media at the location or locations covered by area 1908.

[0427] In accordance with yet another embodiment of the present invention, the stored data retrieval wavefront comprises at least one one-dimensional component, corresponding to at least one one-dimensional part of area 1908 on media 1900. In such a case, the imaging system 1910 is preferably similar to the imaging system described hereinabove with reference to FIG. 10B. It preferably includes a first lens, such as cylindrical lens 1086 (FIG. 10B).

[0428] The first lens preferably produces a one-dimensional Fourier transforrn, performed along an axis extending perpendicularly to the direction of propagation of the data retrieval wavefront, thereby providing at least one one-dimensional component of the transformed wavefront in a dimension perpendicular to direction of propagation. The first lens, such as lens 1086, focuses the stored data retrieval wavefront onto a phase manipulator, such as a single axis displaceable phase manipulator 1087 (FIG. 10B), which is preferably located at the focal plane of lens 1086. The phase manipulator 1087 applies a plurality of different phase changes to each of the at least one one-dimensional components of the transformed wavefront, thereby obtaining at least one one-dimensional component of the plurality of phase changed transformed wavefronts.

[0429] Additionally the imaging system may include a second lens, such as cylindrical lens 1088 (FIG. 10B), arranged so as to image the at least one one-dimensional component of the stored data retrieval wavefront onto a detector 1089, such as a CCD detector. Additionally the plurality of intensity maps are employed by circuitry 1912 to obtain an output indicating al least one and possibly both of the amplitude and phase of the at least one one-dimensional component of the data retrieval wavefront.

[0430] Additionally, in accordance with the foregoing methodology, and as described hereinabove with reference to FIG. 10B, the phase manipulator 1087 preferably comprises a multiple local phase delay element, such as a spatially nonuniform transparent object, typically including several different phase delay regions, each arranged to apply a phase delay to one of the at least one one-dimensional component at a given position of the object along a phase manipulator axis, extending perpendicularly to the direction of propagation of the wavefront and perpendicular to the axis of the transform produced by lens 1086.

[0431] In such a case, there is provided relative movement between the imaging system 1910 and the media 1900 along the phase manipulator axis. This relative movement sequentially matches different phase delay regions with different wavefront components, corresponding to different parts of area 1908 on media 1900, such that preferably each wavefront component passes through each phase delay region of the phase manipulator.

[0432] It is appreciated that the relative movement between the imaging system 1910 and the at least one one-dimensional wavefront component can be obtained by the rotation of media 1900 about its axis, while the imaging system is not moving.

[0433] It is a particular feature of this embodiment, that each of the at least one one-dimensional component of the wavefront is separately processed. Thus, each of the at least one one-dimensional wavefront component, corresponding to a one-dimensional part of area 1908, is focused by a separate portion of the first cylindrical lens of imaging system 1910, is imaged by a corresponding separate portion of the second cylindrical lens and passes through a distinct region of the phase manipulator. The images of each of the one-dimensional parts of area 1908 at the detector incorporated in imaging system 1910 are thus separate and distinct images. It is appreciated that these images may appear on separate detectors or on a monolithic detector.

[0434] In accordance with an embodiment of the present invention, the transform applied to the stored data retrieval wavefront includes an additional Fourier transform. This additional Fourier transform may be performed by the first cylindrical lens of imaging system 1910 or by an additional lens and is operative to minimize cross-talk between different one-dimensional components of the wavefront. In such a case, preferably an additional transform, such as that provided by an additional lens adjacent to the second cylindrical lens, is applied to the phase changed transformed wavefront. In such a case, preferably a further transform is applied to the phase changed transformed wavefront. This further transform may be performed by the second cylindrical lens of imaging system 1910 or by an additional lens.

[0435] Reference is now made to FIG. 20, which is a simplified partially schematic, partially pictorial illustration of a system for 3-dimensional imaging employing the functionality and structure of FIGS. 1A and 1B. As seen in FIG. 20, a beam of radiation, such as light or acoustic energy, is supplied from a radiation source 2000, optionally via a beam expander, onto a beam splitter 2004, which reflects at least part of the radiation onto a 3-dimensional object 2006 to be imaged. The radiation reflected from the object 2006, is a 3-dimensional imaging wavefront, which has an amplitude and a phase, and which contains information about the object 2006. At least part of the radiation incident on the surface of object 2006 is reflected from the object 2006 and transmitted via the beam splitter 2004 and focused via a focusing lens 2008 onto a phase manipulator 2010, which is preferably located at the image plane of radiation source 2000.

[0436] The phase manipulator 2010 may be, for example, a spatial light modulator or a series of different transparent, spatially non-uniform objects. It is appreciated that phase manipulator 2010 can be configured such that a substantial part of the radiation focused thereonto is reflected therefrom. Alternatively the phase manipulator 2010 can be configured such that a substantial part of the radiation focused thereonto is transmitted therethrough.

[0437] A second lens 2012 is arranged so as to image object 2006 onto a detector 2014, such as a CCD detector. Preferably the second lens 2012 is arranged such that the detector 2014 lies in its focal plane. The output of detector 2014, an example of which is a set of intensity maps designated by reference numeral 2015, is preferably supplied to data storage and processing circuitry 2016, which preferably carries out functionality “C” described hereinabove with reference to FIG. 1A, providing an output indicating at least one and possibly both of the phase and amplitude of the 3-dimensional imaging wavefront. This output is preferably further processed to obtain information about the object 2006, such as the 3-dimensional shape of the object.

[0438] In accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 2000 has a narrow wavelength band about a given central wavelength, causing the phase of the radiation reflected from object 2006 to be proportional to geometrical variations in the surface 2006, the proportion being an inverse linear function of the central wavelength of the radiation.

[0439] In accordance with another preferred embodiment of the present invention, the beam of radiation supplied from radiation source 2000 has at least two narrow wavelength bands, each centered about a different wavelength, designated λ₁, . . . , λ_(n). In such a case, the radiation reflected from the object 2006 has at least two wavelength components, each centered around a wavelength λ₁, . . . , λ_(n) and at least two indications of the phase of the 3-dimensional imaging wavefront are obtained. Each such indication corresponds to a different wavelength component of the reflected radiation. These at least two indications may be subsequently combined to enable enhanced imaging of the object 2006, by avoiding 2π ambiguity in the 3-dimensional imaging.

[0440] In accordance with still another preferred embodiment of the present invention, the phase manipulator 2010 applies a plurality of different spatial phase changes to the radiation wavefront reflected from surface 2006 and Fourier transformed by lens 2008. Application of the plurality of different spatial phase changes provides a plurality of differently phase changed transformed wavefronts which may be subsequently detected by detector 2014.

[0441] In accordance with yet another preferred embodiment of the present invention, at least three different spatial phase changes are applied by phase manipulator 2010, resulting in at least three different intensity maps 2015. The at least three intensity maps are employed by the data storage and processing circuitry 2016 to obtain an output indicating at least the phase of the 3-dimensional imaging wavefront. In such a case, the data storage and processing circuitry 2016, carries out functionality “C” described hereinabove with reference to FIG. 1A, preferably in a manner described hereinabove with reference to FIG. 13, where the wavefront being analyzed (FIG. 13) is the 3-dimensional imaging wavefront.

[0442] Additionally, in accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 2000 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the 3-dimensional imaging wavefront and subsequently in the transformed wavefront impinging on phase manipulator 2010. In this case the phase manipulator 2010 may be an object, at least one of whose thickness, refractive index and surface geometry varies spatially. This spatial variance of the phase manipulator generates a different spatial phase change for each of the wavelength components, thereby providing a plurality of differently phase changed transformed wavefronts to be subsequently detected by detector 2014.

[0443] Reference is now made to FIG. 21A, which is a simplified partially schematic, partially pictorial illustration of wavefront analysis functionality operative in accordance with another preferred embodiment of the present invention. The functionality of FIG. 21 A can be summarized as including the following sub-functionalities:

[0444] A. obtaining a plurality of differently amplitude changed transformed wavefronts corresponding to a wavefront being analyzed, which has an amplitude and a phase;

[0445] B. obtaining a plurality of intensity maps of the plurality of amplitude changed transformed wavefronts; and

[0446] C. employing the plurality of intensity maps to obtain an output indicating at least one and possibly both of the phase and the amplitude of the wavefront being analyzed.

[0447] As seen in FIG. 21A, the first sub-functionality, designated “A” may be realized by the following functionalities:

[0448] A wavefront, which may be represented by a plurality of point sources of light, is generally designated by reference numeral 2100. Wavefront 2100 has a phase characteristic which is typically spatially non-uniform, shown as a solid line and indicated generally by reference numeral 2102. Wavefront 2100 also has an amplitude characteristic which is typically spatially non-uniform, shown as a dashed line and indicated generally by reference numeral 2103. Such a wavefront may be obtained in a conventional manner by receiving light from any suitable object, such as by reading an optical disk, for example, a DVD or compact disk 2104.

[0449] A principal purpose of the present invention is to measure the phase characteristic, such as that indicated by reference numeral 2102, which is not readily measured. Another purpose of the present invention is to measure the amplitude characteristic, such as that indicated by reference numeral 2103 in an enhanced manner. A further purpose of the present invention is to measure both the phase characteristic 2102 and the amplitude characteristic 2103. While there exist various techniques for carrying out such measurements, the present invention provides a methodology which is believed to be superior to those presently known, inter alia due to its relative insensitivity to noise.

[0450] A transform, indicated here symbolically by reference numeral 2106, is applied to the wavefront being analyzed 2100, thereby to obtain a transformed wavefront. A preferred transform is a Fourier transform. The resulting transformed wavefront is symbolically indicated by reference numeral 2108.

[0451] A plurality of different amplitude changes, preferably spatial amplitude changes, represented by optical attenuation components 2110, 2112 and 2114 are applied to the transformed wavefront 2108, thereby to obtain a plurality of differently amplitude changed transformed wavefronts, represented by reference numerals 2120, 2122 and 2124 respectively. It is appreciated that the illustrated difference between the individual ones of the plurality of differently amplitude changed transformed wavefronts is that portions of the transformed wavefront are attenuated differently relative to the remainder thereof.

[0452] As seen in FIG. 21A, the second sub-functionality, designated “B”, may be realized by applying a transform, preferably a Fourier transform, to the plurality of differently amplitude changed transformed wavefronts. Alternatively, the sub-functionality B may be realized without the use of a Fourier transform, such as by propagation of the differently amplitude changed transformed wavefronts over an extended space. Finally, functionality B requires detection of the intensity characteristics of plurality of differently amplitude changed transformed wavefronts. The outputs of such detection are the intensity maps, examples of which are designated by reference numerals 2130, 2132 and 2134.

[0453] As seen in FIG. 21A, the third sub-functionality, designated “C” may be realized by the following functionalities:

[0454] expressing, such as by employing a computer 2136, the plurality of intensity maps, such as maps 2130, 2132 and 2134, as at least one mathematical function of phase and amplitude of the wavefront being analyzed and of the plurality of different amplitude changes, wherein at least one and possibly both of the phase and the amplitude are unknown and the plurality of different amplitude changes, typically represented by optical attenuation components 2110, 2112 and 2114 applied to the transformed wavefront 2108, are known; and

[0455] employing, such as by means of the computer 2136, the at least one mathematical function to obtain an indication of at least one and possibly both of the phase and the amplitude of the wavefront being analyzed, here represented by the phase function designated by reference numeral 2138 and the amplitude function designated by reference numeral 2139, which, as can be seen, respectively represent the phase characteristics 2102 and the amplitude characteristics 2103 of the wavefront 2100. In this example, wavefront 2100 may represent the information contained in the compact disk or DVD 2104.

[0456] In accordance with an embodiment of the present invention, the plurality of intensity maps comprises at least four intensity maps. In such a case, employing the plurality of intensity maps to obtain an output indicating at least the phase of the wavefront being analyzed includes employing a plurality of combinations, each of at least three of the plurality of intensity maps, to provide a plurality of indications at least of the phase of the wavefront being analyzed.

[0457] Preferably, the methodology also includes employing the plurality of indications of at least the phase of the wavefront being analyzed to provide an enhanced indication at least of the phase of the wavefront being analyzed.

[0458] Also in accordance with an embodiment of the present invention, the plurality of intensity maps comprises at least four intensity maps. In such a case, employing the plurality of intensity maps to obtain an output indicating at least the amplitude of the wavefront being analyzed includes employing a plurality of combinations, each of at least three of the plurality of intensity maps, to provide a plurality of indications at least of the amplitude of the wavefront being analyzed.

[0459] Preferably, the methodology also includes employing the plurality of indications of at least the amplitude of the wavefront being analyzed to provide an enhanced indication at least of the amplitude of the wavefront being analyzed.

[0460] It is appreciated that in this manner, enhanced indications of both phase and amplitude of the wavefront may be obtained.

[0461] In accordance with a preferred embodiment of the present invention, at least some of the plurality of indications of the amplitude and phase are at least second order indications of the amplitude and phase of the wavefront being analyzed.

[0462] In accordance with one preferred embodiment of the present invention, the plurality of intensity maps are employed to provide an analytical output indicating the amplitude and phase.

[0463] Preferably, the amplitude changed transformed wavefronts are obtained by interference of the wavefront being analyzed along a common optical path.

[0464] In accordance with another preferred embodiment of the present invention, the plurality of intensity maps are employed to obtain an output indicating the phase of the wavefront being analyzed, which is substantially free from halo and shading off distortions, which are characteristic of many of the existing ‘phase-contrast’ methods.

[0465] In accordance with still another embodiment of the present invention, the plurality of intensity maps may be employed to obtain an output indicating the phase of the wavefront being analyzed by combining the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, obtaining at least an output indicative of the phase of the wavefront being analyzed from each of the second plurality of combined intensity maps and combining the outputs to provide an enhanced indication of the phase of the wavefront being analyzed.

[0466] In accordance with yet another embodiment of the present invention, the plurality of intensity maps may be employed to obtain an output indicating amplitude of the wavefront being analyzed by combining the plurality of intensity maps into a second plurality of combined intensity maps, the second plurality being less than the first plurality, obtaining at least an output indicative of the amplitude of the wavefront being analyzed from each of the second plurality of combined intensity maps and combining the outputs to provide an enhanced indication of the amplitude of the wavefront being analyzed.

[0467] Additionally in accordance with a preferred embodiment of the present invention, the foregoing methodology may be employed for obtaining a plurality of differently amplitude changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of amplitude changed transformed wavefronts and employing the plurality of intensity maps to obtain an output of an at least second order indication of phase of the wavefront being analyzed.

[0468] Additionally or alternatively in accordance with a preferred embodiment of the present invention, the foregoing methodology may be employed for obtaining a plurality of differently amplitude changed transformed wavefronts corresponding to a wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of amplitude changed transformed wavefronts and employing the plurality of intensity maps to obtain an output of an at least second order indication of amplitude of the wavefront being analyzed.

[0469] In accordance with yet another embodiment of the present invention, the obtaining of the plurality of differently amplitude changed transformed wavefronts comprises applying a transform to the wavefront being analyzed, thereby to obtain a transformed wavefront, and then applying a plurality of different phase and amplitude changes to the transformed wavefront, where each of these changes can be a phase change, an amplitude change or a combined phase and amplitude change, thereby to obtain a plurality of differently phase and amplitude changed transformed wavefronts.

[0470] In accordance with yet another embodiment of the present invention, a wavefront being analyzed comprises at least two wavelength components. In such a case, obtaining a plurality of intensity maps also includes dividing the amplitude changed transformed wavefronts according to the at least two wavelength components in order to obtain at least two wavelength components of the amplitude changed transformed wavefronts and in order to obtain at least two sets of intensity maps, each set corresponding to a different one of the at least two wavelength components of the amplitude changed transformed wavefronts.

[0471] Subsequently, the plurality of intensity maps are employed to provide an output indicating the amplitude and phase of the wavefront being analyzed by obtaining an output indicative of the phase of the wavefront being analyzed from each of the at least two sets of intensity maps and combining the outputs to provide an enhanced indication of phase of the wavefront being analyzed. In the enhanced indication, there is no 2π ambiguity once the value of the phase exceeds 2π, which conventionally results when detecting phase of a single wavelength wavefront.

[0472] It is appreciated that the wavefront being analyzed may be an acoustic radiation wavefront.

[0473] It is also appreciated that the wavefront being analyzed may be an electromagnetic radiation wavefront, of any suitable wavelength, such as visible light, infrared, ultra-violet and X-ray radiation.

[0474] It is further appreciated that wavefront 2100 may be represented by a relatively small number of point sources and defined over a relatively small spatial region. In such a case, the detection of the intensity characteristics of the plurality of differently amplitude changed transformed wavefronts may be performed by a detector comprising only a single detection pixel or several detection pixels. Additionally, the output indicating at least one and possibly both of the phase and amplitude of the wavefront being analyzed may be provided by computer 2136 in a straight-forward manner.

[0475] In accordance with an embodiment of the present invention, the plurality of different amplitude changes 2110, 2112 and 2114, preferably spatial amplitude changes, are effected by applying a time-varying spatial amplitude change to part of the transformed wavefront 2108.

[0476] In accordance with a preferred embodiment of the present invention, the plurality of different amplitude changes 2110, 2112 and 2114 are effected by applying a spatially uniform, time-varying spatial amplitude change to part of the transformed wavefront 2108.

[0477] In accordance with an embodiment of the present invention, each of the wavefront 2100 and the transformed wavefront 2108 comprises a plurality of different wavelength components. In such a case, the plurality of different spatial amplitude changes may be effected by applying an amplitude change to each of the plurality of different wavelength components of the transformed wavefront. It is appreciated that the amplitude changes may be spatially different or that the amplitude may be differently attenuated for each different wavelength component.

[0478] In accordance with another embodiment of the present invention, each of the wavefront 2100 and the transformed wavefront 2108 comprises a plurality of different polarization components. In such a case, the plurality of different spatial amplitude changes may be effected by applying an amplitude change to each of the plurality of different polarization components of the transformed wavefront. It is appreciated that the amplitude changes may be spatially different or that the amplitude may be differently attenuated for each different polarization component.

[0479] In accordance with another embodiment of the present invention, the transform 2106 applied to the wavefront 2100 is a Fourier transform, the plurality of different spatial amplitude changes comprise at least three different amplitude changes, effected by applying a spatially uniform, time-varying spatial amplitude attenuation to part of the transformed wavefront 2108, and the plurality of intensity maps 2130, 2132 and 2134 comprises at least three intensity maps. In such a case, employing the plurality of intensity maps to obtain an output indicating the amplitude and phase of the wavefront being analyzed preferably includes:

[0480] expressing the wavefront being analyzed 2100 as a first complex function which has an amplitude and phase identical to the amplitude and phase of the wavefront being analyzed;

[0481] expressing the plurality of intensity maps as a function of the first complex function and of a spatial function governing the spatially uniform, time-varying spatial amplitude change;

[0482] defining a second complex function, having an absolute value and a phase, as a convolution of the first complex function and of a Fourier transform of the spatial function governing the spatially uniform, time-varying spatial amplitude attenuation;

[0483] expressing each of the plurality of intensity maps as a third function of:

[0484] the amplitude of the wavefront being analyzed;

[0485] the absolute value of the second complex function;

[0486] a difference between the phase of the wavefront being analyzed and the phase of the second complex function; and

[0487] a known amplitude attenuation produced by one of the at least three different amplitude changes, to each of which one of the at least three intensity maps corresponds;

[0488] solving the third function to obtain the amplitude of the wavefront being analyzed, the absolute value of the second complex function and the difference between the phase of the wavefront being analyzed and the phase of the second complex function;

[0489] solving the second complex function to obtain the phase of the second complex function; and

[0490] obtaining the phase of the wavefront being analyzed by adding the phase of the second complex function to the difference between the phase of the wavefront being analyzed and the phase of the second complex function.

[0491] Reference is now made to FIG. 21B, which is a simplified partially schematic, partially block diagram illustration of a wavefront analysis system suitable for carrying out the functionality of FIG. 21A in accordance with a preferred embodiment of the present invention. As seen in FIG. 21B, a wavefront, here designated by reference numeral 2150 is focused, as by a lens 2152, onto an amplitude attenuator 2154, which is preferably located at the focal plane of lens 2152. The amplitude attenuator 2154 generates amplitude changes, such as amplitude attenuation, and may be, for example, a spatial light modulator or a series of different partially transparent objects.

[0492] A second lens 2156 is arranged so as to image wavefront 2150 onto a detector 2158, such as a CCD detector. Preferably the second lens 2156 is arranged such that the detector 2158 lies in its focal plane. The output of detector 2158 is preferably supplied to data storage and processing circuitry 2160, which preferably carries out functionality “C” described hereinabove with reference to FIG. 21A.

[0493] Reference is now made to FIG. 22, which is a simplified partially schematic, partially pictorial illustration of a system for surface mapping employing the functionality and structure of FIGS. 21A and 21B. As seen in FIG. 22, a beam of radiation, such as light or acoustic energy, is supplied from a radiation source 2200, optionally via a beam expander 2202, onto a beam splitter 2204, which reflects at least part of the radiation onto a surface 2206 to be inspected. The radiation reflected from the inspected surface, is a surface mapping wavefront, which has an amplitude and a phase, and which contains information about the surface 2206. At least part of the radiation incident on surface 2206 is reflected from the surface 2206 and transmitted via the beam splitter 2204 and focused via a focusing lens 2208 onto an amplitude attenuator 2210, which is preferably located at the image plane of radiation source 2200.

[0494] The amplitude attenuator 2210 may be, for example, a spatial light modulator or a series of different partially transparent non-spatially uniform objects. It is appreciated that amplitude attenuator 2210 can be configured such that a substantial part of the radiation focused thereonto is reflected therefrom. Alternatively the amplitude attenuator 2210 can be configured such that a substantial part of the radiation focused thereonto is transmitted therethrough.

[0495] A second lens 2212 is arranged so as to image surface 2206 onto a detector 2214, such as a CCD detector. Preferably the second lens 2212 is arranged such that the detector 2214 lies in its focal plane. The output of detector 2214, an example of which is a set of intensity maps designated by reference numeral 2215, is preferably supplied to data storage and processing circuitry 2216, which preferably carries out functionality “C” described hereinabove with reference to FIG. 21A, providing an output indicating at least one and possibly both of the phase and the amplitude of the surface mapping wavefront. This output is preferably further processed to obtain information about the surface 2206, such as geometrical variations and reflectivity of the surface.

[0496] In accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 2200 has a narrow wavelength band about a given central wavelength, causing the phase of the radiation reflected from surface 2206 to be proportional to geometrical variations in the surface 2206, the proportion being an inverse linear function of the central wavelength of the radiation.

[0497] In accordance with an embodiment of the present invention, the beam of radiation supplied from radiation source 2200 has at least two narrow wavelength bands, each centered about a different wavelength, designated λ₁, . . . , λ_(n). In such a case, the radiation reflected from the surface 2206 has at least two wavelength components, each centered around a wavelength, λ₁, . . . , λ_(n).

[0498] At least two indications of the phase of the surface mapping wavefront are obtained. Each such indication corresponds to a different wavelength component of the reflected radiation. These at least two indications may be subsequently combined to enable enhanced mapping of the surface 2206, by avoiding ambiguity in the mapping, known as 2π ambiguity, when the value of the mapping at a given spatial location in the surface exceeds the value of the mapping at a different spatial location in the surface by the largest of the different wavelengths λ₁, . . . , λ_(n). A proper choice of the wavelengths λ₁, . . . , λ_(n), may lead to elimination of this ambiguity when the difference in values of the mapping at different locations is smaller than the multiplication product of all the wavelengths.

[0499] In accordance with a preferred embodiment of the present invention, the amplitude attenuator 2210 applies a plurality of different spatial amplitude changes to the radiation wavefront reflected from surface 2206 and Fourier transformed by lens 2208. Application of the plurality of different spatial amplitude changes provides a plurality of differently amplitude changed transformed wavefronts which may be subsequently detected by detector 2214.

[0500] In accordance with yet another preferred embodiment of the present invention, at least three different spatial amplitude changes are applied by amplitude attenuator 2210, resulting in at least three different intensity maps 2215. The at least three intensity maps are employed by the data storage and processing circuitry 2216 to obtain an output indicating at least one and possibly both of the phase and amplitude of the surface mapping wavefront. In such a case, the data storage and processing circuitry 2216, carries out functionality “C” described hereinabove with reference to FIG. 21A, where the wavefront being analyzed (FIG. 21A) is the surface mapping wavefront.

[0501] Additionally, in accordance with a preferred embodiment of the present invention, the beam of radiation supplied from radiation source 2200 comprises a plurality of different wavelength components, thereby providing a plurality of wavelength components in the surface mapping wavefront and subsequently in the transformed wavefront impinging on amplitude attenuator 2210. In this case the amplitude attenuator may be an object, at least one of whose reflection and transmission varies spatially. This spatial variance of the amplitude attenuator generates a different spatial amplitude change for each of the wavelength components, thereby providing a plurality of differently amplitude changed transformed wavefronts to be subsequently detected by detector 2214. It is appreciated that the amplitude attenuation generated by attenuator 2210 may be different for each of the different wavelength components.

[0502] In accordance with an embodiment of the present invention, the surface 2206 is a surface of media in which information is encoded by selecting the height of the media at each of a multiplicity of different locations on the media. In such a case, the indications of the amplitude and phase of the surface mapping wavefront provided by data storage and processing circuitry 2216 are employed to obtain the information encoded on the media.

[0503] It is appreciated that other applications, such as those described hereinabove with respect to FIGS. 16-20 may also be provided in accordance with the present invention wherein amplitude attenuation is performed instead of phase manipulation. It is further appreciated that all of the applications described hereinabove with reference to FIGS. 15-20 may also be provided in accordance with the present invention wherein both amplitude attenuation and phase manipulation are performed.

[0504] It will be appreciated by persons skilled in the art that the present invention is not limited by what has been particularly shown and described hereinabove. Rather the present invention includes both combinations and subcombinations of features described hereinabove as well as modifications and variations of such features which would occur to a person of ordinary skill in the art upon reading the foregoing description and which are not in the prior art. 

1. A method of wavefront analysis comprising: obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed which has an amplitude and a phase; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed.
 2. A method of wavefront analysis according to claim 1 and wherein said plurality of intensity maps are employed to provide an analytical output indicating said amplitude and phase.
 3. A method of wavefront analysis according to claim 1 and wherein said plurality of differently phase changed transformed wavefronts are obtained by interference of said wavefront being analyzed along a common optical path.
 4. A method of wavefront analysis according to claim 1 and wherein said plurality of differently phase changed transformed wavefronts are realized in a manner substantially different from performing a delta-function phase change to said wavefront following the transforming thereof.
 5. A method of wavefront analysis according to claim 1 and wherein said plurality of intensity maps are employed to obtain an output indicating said phase which is substantially free from halo and shading off distortions.
 6. A method according to claim 1 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts resulting from at least one of application of spatial phase changes to a transformed wavefront and transforming of a wavefront following application of spatial phase changes thereto.
 7. A method according to claim 1 and wherein obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said wavefront being analyzed thereby to obtain a transformed wavefront; and applying a plurality of different phase changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 8. A method according to claim 7 and wherein said plurality of different phase changes includes spatial phase changes.
 9. A method according to claim 8 and wherein said plurality of different spatial phase changes are effected by applying a time-varying spatial phase change to part of said transformed wavefront.
 10. A method according to claim 8 and wherein said plurality of different spatial phase changes are effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront.
 11. A method according to claim 10, wherein said transform applied to said wavefront being analyzed is a Fourier transform and wherein said obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts includes applying a Fourier transform to said plurality of differently phase changed transformed wavefronts.
 12. A method according to claim 10 and wherein: said transform applied to said wavefront being analyzed is a Fourier transform; said plurality of different spatial phase changes comprises at least three different phase changes; said plurality of intensity maps comprises at least three intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes: expressing said wavefront being analyzed as a first complex function which has an amplitude and phase identical to said amplitude and phase of said wavefront being analyzed; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; defining a second complex function, having an absolute value and a phase, as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes which each correspond to one of said at least three intensity maps; solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said wavefront being analyzed by adding said phase of said second complex function to said difference between said phase of said wavefront being analyzed and said phase of said second complex function.
 13. A method according to claim 12 and wherein said absolute value of said second complex function is obtained by approximating said absolute value to a polynomial of a given degree.
 14. A method according to claim 12 and wherein said phase of said second complex function is obtained by expressing said second complex function as an eigen-value problem where the complex function is an eigen-vector obtained by an iterative process.
 15. A method according to claim 12 and wherein said phase of said second complex function is obtained by functionality including: approximating said Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change to a polynomial; and approximating said second complex function to a polynomial.
 16. A method according to claim 12 and wherein said amplitude of said wavefront being analyzed, said absolute value of said second complex function, and said difference between said phase of said second complex function and said phase of said wavefront being analyzed are obtained by a least-square method, which has increased accuracy as the number of said plurality of intensity maps increases.
 17. A method according to claim 12 and wherein: said plurality of different phase changes comprises at least four different phase changes; said plurality of intensity maps comprises at least four intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes: expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; a known phase delay produced by one of said at least four different phase changes which each correspond to one of said at least four intensity maps; and at least one additional unknown relating to said wavefront analysis, where the number of said at least one additional unknown is no greater than the number by which said plurality intensity maps exceeds three; and solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function, said difference between said phase of said wavefront being analyzed and said phase of said second complex function and said at least one additional unknown.
 18. A method according to claim 12 and wherein said phase changes are chosen as to maximize contrast in said intensity maps and to minimize effects of noise on said phase of said wavefront being analyzed.
 19. A method according to claim 12 and wherein: expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes which each correspond to one of said at least three intensity maps comprises: defining fourth, fifth and sixth complex functions, none of which being a function of any of said plurality of intensity maps or of said time-varying spatial phase change, each of said fourth, fifth and sixth complex functions being a function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; and expressing each of said plurality of intensity maps as a sum of said fourth complex function, said fifth complex function multiplied by the sine of said known phase delay corresponding to each one of said plurality of intensity maps and said sixth complex function multiplied by the cosine of said known phase delay corresponding to each one of said plurality of intensity maps.
 20. A method according to claim 12 and wherein solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function includes: obtaining two solutions for each of said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function, said two solutions being a higher value solution and a lower value solution; combining said two solutions into an enhanced absolute value solution for said absolute value of said second complex function, by choosing at each spatial location either the higher value solution or the lower value solution of said two solutions in a way that said enhanced absolute value solution satisfies said second complex function; and combining said two solutions of said amplitude of said wavefront being analyzed into enhanced amplitude solution, by choosing at each spatial location the higher value solution or the lower value solution of said two solutions of said amplitude in said way that at each location where said higher value solution is chosen for said absolute value solution, said higher value solution is chosen for said amplitude solution and at each location where said lower value solution is chosen for said absolute value solution, said lower value solution is chosen for said amplitude solution; and combining said two solutions of said difference between said phase of said wavefront being analyzed and said phase of said second complex function into an enhanced difference solution, by choosing at each spatial location the higher value solution or the lower value solution of said two solutions of said difference in said way that at each location where said higher value solution is chosen for said absolute value solution, said higher value solution is chosen for said difference solution and at each location where said lower value solution is chosen for said absolute value solution, said lower value solution is chosen for said difference solution.
 21. A method according to claim 10 and wherein said spatially uniform, time-varying spatial phase change is applied to a spatially central part of said transformed wavefront.
 22. A method according to claim 21 and wherein said transform applied to said wavefront being analyzed is a Fourier transform and wherein said obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts includes applying a Fourier transform to said plurality of differently phase changed transformed wavefronts.
 23. A method according to claim 21 and also comprising: adding a phase component comprising relatively high frequency components to said wavefront being analyzed prior to applying said transform thereto in order to increase the high-frequency content of said transformed wavefront prior to said applying said spatially uniform, time-varying spatial phase change to part of said transformed wavefront.
 24. A method according to claim 10 and wherein said spatially uniform, time-varying spatial phase change is applied to a spatially centered generally circular region of said transformed wavefront.
 25. A method according to claim 10 and wherein said spatially uniform, time-varying spatial phase change is applied to approximately one half of said transformed wavefront.
 26. A method according to claim 10 and wherein: said transformed wavefront includes a DC region and a non-DC region; and said spatially uniform, time-varying spatial phase change is applied to at least part of both said DC region and said non-DC region.
 27. A method according to claim 1 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by employing an at least time varying phase change function.
 28. A method according to claim 1 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by applying an at least time varying phase change function to said wavefront being analyzed.
 29. A method according to claim 28 and wherein said at least time varying phase change function is applied to said wavefront being analyzed prior to transforming thereof.
 30. A method according to claim 28 and wherein said at least time varying phase change function is applied to said wavefront being analyzed subsequent to transforming thereof.
 31. A method according to claim 6 and wherein said plura11ty of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by employing an at least time varying phase change function.
 32. A method according to claim 6 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by applying an at least time varying phase change function to said wavefront to be analyzed.
 33. A method according to claim 32 and wherein said at least time varying phase change function is applied to said wavefront to be analyzed prior to transforming thereof.
 34. A method according to claim 33 and wherein said at least time varying phase change function is a spatially uniform spatial function.
 35. A method according to claim 32 and wherein said at least time varying phase change function is applied to said wavefront to be analyzed subsequent to transforming thereof.
 36. A method according to claim 8 and wherein: said transformed wavefront comprises a plurality of different wavelength components; and said plurality of different spatial phase changes are effected by applying a phase change to said plurality of different wavelength components of said transformed wavefront.
 37. A method according to claim 36 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is a time-varying spatial phase change.
 38. A method according to claim 36 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by passing said transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.
 39. A method according to claim 36 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by reflecting said transformed wavefront from a spatially varying surface.
 40. A method according to claim 36 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is selected to be different to a predetermined extent for at least some of said plurality of different wavelength components.
 41. A method according to claim 36 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is identical for at least some of said plurality of different wavelength components.
 42. A method according to claim 40 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by passing said transformed wavefront through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 43. A method according to claim 41 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by passing said transformed wavefront through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 44. A method according to claim 1 and wherein: said wavefront being analyzed comprises a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to said plurality of different wavelength components of said wavefront being analyzed.
 45. A method according to claim 44 and wherein said phase change is applied to said plurality of different wavelength components of said wavefront being analyzed prior to transforming thereof.
 46. A method according to claim 44 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through an object, at least one of whose thickness and refractive index varies spatially.
 47. A method according to claim 46 and wherein: said obtaining a plurality of intensity maps is performed simultaneously for all of said plurality of different wavelength components; and said obtaining a plurality of intensity maps includes dividing said plurality of phase changed transformed wavefronts into separate wavelength components.
 48. A method according to claim 47 and wherein said dividing said plurality of phase changed transformed wavefronts is effected by passing said plurality of phase changed transformed wavefronts through a dispersion element.
 49. A method according to claim 46 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through an object, at least one of whose thickness and refractive index varies spatially, following transforming of said wavefront being analyzed.
 50. A method according to claim 44 and wherein said phase change applied to said plurality of different wavelength components is effected by reflecting said wavefront being analyzed from a spatially varying surface.
 51. A method according to claim 50 and wherein said phase change applied to said plurality of different wavelength components is effected by reflecting said wavefront being analyzed from a spatially varying surface, following transforming of said wavefront being analyzed.
 52. A method according to claim 44 and wherein said phase change applied to said plurality of different wavelength components is selected to be different to a predetermined extent for at least some of said plurality of different wavelength components.
 53. A method according to claim 44 and wherein said phase change applied to said plurality of different wavelength components is identical for at least some of said plurality of different wavelength components.
 54. A method according to claim 52 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 55. A method according to claim 54 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially, following transforming of said wavefront being analyzed.
 56. A method according to claim 53 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 57. A method according to claim 56 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially, following transforming of said wavefront being analyzed.
 58. A method according to claim 1 and wherein: said wavefront being analyzed comprises a plurality of different polarization components; and said plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to said plurality of different polarization components of said wavefront being analyzed prior to transforming thereof.
 59. A method according to claim 8 and wherein: said transformed wavefront comprises a plurality of different polarization components; and said plurality of different spatial phase changes are effected by applying a phase change to said plurality of different polarization components of said transformed wavefront.
 60. A method according to claim 59 and wherein said phase change applied to said plurality of different polarization components of said transformed wavefront is different for at least some of said plurality of different polarization components.
 61. A method according to claim 59 and wherein said phase change applied to said plurality of different polarization components of said transformed wavefront is identical for at least some of said plurality of different polarization components.
 62. A method of wavefront analysis according to claim 7 and wherein obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts includes: applying a transform to said plurality of differently phase changed transformed wavefronts.
 63. A method according to claim 62 and wherein said plurality of phase changed transformed wavefronts are reflected from a reflecting surface so that said transform applied to said plurality of differently phase changed transformed wavefronts is identical to said transform applied to said wavefront being analyzed.
 64. A method according to claim 7 and wherein said transform applied to said wavefront being analyzed is a Fourier transform.
 65. A method according to claim 1 and wherein said plurality of intensity maps are obtained by reflecting said plurality of differently phase changed transformed wavefronts from a reflecting surface so as to transform said plurality of differently phase changed transformed wavefronts.
 66. A method of wavefront analysis according to claim 1 and wherein obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts includes: applying a transform to said plurality of differently phase changed transformed wavefronts.
 67. A method of wavefront analysis according to claim 1 and wherein employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes: expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said wavefront being analyzed; and employing said at least one mathematical function to obtain an output indicating said phase and amplitude.
 68. A method of wavefront analysis according to claim 7 and wherein employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes: expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said wavefront being analyzed and of said plurality of different phase changes, wherein said phase and amplitude are unknowns and said plurality of different phase changes are known; and employing said at least one mathematical function to obtain an output indicating said phase and amplitude.
 69. A method of wavefront analysis according to claim 1 and wherein: said plurality of intensity maps comprises at least four intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes employing a plurality of combinations, each of at least three of said plurality of intensity maps, to provide a plurality of indications of said amplitude and phase of said wavefront being analyzed.
 70. A method of wavefront analysis according to claim 69 and also comprising employing said plurality of indications of said amplitude and phase of said wavefront being analyzed to provide an enhanced indication of said amplitude and phase of said wavefront being analyzed.
 71. A method of wavefront analysis according to claim 69 and wherein at least some of said plurality of indications of said amplitude and phase are at least second order indications of said amplitude and phase of said wavefront being analyzed.
 72. A method according to claim 1 and wherein obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said wavefront being analyzed, thereby to obtain a transformed wavefront; and applying a plurality of different phase and amplitude changes to said transformed wavefront, thereby to obtain a plurality of differently phase and amplitude changed transformed wavefronts.
 73. A method according to claim 72 and wherein: said plurality of different phase and amplitude changes comprises at least three different phase and intensity changes; said plurality of different phase and amplitude changes are effected by applying at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial amplitude change to at least part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes: expressing said wavefront being analyzed as a first complex function which has an amplitude and phase identical to said amplitude and phase of said wavefront being analyzed; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial amplitude change; defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; and a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and said spatial function governing at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial amplitude change, comprising: defining fourth, fifth, sixth and seventh complex functions, none of which being a function of any of said plurality of intensity maps or of said time-varying spatial phase change, each of said fourth, fifth, sixth and seventh complex functions being a function of at least one of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; defining an eighth function of a phase delay and of an amplitude change, both produced by one of said at least three different phase and amplitude changes, corresponding to said at least three intensity maps; and expressing each of said plurality of intensity maps as a sum of said fourth complex function, said fifth complex function multiplied by the absolute value squared of said eighth function, said sixth complex function multiplied by said eighth function and said seventh complex function multiplied by the complex conjugate of said eighth function; solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said wavefront being analyzed by adding said phase of said second complex function to said difference between said phase of said wavefront being analyzed and phase of said second complex function.
 74. A method according to claim 1 and wherein: said wavefront being analyzed comprises at least two wavelength components; said obtaining a plurality of intensity maps also includes dividing said phase changed transformed wavefronts according to said at least two wavelength components in order to obtain at least two wavelength components of said phase changed transformed wavefronts and in order to obtain at least two sets of intensity maps, each set corresponding to a different one of said at least two wavelength components of said phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said wavefront being analyzed includes obtaining an output indicative of the phase of said wavefront being analyzed from each of said at least two sets of intensity maps and combining said outputs to provide an enhanced indication of phase of said wavefront being analyzed, in which enhanced indication, there is no 2π ambiguity.
 75. A method according to claim 1 and wherein said wavefront being analyzed is an acoustic radiation wavefront.
 76. A method according to claim 7 and wherein: said wavefront being analyzed comprises at least one one-dimensional component; said transform applied to said wavefront being analyzed is a one-dimensional Fourier transform, performed in a dimension perpendicular to a direction of propagation of said wavefront being analyzed, thereby to obtain at least one one-dimensional component of said transformed wavefront in said dimension perpendicular to said direction of propagation; said plurality of differently phase changed transformed wavefronts are obtained by applying said plurality of different phase changes to each of said at least one one-dimensional component, thereby obtaining at least one one-dimensional component of said plurality of phase changed transformed wavefronts; and said plurality of intensity maps are employed to obtain an output indicating amplitude and phase of said at least one one-dimensional component of said wavefront being analyzed.
 77. A method according to claim 76 and wherein said plurality of different phase changes is applied to each of said at least one one-dimensional component by providing a relative movement between said wavefront being analyzed and an element, which element generates spatially varying, time-constant phase changes, said relative movement being in an additional dimension which is perpendicular both to said direction of propagation and to said dimension perpendicular to said direction of propagation.
 78. A method according to claim 76 and wherein: said wavefront being analyzed comprises a plurality of different wavelength components; said plurality of different phase changes are applied to said plurality of different wavelength components of each of said plurality of one-dimensional components of said wavefront being analyzed; and said obtaining a plurality of intensity maps includes dividing said plurality of one-dimensional components of said plurality of phase changed transformed wavefronts into separate wavelength components.
 79. A method according to claim 78 and wherein: said dividing said plurality of one-dimensional components of said plurality of phase changed transformed wavefronts into separate wavelength components is achieved by passing said plurality of phase changed transformed wavefronts through a dispersion element.
 80. A method according to claim 76 and wherein said transform applied to said wavefront being analyzed includes an additional Fourier transform to minimize cross-talk between different one-dimensional components of said wavefront being analyzed.
 81. A method of surface mapping comprising: obtaining a surface mapping wavefront having an amplitude and a phase, by reflecting radiation from a surface; and analyzing said surface mapping wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to said surface mapping wavefront; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said surface mapping wavefront.
 82. A method according to claim 81 and wherein said radiation reflected from said surface has a narrow band about a given wavelength, causing said phase of said surface mapping wavefront to be proportional to geometrical variations in said surface, said proportion being an inverse linear function of said wavelength.
 83. A method according to claim 81 and wherein said radiation reflected from said surface has at least two narrow bands, each centered about a different wavelength, providing at least two wavelength components in said surface mapping wavefront and at least two indications of said phase of said surface mapping wavefront, thereby enabling an enhanced mapping of said surface to be obtained by avoiding an ambiguity in the mapping which exceeds the larger of said different wavelengths about which said two narrow bands are centered.
 84. A method according to claim 81 and wherein obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said surface mapping wavefront, thereby to obtain a transformed wavefront; and applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 85. A method according to claim 84 and wherein: said transform applied to said surface mapping wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said surface mapping wavefront includes: expressing said surface mapping wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said surface mapping wavefront; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said surface mapping wavefront; said absolute value of said second complex function; a difference between said phase of said surface mapping wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; solving said third function to obtain said amplitude of said surface mapping wavefront, said absolute value of said second complex function and said difference between said phase of said surface mapping wavefront and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said surface mapping wavefront by adding said phase of said second complex function to said difference between said phase of said surface mapping wavefront and phase of said second complex function.
 86. A method according to claim 81 and wherein: said surface mapping wavefront comprises a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by: transforming said surface mapping wavefront thereby obtaining a transformed wavefront comprising a plurality of different wavelength components; and applying a phase change to said plurality of different wavelength components of said transformed wavefront by passing said transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.
 87. A method of inspecting an object comprising: obtaining an object inspection wavefront which has an amplitude and a phase, by transmitting radiation through said object; and analyzing said object inspection wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to said object inspection wavefront; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said object inspection wavefront.
 88. A method according to claim 87 and wherein when said object is substantially uniform in material and other optical properties, said phase of said object inspection wavefront is proportional to said object thickness.
 89. A method according to claim 87 and wherein when said object is substantially uniform in thickness, said phase of said object inspection wavefront is proportional to optical properties of said object.
 90. A method according to claim 87 and wherein said radiation has at least two narrow bands, each centered about a different wavelength, providing at least two wavelength components in said object inspection wavefront and at least two indications of said phase of said object inspection wavefront, thereby enabling an enhanced mapping of thickness of said object to be inspected by avoiding an ambiguity in the mapping which exceeds the larger of said different wavelengths about which said two narrow bands are centered.
 91. A method according to claim 87 and wherein obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said object inspection wavefront, thereby to obtain a transformed wavefront; and applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 92. A method according to claim 91 and wherein: said transform applied to said object inspection wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said object inspection wavefront includes: expressing said object inspection wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said object inspection wavefront; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said object inspection wavefront; said absolute value of said second complex function; a difference between said phase of said object inspection wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; solving said third function to obtain said amplitude of said object inspection wavefront, said absolute value of said second complex function and said difference between said phase of said object inspection wavefront and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said object inspection wavefront by adding said phase of said second complex function to said difference between said phase of said object inspection wavefront and phase of said second complex function.
 93. A method according to claim 87 and wherein: said object inspection wavefront comprises a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by: transforming said object inspection wavefront thereby obtaining a transformed wavefront comprising a plurality of different wavelength components; and applying a phase change to said plurality of different wavelength components of said transformed wavefront by reflecting said transformed wavefront from a spatially varying surface.
 94. A method of spectral analysis comprising: obtaining a spectral analysis wavefront having an amplitude and a phase, by causing radiation to impinge on an object; analyzing said spectral analysis wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to said spectral analysis wavefront which has an amplitude and a phase; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said spectral analysis wavefront; and employing said output indicating said amplitude and phase to obtain an output indicating spectral content of said radiation.
 95. A method according to claim 94 and wherein obtaining said spectral analysis wavefront is effected by reflecting said radiation from said object.
 96. A method according to claim 94 and wherein obtaining said spectral analysis wavefront is effected by transmitting said radiation through said object.
 97. A method according to claim 94 and wherein when said radiation is substantially of a single wavelength, said phase of said spectral analysis wavefront is inversely proportional to said single wavelength, and is related to at least one of a surface characteristic and thickness of said impinged object.
 98. A method according to claim 94 and wherein employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said spectral analysis wavefront includes: expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said spectral analysis wavefront and of said plurality of different phase changes, wherein at least said phase is unknown and a function generating said plurality of phase changed transformed wavefronts is known; and employing said at least one mathematical function to obtain an output indicating at least said phase.
 99. A method according to claim 94 and wherein obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said spectral analysis wavefront, thereby to obtain a transformed wavefront; and applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 100. A method according to claim 99 and wherein: said transform applied to said spectral analysis wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said spectral analysis wavefront includes: expressing said spectral analysis wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said spectral analysis wavefront; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said spectral analysis wavefront; said absolute value of said second complex function; a difference between said phase of said spectral analysis wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; solving said third function to obtain said amplitude of said spectral analysis wavefront, said absolute value of said second complex function and said difference between said phase of said spectral analysis wavefront and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said spectral analysis wavefront by adding said phase of said second complex function to said difference between said phase of said spectral analysis wavefront and phase of said second complex function.
 101. A method according to claim 94 and wherein: said spectral analysis wavefront comprises a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to said plurality of different wavelength components of said spectral analysis wavefront.
 102. A method of phase change analysis comprising: obtaining a phase change analysis wavefront which has an amplitude and a phase; applying a transform to said phase change analysis wavefront thereby to obtain a transfomied wavefront; applying a plurality of different phase changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts. obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indication of differences between said plurality of different phase changes applied to said transformed phase change analysis wavefront.
 103. A method according to claim 102 and wherein when lateral shifts appear in said plurality of different phase changes, corresponding changes appear in said plurality of intensity maps, said employing results in obtaining an indication of said lateral shifts.
 104. A method according to claim 102 and wherein employing said plurality of intensity maps to obtain an output indication of differences between said plurality of different phase changes applied to said transformed phase change analysis wavefront includes: expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said phase change analysis wavefront and of said plurality of different phase changes, where at least said phase and amplitude are known and said plurality of different phase changes are unknown; and employing said at least one mathematical function to obtain an output indicating said differences between said plurality of different phase changes.
 105. A method of phase change analysis comprising: obtaining a phase change analysis wavefront which has an amplitude and a phase; applying a transform to said phase change analysis wavefront thereby to obtain a transformed wavefront; applying at least one phase change to said transformed wavefront, thereby to obtain at least one phase changed transformed wavefront. obtaining at least one intensity map of said at least one phase changed transformed wavefront; and employing said at least one intensity map to obtain an output indication of said at least one phase change applied to said transformed phase change analysis wavefront.
 106. A method according to claim 105 and wherein said at least one phase change is a phase delay, having a value selected from a plurality of pre-determined values, and said output indication of said at least one phase change includes said value of said phase delay.
 107. A method of stored data retrieval comprising: obtaining a stored data retrieval wavefront which has an amplitude and a phase, by reflecting radiation from media in which information is encoded by selecting the height of the media at each of a multiplicity of different locations on the media; analyzing said stored data retrieval wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to said stored data retrieval wavefront; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an indication of said amplitude and phase of said stored data retrieval wavefront; and employing said indication of said amplitude and phase to obtain said information.
 108. A method according to claim 107 and wherein said obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said stored data retrieval wavefront thereby to obtain a transformed wavefront; and applying a plurality of different phase changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 109. A method according to claim 108 and wherein: said transform applied to said stored data retrieval wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said stored data retrieval wavefront includes: expressing said stored data retrieval wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said stored data retrieval wavefront; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said stored data retrieval wavefront; said absolute value of said second complex function; a difference between said phase of said stored data retrieval wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; solving said third function to obtain said amplitude of said stored data retrieval wavefront, said absolute value of said second complex function and said difference between said phase of said stored data retrieval wavefront and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said stored data retrieval wavefront by adding said phase of said second complex function to said difference between said phase of said stored data retrieval wavefront and phase of said second complex function.
 110. A method according to claim 108 and wherein: said stored data retrieval wavefront comprises at least one one-dimensional component; said transform applied to said data retrieval wavefront is a one-dimensional Fourier transform, performed in a dimension perpendicular to a direction of propagation of said data retrieval wavefront, thereby to obtain at least one one-dimensional component of said transformed wavefront in said dimension perpendicular to said direction of propagation; said plurality of differently phase changed transformed wavefronts are obtained by applying said plurality of different phase changes to each of said at least one one-dimensional component, thereby obtaining at least one one-dimensional component of said plurality of phase changed transformed wavefronts; and said plurality of intensity maps are employed to obtain an output indicating amplitude and phase of said at least one one-dimensional component of said data retrieval wavefront.
 111. A method according to claim 110 and wherein said plurality of different phase changes is applied to each of said at least one one-dimensional component by providing a relative movement between said media and a component generating spatially varying, time-constant phase changes, said relative movement being in a dimension perpendicular to said direction of propagation and to said dimension perpendicular to said direction of propagation.
 112. A method according to claim 107 and wherein said information is encoded on said media whereby: an intensity value is realized by reflection of light from each location on said media to lie within a predetermined range of values, said range corresponding an element of said information stored at said location; and by employing said plurality of intensity maps, multiple intensity values are realized for each location, providing multiple elements of information for each location on said media.
 113. A method according to claim 112 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by applying an at least time varying phase change function to said stored data retrieval wavefront.
 114. A method according to claim 112 and wherein: said stored data retrieval wavefront comprises a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by applying at least one phase change to said plurality of different wavelength components of said stored data retrieval wavefront.
 115. A method according to claim 107 and wherein: said radiation which is reflected from said media comprises a plurality of different wavelength components, resulting in said stored data retrieval wavefront comprising a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by applying a phase change to said plurality of different wavelength components of said stored data retrieval wavefront.
 116. A method according to claim 107 and wherein: said information encoded by selecting the height of the media at each of a multiplicity of different locations on the media is also encoded by selecting the reflectivity of the media at each of a plurality of different locations on the media; and employing said indication of said amplitude and phase to obtain said information includes employing said indication of said phase to obtain said information encoded by selecting the height of the media and employing said indication of said amplitude to obtain said information encoded by selecting the reflectivity of the media.
 117. A method of 3-dimensional imaging comprising: obtaining a 3-dimensional imaging wavefront, which has an amplitude and a phase, by reflecting radiation from an object to be viewed; and analyzing said 3-dimensional imaging wavefront by: obtaining a plurality of differently phase changed transformed wavefronts corresponding to said 3-dimensional imaging wavefront; obtaining a plurality of intensity maps of said plurality of differently phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating said amplitude and phase of said 3-dimensional imaging wavefront.
 118. A method according to claim 117 and wherein said radiation reflected from said object has a narrow band about a given wavelength, causing said phase of said 3-dimensional imaging wavefront to be proportional to geometrical variations in said object, said proportion being an inverse linear function of said wavelength.
 119. A method according to claim 117 and wherein obtaining a plurality of differently phase changed transformed wavefronts comprises: applying a transform to said 3-dimensional imaging wavefront, thereby to obtain a transformed wavefront; and applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 120. A method according to claim 117 and wherein: said 3-dimensional imaging wavefront comprises a plurality of different wavelength components; and said plurality of differently phase changed transformed wavefronts are obtained by: transforming said 3-dimensional imaging wavefront, thereby obtaining a transformed wavefront comprising a plurality of different wavelength components; and applying phase changes to said plurality of different wavelength components of said transformed wavefront by passing said transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.
 121. A method of wavefront analysis comprising: obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating at least phase of said wavefront being analyzed by combining said plurality of intensity maps into a second plurality of combined intensity maps, said second plurality being less than said first plurality, obtaining at least an output indicative of said phase of said wavefront being analyzed from each of said second plurality of combined intensity maps and combining said outputs to provide at least an enhanced indication of phase of said wavefront being analyzed.
 122. A method of wavefront analysis comprising: obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating at least amplitude of said wavefront being analyzed by combining said plurality of intensity maps into a second plurality of combined intensity maps, said second plurality being less than said first plurality, obtaining at least an output indicative of said amplitude of said wavefront being analyzed from each of said second plurality of combined intensity maps and combining said outputs to provide at least an enhanced indication of amplitude of said wavefront being analyzed.
 123. A method of wavefront analysis comprising: obtaining a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed; obtaining a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and employing said plurality of intensity maps to obtain an output indicating at least phase of said wavefront being analyzed by: expressing said plurality of intensity maps as a function of: amplitude of said wavefront being analyzed; phase of said wavefront being analyzed; and a phase change function characterizing said plurality of differently phase changed transformed wavefronts; defining a complex function of: said amplitude of said wavefront being analyzed; said phase of said wavefront being analyzed; and said phase change function characterizing said plurality of differently phase changed transformed wavefronts, said complex function being characterized in that intensity at each location in said plurality of intensity maps is a function predominantly of a value of said complex function at said location and of said amplitude and said phase of said wavefront being analyzed at said location; expressing said complex function as a function of said plurality of intensity maps; and obtaining values for said phase by employing said complex function expressed as a function of said plurality of intensity maps.
 124. A method of wavefront analysis comprising: applying a Fourier transform to a wavefront being analyzed which has an amplitude and a phase thereby to obtain a transformed wavefront; applying a spatially uniform time-varying spatial phase change to part of said transformed wavefront, thereby to obtain at least three differently phase changed transformed wavefronts; applying a second Fourier transform to obtain at least three intensity maps of said at least three phase changed transformed wavefronts; and employing said at least three intensity maps to obtain an output indicating at least one of said phase and said amplitude of said wavefront being analyzed by: expressing said wavefront being analyzed as a first complex function which has an amplitude and phase identical to said amplitude and phase of said wavefront being analyzed; expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; solving said second complex function to obtain said phase of said second complex function; and obtaining said phase of said wavefront being analyzed by adding said phase of said second complex function to said difference between said phase of said wavefront being analyzed and phase of said second complex function.
 125. Apparatus for wavefront analysis comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed which has an amplitude and a phase; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating said amplitude and phase of said wavefront being analyzed.
 126. Apparatus for wavefront analysis according to claim 125 and wherein said intensity map utilizer employs said plurality of intensity maps to provide an analytical output indicating said amplitude and phase.
 127. Apparatus for wavefront analysis according to claim 125 and wherein said wavefront transformer provides said plurality of differently phase changed transformed wavefronts by interference of said wavefront being analyzed along a common optical path.
 128. Apparatus for wavefront analysis according to claim 125 and wherein said plurality of differently phase changed transformed wavefronts are realized in a manner substantially different from performing a delta-function phase change to said wavefront following transforming thereof.
 129. Apparatus for wavefront analysis according to claim 125 and wherein said intensity map utilizer employs said plurality of intensity maps to provide an output indicating said phase which is substantially free from halo and shading off distortions.
 130. Apparatus according to claim 125 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts resulting from at least one of application of spatial phase changes to a transformed wavefront and transforming of a wavefront following application of spatial phase changes thereto.
 131. Apparatus according to claim 125 and wherein said wavefront transformer comprises: a transform applier, applying a transform to said wavefront being analyzed thereby to obtain a transformed wavefront; and a phase change applier, applying a plurality of different phase changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 132. Apparatus according to claim 131 and wherein said plurality of different phase changes includes spatial phase changes.
 133. Apparatus according to claim 132 and wherein said plurality of different spatial phase changes are effected by applying a time-varying spatial phase change to part of said transformed wavefront.
 134. Apparatus according to claim 132 and wherein said plurality of different spatial phase changes are effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront.
 135. Apparatus according to claim 134, wherein said transform applied to said wavefront being analyzed is a Fourier transform and wherein said intensity map generator includes a Fourier transform applier which applies a Fourier transform to said plurality of differently phase changed transformed wavefronts.
 136. Apparatus according to claim 134 and wherein: said transform applied to said wavefront being analyzed is a Fourier transform; said plurality of different spatial phase changes comprises at least three different phase changes; said plurality of intensity maps comprises at least three intensity maps; and said intensity map utilizer includes: a wavefront expresser, expressing said wavefront being analyzed as a first complex function which has an amplitude and phase identical to said amplitude and phase of said wavefront being analyzed; a first intensity map expresser, expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; a complex function definer, defining a second complex function, having an absolute value and a phase, as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser, expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes which each correspond to one of said at least three intensity maps; a first function solver, solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said wavefront being analyzed by adding said phase of said second complex function to said difference between said phase of said wavefront being analyzed and said phase of said second complex function.
 137. Apparatus according to claim 136 and wherein said first function solver is operative to obtain said absolute value of said second complex function by approximating said absolute value to a polynomial of a given degree.
 138. Apparatus according to claim 136 and wherein said phase obtainer is operative to obtain the phase of said second complex function by expressing said second complex function as an eigen-value problem where the complex function is an eigen-vector obtained by an iterative process.
 139. Apparatus according to claim 136 and wherein said second function solver is operative to obtain the phase of said second complex function by employing functionality including: first approximation functionality approximating said Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change to a polynomial; and second approximation functionality approximating said second complex function to a polynomial.
 140. Apparatus according to claim 136 and wherein said first function solver is operative to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function, and said difference between said phase of said second complex function and said phase of said wavefront being analyzed by a least-square method, which has increased accuracy as the number of said plurality of intensity maps increases.
 141. Apparatus according to claim 136 and wherein: said plurality of different phase changes comprises at least four different phase changes; said plurality of intensity maps comprises at least four intensity maps; and said intensity map utilizer includes: intensity map expressing functionality, expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; a known phase delay produced by one of said at least four different phase changes which each correspond to one of said at least four intensity maps; and at least one additional unknown relating to said wavefront analysis, where the number of said at least one additional unknown is no greater than the number by which said plurality intensity maps exceeds three; and a function solver, solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function, said difference between said phase of said wavefront being analyzed and said phase of said second complex function and said at least one additional unknown.
 142. Apparatus according to claim 136 and wherein said phase changes are chosen as to maximize contrast in said intensity maps and to minimize effects of noise on said phase of said wavefront being analyzed.
 143. Apparatus according to claim 136 and wherein: said second intensity map expresser is operative to express each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes which each correspond to one of said at least three intensity maps comprises: a second complex function definer, defining fourth, fifth and sixth complex functions, none of which being a function of any of said plurality of intensity maps or of said time-varying spatial phase change, each of said fourth, fifth and sixth complex functions being a function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a third intensity map expresser, expressing each of said plurality of intensity maps as a sum of said fourth complex function, said fifth complex function multiplied by the sine of said known phase delay corresponding to each one of said plurality of intensity maps and said sixth complex function multiplied by the cosine of said known phase delay corresponding to each one of said plurality of intensity maps.
 144. Apparatus according to claim 135 and wherein said first function solver includes: function solving functionality, obtaining two solutions for each of said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function, said two solutions being a higher value solution and a lower value solution; first combining functionality, combining said two solutions into an enhanced absolute value solution for said absolute value of said second complex function, by choosing at each spatial location either the higher value solution or the lower value solution of said two solutions in a way that said enhanced absolute value solution satisfies said second complex function; and second combining functionality, combining said two solutions of said amplitude of said wavefront being analyzed into enhanced amplitude solution, by choosing at each spatial location the higher value solution or the lower value solution of said two solutions of said amplitude in said way that at each location where said higher value solution is chosen for said absolute value solution, said higher value solution is chosen for said amplitude solution and at each location where said lower value solution is chosen for said absolute value solution, said lower value solution is chosen for said amplitude solution; and third combining functionality, combining said two solutions of said difference between said phase of said wavefront being analyzed and said phase of said second complex function into an enhanced difference solution, by choosing at each spatial location the higher value solution or the lower value solution of said two solutions of said difference in said way that at each location where said higher value solution is chosen for said absolute value solution, said higher value solution is chosen for said difference solution and at each location where said lower value solution is chosen for said absolute value solution, said lower value solution is chosen for said difference solution.
 145. Apparatus according to claim 134 and wherein said spatially uniform, time-varying spatial phase change is applied to a spatially central part of said transformed wavefront.
 146. Apparatus according to claim 145, wherein said transform applied to said wavefront being analyzed is a Fourier transform and wherein said intensity map generator includes a Fourier transform applier which applies a Fourier transform to said plurality of differently phase changed transformed wavefronts.
 147. Apparatus according to claim 145, and also comprising: a phase adder operative to add a phase component comprising relatively high frequency components to said wavefront being analyzed prior to applying said transform thereto in order to increase the high-frequency content of said transformed wavefront prior to said applying said spatially uniform, time-varying spatial phase change to part of said transformed wavefront.
 148. Apparatus according to claim 134 and wherein said spatially uniform, time-varying spatial phase change is applied to a spatially centered generally circular region of said transformed wavefront.
 149. Apparatus according to claim 134 and wherein said spatially uniform, time-varying spatial phase change is applied to approximately one half of said transformed wavefront.
 150. Apparatus according to claim 134 and wherein: said transformed wavefront includes a DC region and a non-DC region; and said spatially uniform, time-varying spatial phase change is applied to at least part of both said DC region and said non-DC region.
 151. Apparatus according to claim 125 and wherein said wavefront transformer comprises a phase changer operative to change the phase of a plurality of wavefronts by employing an at least time varying phase change function.
 152. Apparatus according to claim 125 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by a phase changer, operative to apply an at least time varying phase change function to said wavefront being analyzed.
 153. Apparatus according to claim 152 and wherein said phase changer is operative to provide an at least time varying phase change function to said wavefront being analyzed prior to transforming thereof.
 154. Apparatus according to claim 152 and wherein said phase changer is operative to provide an at least time varying phase change function to said wavefront being analyzed subsequent to transforming thereof.
 155. Apparatus according to claim 130 and wherein said wavefront transformer comprises a phase changer operative to change the phase of a plurality of wavefronts by employing an at least time varying phase change function.
 156. Apparatus according to claim 130 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by a phase changer, operative to apply an at least time varying phase change function to said wavefront to be analyzed.
 157. Apparatus according to claim 156 and wherein said phase changer is operative to provide an at least time varying phase change function to said wavefront to be analyzed prior to transforming thereof.
 158. Apparatus according to claim 157, and wherein said at least time varying phase change function is a spatially uniform spatial function.
 159. Apparatus according to claim 156 and wherein said phase changer is operative to provide an at least time varying phase change function to said wavefront to be analyzed subsequent to transforming thereof.
 160. Apparatus according to claim 132 and wherein: said transformed wavefront comprises a plurality of different wavelength components; and said plurality of different spatial phase changes are effected by applying a phase change to said plurality of different wavelength components of said transformed wavefront.
 161. Apparatus according to claim 160, and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is a time-varying spatial phase change.
 162. Apparatus according to claim 160 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by passing said transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.
 163. Apparatus according to claim 160 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by reflecting said transformed wavefront from a spatially varying surface.
 164. Apparatus according to claim 160 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is selected to be different to a predetermined extent for at least some of said plurality of different wavelength components.
 165. Apparatus according to claim 160 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is identical for at least some of said plurality of different wavelength components.
 166. Apparatus according to claim 164 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by passing said transformed wavefront through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 167. Apparatus according to claim 165 and wherein said phase change applied to said plurality of different wavelength components of said transformed wavefront is effected by passing said transformed wavefront through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 168. Apparatus according to claim 125 and wherein: said wavefront being analyzed comprises a plurality of different wavelength components; and said wavefront transformer is operative to apply a phase change to said plurality of different wavelength components of said wavefront being analyzed, thereby obtaining said plurality of differently phase changed transformed wavefronts.
 169. Apparatus according to claim 168 and wherein said wavefront transformer is operative to apply said phase change to said plurality of different wavelength components of said wavefront being analyzed prior to transforming thereof.
 170. Apparatus according to claim 168 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through an object, at least one of whose thickness and refractive index varies spatially.
 171. Apparatus according to claim 170, and wherein: said intensity map generator is operative to obtain said plurality of intensity maps simultaneously for all of said plurality of different wavelength components; and said intensity map generator includes a wavelength divider, dividing said plurality of phase changed transformed wavefronts into separate wavelength components.
 172. Apparatus according to claim 171 and wherein said wavelength divider includes a dispersion element, dividing said plurality of phase changed transformed wavefronts passing therethrough into separate wavelength components.
 173. Apparatus according to claim 170 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through an object, at least one of whose thickness and refractive index varies spatially, following transforming of said wavefront being analyzed.
 174. Apparatus according to claim 168 and wherein said phase change applied to said plurality of different wavelength components is effected by reflecting said wavefront being analyzed from a spatially varying surface.
 175. Apparatus according to claim 174 and wherein said phase change applied to said plurality of different wavelength components is effected by reflecting said wavefront being analyzed from a spatially varying surface, following transforming of said wavefront being analyzed.
 176. Apparatus according to claim 168 and wherein said phase change applied to said plurality of different wavelength components is selected to be different to a predetermined extent for at least some of said plurality of different wavelength components.
 177. Apparatus according to claim 168 and wherein said phase change applied to said plurality of different wavelength components is identical for at least some of said plurality of different wavelength components.
 178. Apparatus according to claim 176 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 179. Apparatus according to claim 178 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially, following transforming of said wavefront being analyzed.
 180. Apparatus according to claim 177 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially.
 181. Apparatus according to claim 180 and wherein said phase change applied to said plurality of different wavelength components is effected by passing said wavefront being analyzed through a plurality of objects, each characterized in that at least one of its thickness and refractive index varies spatially, following transforming of said wavefront being analyzed.
 182. Apparatus according to claim 125, and wherein: said wavefront being analyzed comprises a plurality of different polarization components; and said wavefront transformer is operative to apply a phase change to said plurality of different polarization components of said wavefront being analyzed prior to transforming thereof, thereby obtaining said plurality of differently phase changed transformed wavefronts.
 183. Apparatus according to claim 132 and wherein: said transformed wavefront comprises a plurality of different polarization components; and said plurality of different spatial phase changes are effected by applying a phase change to said plurality of different polarization components of said transformed wavefront.
 184. Apparatus according to claim 183 and wherein said phase change applied to said plurality of different polarization components of said transformed wavefront is different for at least some of said plurality of different polarization components.
 185. Apparatus according to claim 183 and wherein said phase change applied to said plurality of different polarization components of said transformed wavefront is identical for at least some of said plurality of different polarization components.
 186. Apparatus for wavefront analysis according to claim 131 and wherein said intensity map generator includes: a second transform applier, applying a transform to said plurality of differently phase changed transformed wavefronts.
 187. Apparatus according to claim 186 and wherein said plurality of phase changed transformed wavefronts are reflected from a reflecting surface so that said transform applier, applying a transform to said wavefront being analyzed and said second transform applier, applying a transform to said plurality of differently phase changed transformed wavefronts, are the same element.
 188. Apparatus according to claim 131 and wherein said transform applier applies a Fourier transform to said wavefront being analyzed.
 189. Apparatus according to claim 125 and wherein said intensity map generator includes a reflecting surface, reflecting said plurality of differently phase changed transformed wavefronts so as to transform said plurality of differently phase changed transformed wavefronts.
 190. Apparatus for wavefront analysis according to claim 125 and wherein said intensity map generator includes a transform applier, applying a transform to said plurality of differently phase changed transformed wavefronts.
 191. Apparatus for wavefront analysis according to claim 125 and wherein said intensity map utilizer includes: an intensity map expresser, expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said wavefront being analyzed; and a function solver, employing said at least one mathematical function to obtain an output indicating said phase and amplitude.
 192. Apparatus for wavefront analysis according to claim 131 and wherein said intensity map utilizer includes: an intensity map expresser, expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said wavefront being analyzed and of said plurality of different phase changes, wherein said phase and amplitude are unknowns and said plurality of different phase changes are known; and a function solver, employing said at least one mathematical function to obtain an output indicating said phase and amplitude.
 193. Apparatus for wavefront analysis according to claim 125 and wherein: said plurality of intensity maps comprises at least four intensity maps; and said intensity map utilizer includes an indication provider, employing a plurality of combinations, each of at least three of said plurality of intensity maps, to provide a plurality of indications of said amplitude and phase of said wavefront being analyzed.
 194. Apparatus for wavefront analysis according to claim 193 and wherein said indication provider also includes an enhanced indication provider, employing said plurality of indications of said amplitude and phase of said wavefront being analyzed to provide an enhanced indication of said amplitude and phase of said wavefront being analyzed.
 195. Apparatus for wavefront analysis according to claim 193 and wherein at least some of said plurality of indications of said amplitude and phase are at least second order indications of said amplitude and phase of said wavefront being analyzed.
 196. Apparatus according to claim 125 and wherein said wavefront transformer comprises: a transform applier, applying a transform to said wavefront being analyzed, thereby to obtain a transformed wavefront; and a phase and intensity change applier, applying a plurality of different phase and intensity changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 197. Apparatus according to claim 196 and wherein: said plurality of different phase and intensity changes comprises at least three different phase and intensity changes; said phase and intensity change applier is operative by applying at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial intensity change to at least part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and said intensity map utilizer includes: a wavefront expresser, expressing said wavefront being analyzed as a first complex function which has an amplitude and phase identical to said amplitude and phase of said wavefront being analyzed; a first intensity map expresser, expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial intensity change; a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser, expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; and a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and said spatial function governing at least one of a spatially uniform, time-varying spatial phase change and a spatially uniform, time-varying spatial intensity change, comprising: a second complex function definer, defining fourth, fifth, sixth and seventh complex functions, none of which being a function of any of said plurality of intensity maps or of said time-varying spatial phase change, each of said fourth, fifth, sixth and seventh complex functions being a function of at least one of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; a third function definer, defining an eighth function of a phase delay and of an intensity change, both produced by one of said at least three different phase and intensity changes, corresponding to said at least three intensity maps; and a third intensity map expresser, expressing each of said plurality of intensity maps as a sum of said fourth complex function, said fifth complex function multiplied by the absolute value squared of said eighth function, said sixth complex function multiplied by said eighth function and said seventh complex function multiplied by the complex conjugate of said eighth function; a first function solver, solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said wavefront being analyzed by adding said phase of said second complex function to said difference between said phase of said wavefront being analyzed and said phase of said second complex function.
 198. Apparatus according to claim 125, and wherein: said wavefront being analyzed comprises at least two wavelength components; said intensity map generator also includes a wavefront divider, dividing said phase changed transformed wavefronts according to said at least two wavelength components thereby obtaining at least two wavelength components of said phase changed transformed wavefronts and subsequently obtaining at least two sets of intensity maps, each set corresponding to a different one of said at least two wavelength components of said phase changed transformed wavefronts; and said intensity map utilizer includes a phase obtainer, obtaining an output indicative of the phase of said wavefront being analyzed from each of said at least two sets of intensity maps and combining said outputs to provide an enhanced indication of phase of said wavefront being analyzed, in which enhanced indication, there is no 2π ambiguity.
 199. Apparatus according to claim 125, and wherein said wavefront being analyzed is an acoustic radiation wavefront.
 200. Apparatus according to claim 131 and wherein: said wavefront being analyzed comprises at least one one-dimensional component; said transform applier is operative to perform a one-dimensional Fourier transform to said wavefront being analyzed is, performed in a dimension perpendicular to a direction of propagation of said wavefront being analyzed, thereby to obtain at least one one-dimensional component of said transformed wavefront in said dimension perpendicular to said direction of propagation; said phase change applier is operative to apply a plurality of different phase changes to each of said at least one one-dimensional component, thereby obtaining at least one one-dimensional component of said plurality of phase changed transformed wavefronts; and said intensity map utilizer is operative to obtain an output indicating amplitude and phase of said at least one one-dimensional component of said wavefront being analyzed.
 201. Apparatus according to claim 200 and wherein said phase change applier comprises a movement generator, providing a relative movement between said wavefront being analyzed and an element, which element generates spatially varying, time-constant phase changes, said relative movement being in an additional dimension which is perpendicular both to said direction of propagation and to said dimension perpendicular to said direction of propagation.
 202. Apparatus according to claim 200 and wherein: said wavefront being analyzed comprises a plurality of different wavelength components; said phase change applier is operative to apply a plurality of different phase changes to said plurality of different wavelength components of each of said plurality of one-dimensional components of said wavefront being analyzed; and said intensity map generator includes a wavelength divider, dividing said plurality of one-dimensional components of said plurality of phase changed transformed wavefronts into separate wavelength components.
 203. Apparatus according to claim 202 and wherein: said wavelength divider includes a dispersion element, dividing said plurality of one-dimensional components of said plurality of phase changed transformed wavefronts passing therethrough into separate wavelength components.
 204. Apparatus according to claim 200 and wherein said transform applier includes an additional transform applier, operative to perform an additional Fourier transform to minimize cross-talk between different one-dimensional components of said wavefront being analyzed.
 205. Apparatus for surface mapping comprising: a wavefront obtainer operative to obtain a surface mapping wavefront having an amplitude and a phase, by reflecting radiation from a surface; and a wavefront analyzer, analyzing said surface mapping wavefront and comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to said surface mapping wavefront; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating said amplitude and phase of said surface mapping wavefront.
 206. Apparatus according to claim 205, and wherein said radiation reflected from said surface has a narrow band about a given wavelength, causing said phase of said surface mapping wavefront to be proportional to geometrical variations in said surface, said proportion being an inverse linear function of said wavelength.
 207. Apparatus according to claim 205, and wherein said radiation reflected from said surface has at least two narrow bands, each centered about a different wavelength, providing at least two wavelength components in said surface mapping wavefront and at least two indications of said phase of said surface mapping wavefront, thereby enabling an enhanced mapping of said surface to be obtained by avoiding an ambiguity in the mapping which exceeds the larger of said different wavelengths about which said two narrow bands are centered.
 208. Apparatus according to claim 205 and wherein said wavefront transformer comprises: a transform applier, applying a transform to said surface mapping wavefront, thereby to obtain a transformed wavefront; and a phase change applier, applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 209. Apparatus according to claim 208, and wherein: said transform applied to said surface mapping wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and said intensity map utilizer includes: a wavefront expresser, expressing said surface mapping wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said surface mapping wavefront; a first intensity map expresser, expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser, expressing each of said plurality of intensity maps as a third function of: said amplitude of said surface mapping wavefront; said absolute value of said second complex function; a difference between said phase of said surface mapping wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; a first function solver, solving said third function to obtain said amplitude of said surface mapping wavefront, said absolute value of said second complex function and said difference between said phase of said surface mapping wavefront and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said surface mapping wavefront by adding said phase of said second complex function to said difference between said phase of said surface mapping wavefront and said phase of said second complex function.
 210. Apparatus according to claim 205 and wherein: said surface mapping wavefront comprises a plurality of different wavelength components; and said wavefront transformer includes: a transform applier, applying a transform to said surface mapping wavefront, thereby to obtain a transformed wavefront comprising a plurality of different wavelength components; and a phase change applier, applying a phase change to said plurality of different wavelength components of said transformed wavefront by passing said transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.
 211. Apparatus for inspecting an object comprising: a wavefront obtainer operative to obtain an object inspection wavefront which has an amplitude and a phase, by transmitting radiation through said object; and a wavefront analyzer, analyzing said object inspection wavefront and comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to said object inspection wavefront; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating said amplitude and phase of said object inspection wavefront.
 212. Apparatus according to claim 211, and wherein when said object is substantially uniform in material and other optical properties, said phase of said object inspection wavefront is proportional to said object thickness.
 213. Apparatus according to claim 211, and wherein when said object is substantially uniform in thickness, said phase of said object inspection wavefront is proportional to optical properties of said object.
 214. Apparatus according to claim 211, and wherein said radiation has at least two narrow bands, each centered about a different wavelength, providing at least two wavelength components in said object inspection wavefront and at least two indications of said phase of said object inspection wavefront, thereby enabling an enhanced mapping of thickness of said object to be inspected by avoiding an ambiguity in the mapping which exceeds the larger of said different wavelengths about which said two narrow bands are centered.
 215. Apparatus according to claim 211, and wherein said wavefront transformer comprises: a transform applier, applying a transform to said object inspection wavefront, thereby to obtain a transformed wavefront; and a phase change applier, applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 216. Apparatus according to claim 215, and wherein: said transform applied to said object inspection wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and said intensity map utilizer includes: a wavefront expresser, expressing said object inspection wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said object inspection wavefront; a first intensity map expresser, expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser, expressing each of said plurality of intensity maps as a third function of: said amplitude of said object inspection wavefront; said absolute value of said second complex function; a difference between said phase of said object inspection wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; a first function solver, solving said third function to obtain said amplitude of said object inspection wavefront, said absolute value of said second complex function and said difference between said phase of said object inspection wavefront and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said object inspection wavefront by adding said phase of said second complex function to said difference between said phase of said object inspection wavefront and said phase of said second complex function.
 217. Apparatus according to claim 211 and wherein: said object inspection wavefront comprises a plurality of different wavelength components; and said wavefront transformer includes: a transform applier, applying a transform to said object inspection wavefront thereby obtaining a transformed wavefront comprising a plurality of different wavelength components; and a phase change applier, applying a phase change to said plurality of different wavelength components of said transformed wavefront by reflecting said transformed wavefront from a spatially varying surface.
 218. Apparatus for spectral analysis comprising: a wavefront obtainer operative to obtain a spectral analysis wavefront having an amplitude and a phase, by causing radiation to impinge on an object; a wavefront analyzer, analyzing said spectral analysis wavefront and comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to said spectral analysis wavefront which has an amplitude and a phase; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating said amplitude and phase of said spectral analysis wavefront; and a phase and amplitude utilizer, employing said output indicating said amplitude and phase to obtain an output indicating spectral content of said radiation.
 219. Apparatus for spectral analysis according to claim 218 and wherein said wavefront obtainer is operative to obtain said spectral analysis wavefront by reflecting said radiation from said object.
 220. Apparatus for spectral analysis according to claim 218 and wherein said wavefront obtainer is operative to obtain said spectral analysis wavefront by transmitting said radiation through said object.
 221. Apparatus for spectral analysis according to claim 218 and wherein when said radiation is substantially of a single wavelength, said phase of said spectral analysis wavefront is inversely proportional to said single wavelength, and is related to at least one of a surface characteristic and thickness of said impinged object.
 222. Apparatus for spectral analysis according to claim 218 and wherein said intensity map utilizer includes: an intensity map expresser, expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said spectral analysis wavefront and of said plurality of different phase changes, wherein at least said phase is unknown and a function generating said plurality of phase changed transformed wavefronts is known; and a function solver, employing said at least one mathematical function to obtain an output indicating at least said phase.
 223. Apparatus according to claim 218 and wherein said wavefront transformer comprises: a transform applier, applying a transform to said spectral analysis wavefront, thereby to obtain a transformed wavefront; and a phase change applier, applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 224. Apparatus according to claim 223, and wherein: said transform applied to said spectral analysis wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and said intensity map utilizer includes: a wavefront expresser, expressing said spectral analysis wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said spectral analysis wavefront; a first intensity map expresser, expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser, expressing each of said plurality of intensity maps as a third function of: said amplitude of said spectral analysis wavefront; said absolute value of said second complex function; a difference between said phase of said spectral analysis wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; a first function solver, solving said third function to obtain said amplitude of said spectral analysis wavefront, said absolute value of said second complex function and said difference between said phase of said spectral analysis wavefront and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said spectral analysis wavefront by adding said phase of said second complex function to said difference between said phase of said spectral analysis wavefront and said phase of said second complex function.
 225. Apparatus according to claim 218 and wherein: said spectral analysis wavefront comprises a plurality of different wavelength components; and said wavefront transformer is operative to apply a phase change to said plurality of different wavelength components of said spectral analysis wavefront.
 226. Apparatus for phase change analysis comprising: a wavefront obtainer, operative to obtain a phase change analysis wavefront which has an amplitude and a phase; a transform applier, applying a transform to said phase change analysis wavefront thereby to obtain a transformed wavefront; a phase change applier, applying a plurality of different phase changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indication of differences between said plurality of different phase changes applied to said transformed phase change analysis wavefront.
 227. Apparatus according to claim 226, and wherein when lateral shifts appear in said plurality of different phase changes, corresponding changes appear in said plurality of intensity maps, said intensity map utilizer obtains an indication of said lateral shifts.
 228. Apparatus according to claim 226, and wherein said intensity map utilizer includes: an intensity map expresser, expressing said plurality of intensity maps as at least one mathematical function of phase and amplitude of said phase change analysis wavefront and of said plurality of different phase changes, where at least said phase and amplitude are known and said plurality of different phase changes are unknown; and a function solver, employing said at least one mathematical function to obtain an output indicating said differences between said plurality of different phase changes.
 229. Apparatus for phase change analysis comprising: a wavefront obtainer, operative to obtain a phase change analysis wavefront which has an amplitude and a phase; a transform applier, applying a transform to said phase change analysis wavefront thereby to obtain a transformed wavefront; a phase change applier, applying at least one phase change to said transformed wavefront, thereby to obtain at least one phase changed transformed wavefront. an intensity map generator operative to obtain at least one intensity map of said at least one phase changed transformed wavefront; and an intensity map utilizer, employing said at least one intensity map to obtain an output indication of said at least one phase change applied to said transformed phase change analysis wavefront.
 230. Apparatus according to claim 229, and wherein said at least one phase change is a phase delay, having a value selected from a plurality of pre-determined values, and said output indication of said at least one phase change includes said value of said phase delay.
 231. Apparatus for stored data retrieval comprising: a wavefront obtainer operative to obtain a stored data retrieval wavefront which has an amplitude and a phase, by reflecting radiation from media in which information is encoded by selecting the height of the media at each of a multiplicity of different locations on the media; a wavefront analyzer, analyzing said stored data retrieval wavefront and comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to said stored data retrieval wavefront; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an indication of said amplitude and phase of said stored data retrieval wavefront; and a phase and amplitude utilizer, employing said indication of said amplitude and phase to obtain said information.
 232. Apparatus according to claim 231, and wherein said wavefront transformer comprises: a transform applier, applying a transform to said stored data retrieval wavefront thereby to obtain a transformed wavefront; and a phase change applier, applying a plurality of different phase changes to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 233. Apparatus according to claim 232, and wherein: said transform applied to said stored data retrieval wavefront is a Fourier transform; said plurality of different phase changes comprises at least three different phase changes, effected by applying a spatially uniform, time-varying spatial phase change to part of said transformed wavefront; said plurality of intensity maps comprises at least three intensity maps; and said intensity map utilizer includes: a wavefront expresser, expressing said stored data retrieval wavefront as a first complex function which has an amplitude and phase identical to said amplitude and phase of said stored data retrieval wavefront; a first intensity map expresser, expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser expressing each of said plurality of intensity maps as a third function of: said amplitude of said stored data retrieval wavefront; said absolute value of said second complex function; a difference between said phase of said stored data retrieval wavefront and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; a first function solver, solving said third function to obtain said amplitude of said stored data retrieval wavefront, said absolute value of said second complex function and said difference between said phase of said stored data retrieval wavefront and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said stored data retrieval wavefront by adding said phase of said second complex function to said difference between said phase of said stored data retrieval wavefront and said phase of said second complex function.
 234. Apparatus according to claim 232 and wherein: said stored data retrieval wavefront comprises at least one one-dimensional component; said transform applier is operative to perform a one-dimensional Fourier transform to said data retrieval wavefront, performed in a dimension perpendicular to a direction of propagation of said data retrieval wavefront, thereby to obtain at least one one-dimensional component of said transformed wavefront in said dimension perpendicular to said direction of propagation; said phase change applier is operative to apply a plurality of different phase changes to each of said at least one one-dimensional component, thereby obtaining at least one one-dimensional component of said plurality of phase changed transformed wavefronts; and said intensity map utilizer is operative to obtain an output indicating amplitude and phase of said at least one one-dimensional component of said data retrieval wavefront.
 235. Apparatus according to claim 234 and wherein said phase change applier comprises a movement generator, providing a relative movement between said media and a component generating spatially varying, time-constant phase changes, said relative movement being in a dimension perpendicular to said direction of propagation and to said dimension perpendicular to said direction of propagation.
 236. Apparatus according to claim 231, and wherein said information is encoded on said media whereby: an intensity value is realized by reflection of light from each location on said media to lie within a predetermined range of values, said range corresponding an element of said information stored at said location; and said intensity map utilizer employs said plurality of intensity maps to realize multiple intensity values for each location, providing multiple elements of information for each location on said media.
 237. Apparatus according to claim 236 and wherein said plurality of differently phase changed transformed wavefronts comprise a plurality of wavefronts whose phase has been changed by a phase changer, operative to apply an at least time varying phase change function to said stored data retrieval wavefront.
 238. Apparatus according to claim 236, and wherein: said stored data retrieval wavefront comprises a plurality of different wavelength components; and said wavefront transformer is operative to apply at least one phase change to said plurality of different wavelength components of said stored data retrieval wavefront.
 239. Apparatus according to claim 231, and wherein: said radiation which is reflected from said media comprises a plurality of different wavelength components, resulting in said stored data retrieval wavefront comprising a plurality of different wavelength components; and said wavefront transformer is operative to apply a phase change to said plurality of different wavelength components of said stored data retrieval wavefront.
 240. Apparatus according to claim 231, and wherein: said information encoded by selecting the height of the media at each of a multiplicity of different locations on the media is also encoded by selecting the reflectivity of the media at each of a plurality of different locations on the media; and said phase and amplitude utilizer includes a phase utilizer, employing said indication of said phase to obtain said information encoded by selecting the height of the media and an amplitude utilizer, employing said indication of said amplitude to obtain said information encoded by selecting the reflectivity of the media.
 241. Apparatus for 3-dimensional imaging comprising: a wavefront obtainer operative to obtain a 3-dimensional imaging wavefront, which has an amplitude and a phase, by reflecting radiation from an object to be viewed; and a wavefront analyzer, analyzing said 3-dimensional imaging wavefront and comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to said 3-dimensional imaging wavefront; an intensity map generator operative to provide a plurality of intensity maps of said plurality of differently phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating said amplitude and phase of said 3-dimensional imaging wavefront.
 242. Apparatus according to claim 241, and wherein said radiation reflected from said object has a narrow band about a given wavelength, causing said phase of said 3-dimensional imaging wavefront to be proportional to geometrical variations in said object, said proportion being an inverse linear function of said wavelength.
 243. Apparatus according to claim 241 and wherein said wavefront transformer comprises: a transform applier, applying a transform to said 3-dimensional imaging wavefront, thereby to obtain a transformed wavefront; and a phase change applier, applying a plurality of different phase changes, including spatial phase changes, to said transformed wavefront, thereby to obtain a plurality of differently phase changed transformed wavefronts.
 244. Apparatus according to claim 241 and wherein: said 3-dimensional imaging wavefront comprises a plurality of different wavelength components; and said wavefront transformer includes: a transform applier, applying a transform to said 3-dimensional imaging wavefront, thereby to obtain a transformed wavefront comprising a plurality of different wavelength components; and a phase change applier, applying phase changes to said plurality of different wavelength components of said transformed wavefront by passing said transformed wavefront through an object, at least one of whose thickness and refractive index varies spatially.
 245. Apparatus for wavefront analysis comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating at least phase of said wavefront being analyzed, comprising: an intensity combiner operative to combine said plurality of intensity maps into a second plurality of combined intensity maps, said second plurality being less than said first plurality; an indication provider operative to provide at least an output indicative of said phase of said wavefront being analyzed from each of said second plurality of combined intensity maps; and an enhanced indication provider, combining said outputs to provide at least an enhanced indication of phase of said wavefront being analyzed.
 246. Apparatus for wavefront analysis comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating at least amplitude of said wavefront being analyzed, comprising: an intensity combiner operative to combine said plurality of intensity maps into a second plurality of combined intensity maps, said second plurality being less than said first plurality; an indication provider operative to provide at least an output indicative of said amplitude of said wavefront being analyzed from each of said second plurality of combined intensity maps; and an enhanced indication provider, combining said outputs to provide at least an enhanced indication of amplitude of said wavefront being analyzed.
 247. Apparatus for wavefront analysis comprising: a wavefront transformer operative to provide a plurality of differently phase changed transformed wavefronts corresponding to a wavefront being analyzed; an intensity map generator operative to provide a plurality of intensity maps of said plurality of phase changed transformed wavefronts; and an intensity map utilizer, employing said plurality of intensity maps to provide an output indicating at least phase of said wavefront being analyzed, comprising: an intensity map expresser, expressing said plurality of intensity maps as a function of: amplitude of said wavefront being analyzed; phase of said wavefront being analyzed; and a phase change function characterizing said plurality of differently phase changed transformed wavefronts; a complex function definer, defining a complex function of: said amplitude of said wavefront being analyzed; said phase of said wavefront being analyzed; and said phase change function characterizing said plurality of differently phase changed transformed wavefronts, said complex function being characterized in that intensity at each location in said plurality of intensity maps is a function predominantly of a value of said complex function at said location and of said amplitude and said phase of said wavefront being analyzed at said location; a complex function expresser, expressing said complex function as a function of said plurality of intensity maps; and a phase obtainer, obtaining values for said phase by employing said complex function expressed as a function of said plurality of intensity maps.
 248. Apparatus for wavefront analysis comprising: a first transform applier, applying a Fourier transform to a wavefront being analyzed which has an amplitude and a phase thereby to obtain a transformed wavefront; a phase change applier, applying a spatially uniform time-varying spatial phase change to part of said transformed wavefront, thereby to obtain at least three differently phase changed transformed wavefronts; a second transform applier, applying a second Fourier transform to said at least three differently phase changed transformed wavefronts, thereby obtaining at least three intensity maps; and an intensity map utilizer, employing said at least three intensity maps to provide an output indicating at least one of said phase and said amplitude of said wavefront being analyzed and comprising: a wavefront expresser, expressing said wavefront being analyzed as a first complex function which has an amplitude and phase identical to said amplitude and phase of said wavefront being analyzed; a first intensity map expresser. expressing said plurality of intensity maps as a function of said first complex function and of a spatial function governing said spatially uniform, time-varying spatial phase change; a complex function definer, defining a second complex function having an absolute value and a phase as a convolution of said first complex function and of a Fourier transform of said spatial function governing said spatially uniform, time-varying spatial phase change; a second intensity map expresser, expressing each of said plurality of intensity maps as a third function of: said amplitude of said wavefront being analyzed; said absolute value of said second complex function; a difference between said phase of said wavefront being analyzed and said phase of said second complex function; and a known phase delay produced by one of said at least three different phase changes, which each correspond to one of said at least three intensity maps; a first function solver, solving said third function to obtain said amplitude of said wavefront being analyzed, said absolute value of said second complex function and said difference between said phase of said wavefront being analyzed and said phase of said second complex function; a second function solver, solving said second complex function to obtain said phase of said second complex function; and a phase obtainer, obtaining said phase of said wavefront being analyzed by adding said phase of said second complex function to said difference between said phase of said wavefront being analyzed and said phase of said second complex function. 